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Abstract 



Within the framework of lattice QCD we investigate different aspects of 
QCD in Landau gauge using Monte Carlo simulations. In particular, we 
focus on the low momentum behavior of gluon and ghost propagators. The 
gauge group is that of QCD, namely SU{3). For our study of the lattice 
gluodynamic, simulations were performed on several lattice sizes ranging from 
12^ to 48"^ at the three values of the inverse coupling constant f3 = 5.8, 6.0 
and 6.2. 

Different systematic effects on the gluon and ghost propagators are stud- 
ied. We demonstrate that the ghost dressing function systematically depends 
on the choice of Gribov copies at low momentum, while the influence on the 
gluon dressing function is not resolvable. Also the eigenvalue distribution of 
the Faddeev-Popov operator is sensitive to Gribov copies. 

We show that the influence of dynamical Wilson fermions on the ghost 
propagator is negligible at the momenta available to us. For this we have 
used gauge configurations which were generated with two dynamical flavors 
of clover-improved W'ilson fermions. On the contrary, fermions affect the 
gluon propagator at large and intermediate momenta, in particular where 
the gluon propagator exposes its characteristic enhancement compared to 
the free propagator. 

We also analyze data for both propagators obtained on asymmetric lat- 
tices. By comparing these results with data obtained on symmetric lattices, 
we flnd that both the gluon and the ghost propagator suffer from systematic 
effects at the lowest on-axis momenta available on asymmetric lattices. 

We compare our data with the infrared exponents predicted in studies 
of truncated systems of Dyson-Schwinger equations for the gluon and ghost 
propagators. We cannot conflrm neither the values for both exponents nor 
the relation which is proposed to hold between them. In any case, we demon- 
strate that the infrared behavior of gluon and ghost propagators, as found 
in this thesis, is consistent with different criteria for conflnement. In fact, 
we verify that our data of the ghost propagator and also of the Kugo-Ojima 
conflnement parameter satisfy the Kugo-Ojima conflnement criterion. The 
Gribov-Zwanziger horizon condition is satisfled by the ghost propagator. Also 
the gluon propagator seems to vanish in the zero-momentum limit. However, 
we cannot judge without doubt on the existence of an infrared vanishing 
gluon propagator. Furthermore, explicit violation of reflection positivity by 
the transverse gluon propagator is shown for the quenched and unquenched 
case of SU (3) gauge theory. 



The running coupling constant given as a renormalization-group-invariant 
combination of the gluon and ghost dressing functions does not expose a fi- 
nite infrared fixed point. Rather the data are in favor of an infrared vanishing 
coupUng constant. This behavior does not change if the Gribov ambiguity 
or unquenching effects are taken into account. We also report on a first non- 
perturbative computation of the SU(3) ghost-gluon- vertex renormalization 
constant. We find that it deviates only weakly from being constant in the 
momentum subtraction scheme considered here. 

We present results of an investigation of the spectral properties of the 
Faddeev-Popov operator at /3 = 5.8 and 6.2 using the lattice sizes 12"', 16^ 
and 24"^. For this we have calculated the low-lying eigenvalues and eigen- 
modes of the Faddeev-Popov operator. The larger the volume the more 
eigenvalues are found accumulated close to zero. Using the eigenmodes for 
a spectral representation of the ghost propagator it turns out that for our 
smallest lattice only 200 eigenvalues and eigenmodes are sufficient to saturate 
the ghost propagator at lowest momentum. We associate exceptionally large 
values occurring occasionally in the Monte Carlo history of the ghost propa- 
gator at larger /3 to extraordinary contributions of the low-lying eigenmodes. 



Keywords: 

Gluon and ghost propagators, lattice QCD, Landau gauge, confinement 



iv 



Zusammenfassung 



Diese Arbeit untersucht im Rahmen der Gittereichtheorie verschiedene As- 
pekte der QCD in der Landau-Eichung, insbesondere solche, die mit den 
Gluon- und Geist-Propagatoren zusammenhangen. Die Eichgruppe ist die 
der QCD, SU{3), und wir untersuchen die Propagatoren bei kleinen Im- 
pulsen. Fiir unsere Untersuchungen der reinen Gluodynamik haben wir 
zahlreiche Monte-Carlo Simulationen auf diversen Gittergrossen durchge- 
fiihrt. Die Gittergrossen variieren im Bereich von 12^ bis 48^. Als inverse 
Kopplungskonstanten haben wir die Werte (3 = 5.8, 6.0 und 6.2 gewahlt. 

Wir analysieren den Einfluss unterschiedlicher systematischer Effekte auf 
das Niedrigimpulsverhalten der Gluon- und Geist-Propagatoren. Wir zeigen, 
dass der Formfaktor des Geist-Propagators bei kleinen Impulsen systematisch 
von der Wahl der Eichkopien (Gribov-Kopien) abhangt. Hingegen konnen wir 
einen solchen Einfluss auf den Gluon-Propagator nicht feststellen. Ebenfalls 
wird die Verteilung der kleinsten Eigenwerte des Faddeev-Popov-Operators 
durch die Wahl der Gribov-Kopien beeinflusst. 

Wir zeigen aufierdem, dass der Einfluss dynamischer Wilson-Fermionen 
auf den Geist-Propagator fiir die untersuchten Impulse vernachlassigbar ist. 
Dazu haben wir Eichkonflgurationen betrachtet, die mit einer Nf — 2 clover- 
verbesserten Wirkung erzeugt worden sind. Fiir den Gluon-Propagator kon- 
nen wir jedoch einen deutlichen Einfluss fiir gro£ie und mittlere Impulse fest- 
stellen, insbesondere in dem Impulsbereich, wo der Gluon-Propagator im 
Vergleich zum freien Fall seine charakteristische Erhohung aufweist. 

Zusatzlich wurden beide Propagatoren auf asymmetrischen Gittern gemessen. 
Der Vergleich dieser Daten mit denen, die auf symmetrischen Gittern gewon- 
nen wurden, zeigt, dass die Asymmetrie deutliche systematische Effekte im 
Bereich kleiner Impulse verursacht. Besonders deutlich wird das fiir die Dat- 
en, die bei Impulsen in Richtung der elongierten Gitterlange gemessen worden 
sind. 

Weiterhin vergleichen wir unsere Daten mit den Infrarot-Exponenten, die 
in Studien von abgeschnittenen (tnmcated) Systemen von Dyson-Schwinger- 
Gleichungen fiir den Gluon- und Geist-Propagator vorhergesagt wurden. Im 
Rahmen unserer Messungen konnen wir weder die Werte der Exponenten 
noch die vorhergesagte Beziehung zwischen beiden bestatigen. In jedem 
Falle konnen wir aber zeigen, dass das in dieser Arbeit gefundene Nied- 
rigimpulsverhalten im Einklang mit verschiedenen Kriterien fiir ConJi,nement 
(Einschluss von Farbladungen) ist. Wir zeigen, dass unsere Daten sowohl fiir 
den Geist-Propagator als auch fiir den Kugo-Ojima-Conflnement-Parameter 



das Kugo-Ojima-Confinement-Kriterium erfiillen. Aufterdem ist die Gribov- 
Zwanziger-Horizontbedingung fiir den Geist-Propagator erfullt. Der Gluon- 
Propagator scheint im Grenzfall verschwindender Impulse zu Null zu streben. 
Dennoch konnen wir nicht endgiiltig dariiber urteilen, ob dies der Fall ist. 
Wir zeigen zusatzlich, dass der transversale Gluon-Propagator explizit die 
Reflektions-Positivitat verletzt. Das gilt sowohl mit als auch ohne den Ein- 
fluss dynamischer Fermionen. 

Wir berechnen die laufende (efFektive) Kopplung, die sich als 
eine renormierungsgruppeninvariante Kombination der Gluon- und Geist- 
Formfaktoren ergibt. Unsere Ergebnisse zeigen deutlich, dass im Bere- 
ich kleiner Impulse die laufende Kopplung kleiner wird und so vermut- 
lich kein endlicher Infrarot-Fixpunkt im Grenzfall Impuls Null angestrebt 
wird. Dieses Verhalten ist unabhangig vom Einfluss der Gribov-Kopien 
oder von der Hinzunahme dynamischer Fermionen. Wir prasentieren aufser- 
dem eine erste nichtstorungstheoretische Berechnung der Renormierungskon- 
stante des SU{3) Ghost-Gluon- Vertex. Wir zeigen, dass in dem untersuchten 
Renormierungsschema keine wesentliche Abweichung von einem konstanten 
Verhalten gefunden wird. 

Wir berichten aufserdem iiber Untersuchungen zu spektralen Eigen- 
schaften des Faddeev-Popov-Operators bei (3 = 5.8 and 6.2. Dazu haben 
wir eine Reihe der kleinsten Eigenwerte und Eigenvektoren dieses Opera- 
tors auf den Gittergroften 12^, 16^ und 24^ berechnet. Wir sehen, dass sich 
umso melir Eigenwerte nahe Null konzentrieren, je grower das physikalische 
Volumen ist. Anhand einer spektralen Entwicklung des Geist-Propagators 
konnen wir zeigen, dass fiir unser kleinstes Gitter ca. 200 Eigenwerte und 
Eigenvektoren geniigen, um den Wert des Geist-Propagators beim kleinsten 
Impuls zu reproduzieren. 

Wir zeigen ferner, dass die selten auftretenden, exzeptionell grol?,en Mess- 
werte, die fiir den Geist-Propagator im Verlauf der Monte-Carlo Simula- 
tion bei grofseren /3 Werten gefunden werden, durch aufierordentlich starke 
Beitrage der niedrigsten Eigenmoden zu den entsprechenden Fourierkompo- 
nenten hervorgerufen werden. 



Schlagworter: 

Gluon- und Geist-Propagatoren, Gitter-QCD, Landau-Eichung, 
Confinement 



vi 



Contents 



Abstract iii 

Zusammenfassung (german) v 

Introduction 1 

1 The various colors of QCD 7 

1.1 Quantization of QCD 7 

1.1.1 The classical QCD Lagrangian 7 

1.1.2 Functional-integral quantization of QCD 10 

1.1.3 The Faddeev-Popov method 12 

1.1.4 An effective Lagrangian density in covariant gauge ... 13 

1.1.5 The BRST formalism 16 

1.2 Regularization and renormalization 19 

1.2.1 Regularization 19 

1.2.2 Renormalization 20 

1.2.3 The MOM scheme 21 

1.3 The renormalization group 23 

1.3.1 The renormalization group equation 23 

1.3.2 Perturbative expansion of the /3-function 25 

1.3.3 The running coupling constant 26 

1.3.4 The anomalous dimension 28 

2 Infrared QCD and criteria for confinement 31 

2.1 Nonperturbative approaches to QCD 31 

2.1.1 Brief remarks on nonperturbative methods 32 

2.1.2 The problem of Gribov copies 32 

2.1.3 Nonperturbative quantization in Landau gauge .... 33 

2.2 The Dyson Schwinger equations of QCD 35 

2.2.1 Infrared behavior of ghost and gluon propagators in 
Landau gauge 36 

2.2.2 A nonperturbative running coupling constant 39 

2.2.3 The finiteness of the ghost-gluon vertex 40 

2.3 Criteria for confinement in linear covariant gauges 41 

2.3.1 The Kugo-Ojima confinement scenario 42 

2.3.2 The Gribov-Zwanziger horizon condition 48 



vii 



VIII 



CONTENTS 



2.3.3 Violation of reflection positivity as a criterion for con- 

flnement 49 

3 QCD Green's functions in lattice Landau gauge 53 

3.1 Basics of lattice QCD 53 

3.1.1 The lattice and its flelds 54 

3.1.2 The Wilson action with clover-improved fermions ... 55 

3.1.3 Vacuum expectation values from MC simulations ... 57 

3.2 The Landau gauge on the lattice 58 

3.2.1 The gauge functional 59 

3.2.2 The Faddeev-Popov operator 60 

3.2.3 Defining F, Q and A on the lattice 61 

3.3 Lattice definition of our observables 62 

3.3.1 The (inverse) FP operator in momentum space 62 

3.3.2 The gluon and ghost propagator 64 

3.3.3 The Kugo-Ojima confinement parameter 67 

3.3.4 Renormalization of propagators 68 

3.3.5 The ghost-gluon-vertex renormalization constant .... 68 

4 Results for lattice QCD Green's functions 71 

4.1 General prerequisites 71 

4.1.1 Specification of our lattice samples 71 

4.1.2 Gauge fixing 72 

4.1.3 The f c-bc strategy 72 

4.1.4 Selection of momenta 74 

4.1.5 Mapping to physical units 75 

4.2 Systematic effects on gluon and ghost propagators at low mo- 
mentum 75 

4.2.1 Finite volume and discretization effects 75 

4.2.2 Asymmetric lattices cause strong systematic errors . . 77 

4.2.3 Facing the problem of Gribov copies 80 

4.3 The infrared behavior of gluon and ghost propagator in the 
quenched and unquenched case 85 

4.3.1 The gluon propagator 85 

4.3.2 The ghost propagator 88 

4.4 The running coupling and the ghost-gluon vertex 92 

4.4.1 Results for the running coupling constant 92 

4.4.2 The vertex renormalization constant 95 



CONTENTS IX 



5 Confinement criteria under the lattice microscope 99 

5.1 Is the Gribov-Zwanziger horizon condition satisfied? 99 

5.2 The Kugo-Ojima confinement parameter 102 

5.2.1 Expected infrared behavior 102 

5.2.2 Explicit lattice data for the function ^(g^) 103 

5.3 The gluon propagator explicitly violates refiection positivity . 107 

6 Spectral properties of the FP operator 111 

6.1 Specification of lattice samples Ill 

6.2 The low-lying eigenvalue spectrum 112 

6.2.1 The lowest and second lowest eigenvalues 112 

6.2.2 An estimate for the density of low-lying eigenvalues . . 115 

6.3 Eigenmode expansion of the ghost propagator 116 

6.4 The problem of exceptional configurations 117 

6.5 Localization properties of low-lying eigenvectors 120 

Conclusions and outlook 123 

A Some details on algorithms and performance 131 

A.l A note on the algorithms used 131 

A. 2 Experience report on lattice Landau gauge fixing 132 

A. 2.1 Over-relaxation versus Fourier-accelerated gauge-fixing 132 

A. 2. 2 A way to preselect best gauge copies 134 

A. 3 Speeding up the inversion of the FP operator 136 

Bibliography 139 

Acknowledgments 161 



X CONTENTS 



Introduction 



At present we are reasonably confident that the physics of strong in- 
teraction, i.e. the rich field of hadron physics, is completely described 
by a quantized nonabelian gauge field theory which is based on the gauge 
group of SU{3) color symmetry. This theory is called Quantum Chromody- 
namics (QCD). Its fundamental constituents are quarks and gluons. Quarks 
are spin 1/2 fermion fields carrying fractional electric charge and the gluons 
are nonabelian spin 1 gauge fields which interact with the quarks as well as 
among themselves. Due to its nonabelian nature the renormalization group 
tells us that QCD is asymptotically free at large Euclidean momentum. In 
this regime perturbative QCD is relevant and theoretical predictions have 
been successfully confronted with experiments. The experimental successes 
of QCD and the partial progress towards a full understanding of the theory 
form the basis for our present belief that QCD is the right theory describing 
all strong interaction physics. 

Beyond perturbation theory, however, QCD is still not completely un- 
derstood, even though — as far as we know — it is not in conflict with 
any existing phenomenology of the strong interaction. Note that in contrast 
to QED the elementary fields in QCD, the quarks and gluons, do not de- 
scribe existing particles and thus a particle interpretation in QCD has to 
be completely divorced from its elementary degrees of freedom. According 
to QCD all strongly interacting particles, the hadrons, are colorless bound 
states of quarks. This phenomenon is called confinement, but the mechanism 
which confines quarks and gluons has to be established yet from first princi- 
ples. Moreover, due to the complexity of QCD, a full description of hadronic 
states and processes directly in terms of QCD presents an exciting challenge 
since many years. 

Many hadronic features have been investigated in the framework of phe- 
nomenological models (see e.g. [VW91; Kle92; ERV94]) which mimic the es- 
sential properties of QCD, namely asymptotic freedom at short distance and 
confinement at large distances. This approach represents a rather practical 
point of view and is sufficient if one is just interested in the effective theory 
of hadrons at low energies. But if QCD is the theory of strong interactions a 
coherent description directly based on the dynamics of confined quarks and 
gluons should be possible. 
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Introduction 



For such a description a complete picture for all propagators and vertex 
functions of QCD should be available. These Green's functions may then 
serve as input into bound state calculations based on the Bethe-Salpeter 
equations for mesons or the Faddeev equations for baryons. But also from 
a purely theoretical point of view a consistent picture of all QCD Green's 
functions is interesting. In particular, their infrared momentum behavior pro- 
vides insight into the mechanism of quark and gluon confinement [AvSOl]. To 
give just one example: The realization of the Kugo-Ojima confinement sce- 
nario [K079; Kug9o] in QCD in covariant gauges is encoded in the infrared 
behavior of the ghost 2-point function. Therefore, the investigation of QCD 
Green's function at low momentum is important for a coherent description of 
hadronic states and processes and also for an understanding of confinement. 

The infrared momentum region corresponds to strong coupling rather 
than weak coupling and hence perturbation theory is of no avail in studying 
QCD at low momentum. Genuinely nonperturbative approaches have to be 
used to explore QCD in this area. The Euclidean space, discretized, lattice 
gauge theory provides one possibility to study nonperturbative aspects of 
QCD by using Monte Carlo (MC) simulations. Another approach is given 
by solving trimcated systems of the Dyson-Schwinger equations (DSEs) of 
QCD. The DSEs are infinite towers of coupled nonlinear integral equations 
relating different Green's functions of QCD to each other. They are directly 
derived from a generating functional whose existence beyond perturbation 
theory still has to be assumed. In any case, studying DSEs involves the 
introduction of a gauge condition which is not necessary in the standard 
lattice approach to QCD. 

DSE studies have been performed in recent years with growing intensity 
(see [RW94; RSOO; AvSOl; MR03] for an overview). In particular, for the case 
of Landau gauge it has been shown [vSAH97; vSHA98] that contributions of 
ghost fields are crucial for a consistent description of the infrared behavior 
of Landau gauge giuodynamics. In former studies [Man79; ADJS81; AJS82; 
BP89], ghost fields have always been neglected. 

Different truncations have been employed since then to study the in- 
frared behavior of gluon, ghost and quark propagators and the correspond- 
ing vertex functions. Truncations are essential to manage the infinite towers 
of DSEs. The solutions presented first in [vSAH97; vSHA98] and later in 
[AB98a; AB98b; BloOl; Blo02] and [FAR02; FA02; Fis03] all favor the pic- 
ture of an infrared diverging ghost propagator being intimately connected 
with an infrared vanishing gluon propagator. In fact, both propagators are 
proposed to follow power laws at low momentum with intertwined infrared 
exponents [LvS02; Zwa02]. Such an infrared behavior is in agreement with 
the Gribov-Zwanziger horizon condition [Gri78; Zwa94; Zwa02; Zwa04] as 
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well as with the Kugo-Ojima confinement criterion [K079; Kug95]. Note 
that their satisfaction is crucial for the realization of confinement in QCD 
in Landau gauge. Unquenching effects on the infrared behavior are found 
to be small [FA03]. Moreover, dynamical chiral symmetry breaking and 
gluon confinement have been confirmed from solutions of truncated DSEs 
[AvSOl; FA03; ADFM04]. 

Most of these DSE studies are done in Landau gauge. In this gauge, 
the ghost-gluon vertex was shown to not suffer from ultraviolet divergences 
at any order in perturbation theory [Tay71; MP 78]. Assuming this to hold 
beyond perturbation theory, it allows for a definition of a nonperturbative 
running coupling constant that is solely given in terms of the gluon and ghost 
propagators and has a finite infrared fixed point, provided the mentioned in- 
frared power laws hold [vSAH97; vSHA98]. In a recent DSE study [AFLE05] 
of vertex functions, the infrared fixed point has been confirmed, too. It has 
also been shown that this coupling constant enters directly the kernels of the 
DSEs for the gluon, ghost and quark propagators [BloOl; Blo02]. 

Even though Monte Carlo simulations of lattice QCD provide an alter- 
native possibility to study QCD at a nonperturbative level, at present, they 
cannot compete with the DSE approach concerning the accessible region of 
low momenta. However, lattice QCD is a first principle approach to QCD 
that does not require us to "simplify" the theory. Unlike truncations of DSEs, 
the approximations involved in lattice QCD are systematically removable. 
This possibility of controlling the systematic errors makes this approach in- 
valuable [BHL+Oo]. Therefore, lattice simulations may provide an indepen- 
dent check whether the results obtained in the DSE approach are realized 
in lattice QCD, at least in the region of momenta available at present. Fur- 
thermore, lattice QCD enables us to study different models for confinement 
(see e.g. [Gre03]) by mutilating the theory such that confinement is explic- 
itly lost. For example, removing vortices changes the infrared behavior of 
the lattice ghost propagator in Landau gauge such that it does not satisfy 
anymore aforementioned criteria for confinement [GLR04]. 

In recent years, different groups have investigated different aspects of lat- 
tice Landau gauge QCD. Some have studied the gauge group SU{2), others 
SU{3). In particular, the Adelaide Group has provided an impressive ac- 
count on numerical data for the SU{3) gluon [LSWP98; LSWP99; BBLWOO; 
BBL+01] and quark propagators [SWOl; SLWOl; BBL+02; BHW02; ZBL+04; 
Z+05; B+05d; P+06] and for the quark-gluon vertex [SK02; SBK+03]. Their 
data are based on quenched and unquenched SU{3) gauge configurations 
where the latter were generated with the AsqTad quark action by the MILC 
collaboration. 
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For the SU{3) ghost propagator there were not so many data available 
until a few years ago, even though this propagator is expected to be related to 
the gluon propagator as mentioned above. The first lattice study of the SU (2) 
and SU (3) ghost propagators in Landau gauge was given in [SS96] and there 
were several studies in Landau gauge which have confirmed the anticipated 
behavior for the case of ,5[/(2) [BCLM04; LRG02; GLR04; BCLM03]. Similar 
investigations for the SU{3) case at even lower momenta were not available 
at that time. 

During the last three years, we have tried to bridge this gap by investigat- 
ing the SU (3) ghost and gluon propagators (and related objects) on quenched 
and unquenched SU (3) gauge configurations. Our set of unquenched configu- 
rations were generated with clover-improved Wilson fermions by the QCDSF 
collaboration. We have found that for momenta lower than used in the SU{2) 
studies (see above) qualitative differences to the anticipated infrared behavior 
of ghost and gluon propagators and of the running coupling constant appear. 

At the same time other groups have performed similar investigations fo- 
cussing on different interesting aspects. See, for example, [FN04a; FN04b; 
FN05a; FN06b; FN06a] for investigations of the gluon, ghost and quark 
propagators and of the running coupling constant using quenched and un- 
quenched configurations (provided by the MILC collaboration). Studies of 
the SU (3) gluon propagator at very low momentum can be found in [OS05b; 
SO05b; OS05a; SO05a]. There the infrared exponent has been determined 
using lattices much elongated in time direction. A study of the SU{3) ghost 
propagator at large momentum can be found in [B^Ooc]. 

Furthermore, recent DSE studies [FAR02; FA02; FGA06; FP06] show that 
the infrared behavior of the gluon and ghost dressing functions and of the 
running coupling constant is changed on a torus. In particular, the running 
coupling decreases at low momenta. These findings agree with lattice data as 
shown in this thesis, but they contradict results obtained from DSE studies in 
the continuum. It is still unknown what is the reason for this disagreement. 
Note that a solution to this problem has been proposed in [B+OSb; B"'"06a]. 

It is the intention of this thesis to give a summary of our results obtained 
within the last three years. Some were already published, others are just 
finished and being written up. 

We have structured this thesis as follows: In the first chapter we intro- 
duce the path integral formulation of QCD and discuss the special problems 
related to the necessity of gauge fixing the action. A brief introduction to 
the BRST formalism is given and the renormalization program is recalled. 
In Chapt. 2 we discuss some aspects of nonperturbative QCD and intro- 
duce criteria for confinement which are available for QCD in Landau gauge. 
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The lattice formulation of QCD in this gauge and a definition of all observ- 
ables analyzed in this thesis is given in Chapt. 3. We present our results for 
the ghost and gluon propagators in Chapt. 4. Different systematic effects 
are analyzed and their influence on the infrared behavior of the giuon and 
ghost propagators is discussed. After this, we show that our data for both 
propagators satisfy necessary criteria for confinement. In Chapt. 6 spectral 
properties of the FP operator are analyzed. Finally, we draw our conclusions 
and give an outlook. The appendix contains some notes on algorithms and 
performance. In particular, we compare two popular gauge-fixing algorithms 
and show that the final ranking of gauge functional values is already visible 
at an intermediate iteration state. We also demonstrate how the inversion of 
the FP operator can be accelerated considerably. 
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Chapter 1 



The various colors of QCD 



This chapter briefly reviews the EucHdean formulation of QCD in the con- 
tinuum, mainly in order to fix notations used subsequently. Starting with 
the classical Lagrangian density and its quantization, the problems encountered by 
fixing to covariant gauges are discussed. A short summary of the BRST formalism 
and renormalization is given. 



1.1 Quantization of QCD 

The success of quark-models in describing hadrons as bound states of quarks, 
but also of quark-parton models in deep-inelastic lepton-hadron scattering, to 
name but a few, suggests that the strong interaction should be described by 
a theory v^^here the color symmetry of each quark flavor is a gauge symmetry 
and which is also asymptotically free at high energy-momentum transfers or 
short distances. Since asymptotic freedom is inherent in non-abelian gauge 
theories, and experiments like, for instances, the pion-decay 7r° — > 27 suggest 
the gauge group to be SU{3), we are reasonably confident at present that the 
strong interaction is completely described by a quantized non-Abelian gauge 
field theory based on the SU (3) gauge group. 

1.1.1 The classical QCD Lagrangian 

A common way to setup a quantum field theory is to define first a Lagrangian 
density 

£ = . . . , $„(a;); 9^<l>i(a;), . . . , d^<^n{x)] (1.1) 

that is a functional of several fields . . . , $„(a;) and their derivatives 

necessary to host (in a consistent way) all the features and symmetries ob- 
served in experiments. This Lagrangian density or its space-time integral, 
the action 

Sm = [ d^xC, (1.2) 
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Chapter 1 The various colors of QCD 



is then used subsequently for a quantization of the theory choosing one of the 
well-known quantization methods, namely the Canonical operator formalism, 
the Stochastic formalism or the Functional-integral formalism [Mut98]. 

The most general form of the QCD Lagrangian density that not only 
accommodates all those mentioned properties of QCD, but also is renormal- 
izable in any order of perturbation theory can be written (in Euclidean space) 
as ^ 

Anv = If;^F^''''^ - i^i^^D, - mo)V' . (1.3) 

Conceived in general terms, this Lagrangian density describes the interaction 
of the quark and antiquark fields, ip and ip, with the self-interacting gluon or 
gauge fields ^4^ = A^(a;)T". The latter are hidden in both the definition of 
the field-strength tensor 

F^. = d,Al - d^Al - gof^'^A'^At (1.4) 

(here in the adjoint representation, i.e. a = 1, . . . , A^"^ — 1) and the covariant 
derivative 

D'J ^ d.S"' + igoA'^^iT'^f (1.5) 

given in the fundamental representation (i.e. k,l = 1, . . . , Nc) of the Lie 
group SU{Nc — 3) with the eight hermitian generators T". Beside being 
hermitian, these generators satisfy Tr (T'^T^) = 5"V2 and [T^T^] = i/^^^T^ 
where /"^'^ are the structure constants of the Lie algebra su(3). The bare 
coupling constant is labeled gg- For sake of completeness, we also remind 
on the covariant derivative in the adjoint representation: 

< = + gor'^^M- (1-6) 

The quark fields 

and antiquark fields 'ip = ip'^Joj of flavor / = 1, . . . , Nf are anti-commuting 
spinor fields that transform under the fundamental representation of the 
SU{Nc = 3) color group, i.e. the color index runs over / = 1, . . . ,Nc. The 
Dirac matrices 7^ act upon the spinor indices a = 1, . . . , 4 of the quark fields. 
The bare mass mo is a free parameter (for each flavor) of the theory as is gg. 

^Here and in the following, a sum over repeated indices is understood if not otherwise 
stated. We will see later that there is always the freedom to have multiplicative renormal- 
ization constants or to add BRST-exact terms, like e.g. gauge-fixing and ghost terms in 
covariant gauges to this density. 
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By definition, the Lagrangian density in Eq. (1.3) is invariant under local 
SU{2>) gauge transformations 

— Qujd^gl, (1.7a) 
So 

(1.7b) 
(1.7c) 

of gluon, quark and antiquark fields. Here is an element of the group 
SU{2>). It can be parameterized by a set of real- valued functions uj^{x), i.e. 

= gUx) = e-^o--^-)^" G SU{3). (1.8) 

In subsequent discussions we will frequently refer to the infinitesimal form 
of those local transformations. What is usually meant by that notion is the 
following. If the field $fc = {tp, ip, A} transforms imder a local gauge transfor- 
mation $/; — > as given in Eq. (1.7) then the corresponding infinitesimal 
transformation is defined by [Col84]: 



OUJh 



uj=0 



Using Eq. (1.8) the infinitesimal local gauge transformations of the gluon and 
fermion fields take the form: 

5^Al = d,u'' + g,r''^uj'Al^Dfu' (1.9a) 

^u.i' = -tgouj^ry (1.9b) 

The invariance of the Lagrangian density Ci^^ under local gauge trans- 
formations, causes some extra difficulties for the quantization using either 
the functional-integral or the canonical formalism. For example, the defi- 
nition of a functional-integral over gauge fields in the continuum requires a 
gauge condition to be introduced. As a consequence additional terms are 
added to i2inv In the resulting Lagrangian density, £efr, the gauge invari- 
ance is explicitly lost, but its particular form — it is BRST invariant (see 
below) — guarantees that expectation values of gauge-invariant observables 
are actually independent of the gauge condition used. 

We note in passing that on the lattice such a gauge condition is super- 
fluous, as long as gauge-invariant observables are studied. Therefore, the 
gauge-invariant action 

5i„v= I d^xC,^, (1.10) 

is sufficient for a lattice discretization. See Chapt. 3 for a particular lattice 
discretization as used in this study. 
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1.1.2 Functional-integral quantization of QCD 

So far the theory is a classical field theory. To quantize it one chooses one of 
the well-known quantization methods, namely the Canonical operator formal- 
ism, the Stochastic formalism or the Functional-integral formalism. Indeed, 
all three methods should lead to the same physical predictions. However, the 
choice depends on the feasibility of the method for a particular topic. 

A quantum field theory is completely characterized by the infinite hierar- 
chy of n-point functions or Green's functions. These are correlation functions 
of the fields ^i{x) and the three mentioned formalisms differ in how Green's 
functions are calculated. For example, in the canonical approach the fields 
are regarded as operators for which canonical commutation relations hold. 
The Green's functions are calculated as vacuum expectation values of time 
ordered products of those operators. The stochastic formalism introduced by 
Parisi and Wu [PW81] starts from the classical equation of motion. The 
fields are regarded as stochastic variables. See [DH87] for a comprehensive 
account on that subject. 

The Functional-integral approach was introduced by Feynman [Fey48]. 
There the fields are taken to be c-numbers and the Lagrangian density takes 
its classical form. The Green's functions are given by functional integrations 
of products of fields over all of their (weighted) possible functional forms. 
The present study focuses on the lattice regularization of QCD in Euclidean 
space. Since this approach relies on the functional integral formalism we 
demonstrate briefiy the general concept.^ 

Functional-integral formalism: Illustration of the general concept 

The functional-integral formalism introduces generating functionals Z, W, 
and r which generate, respectively, the full, connected and one-particle irre- 
ducible (IPI) Green's functions. To get acquainted with the general concept 
let us assume that for the generic Lagrangian density (Eq. (1.1)) of n different 
fields $i the generating functional 



for the full Green's functions can be defined, i.e. there exist a well-defined 
measure [V^]. Then a full Green's function (^"^(xi) • • •<l>^"(x„)) is given by 

^Note that due to the work of KuGO and Ojima [K079] a consistent quantization of 
non-abehan gauge fields is also available in the covariant canonical operator formalism 
[Mut98]. Some of their results, namely the Kugo-Ojima confinement scenario will also be 
investigated in this study. For the covariant canonical operator formalism see also the 
book by Nakanishi and Ojima [NO90]. 




(1.11) 



1.1 Quantization of QCD 



11 



functional derivatives with respect to the sources Ji{x), i.e. 



j"i,...,J^"=0 



Jr(a;i)---J^(a;„ 

Together with Eq. (1.11) and the generic action S[^] (Eq. (1.2)) this yields 

($r {xi)---K"M) = ^J [m $r (^o ■ ■ ■ m e-^[*i . 

Gauge orbits and gauge conditions 

For a quantization of QCD within the functional formalism it is necessary to 
define the generating functional Z[J] that generates all the Green's functions 
of the theory. In particular, the definition of a path-integral over gluon fields 
needs special care, because it is ill-defined if done naively. 

In fact, choosing a particular gauge field \{x) there are infinitely many 
others '^A^ which are related to this by local gauge transformations as defined 
in Eq. (1.7a). The set of all those is usually referred to as the gauge orbit of 
^A^, because each element "^A^ in the orbit is obtained by acting upon 
with a local gauge transformation guj{x). The Lagrangian Ci^y is invariant 
under such a transformation by definition and so all (infinite) elements of 
one particular orbit give rise to the same value of C-mv This spoils a naive 
integration over all gluon fields, because an integral of kind 

J [DA] e-^'- = I [D^A] e-^'- J [Du] 

is divergent. Here Si^y denotes the action in Eq. (1.10), but for simplicity we 
have dropped fermionic fields. The integration over the gluon fields must be 
defined such that it restricts to gauge-inequivalent configurations, i.e. they 
must belong to different gauge orbits. 

This can be achieved by choosing a gauge condition 



J^[^A;x] 



(1.12) 



at each point x in space-time. If this condition is satisfied for only one 
representative on each gauge orbit, i.e. the solution lD is unique, then it is 
called an ideal gauge condition [Wil03]. The set of those representatives is 
called the fundamental modular region A. It is a hypersurface defined by 
Eq. (1.12) in the space of all gauge fields. If we can define an integration 
over this region, the integral 



/ [DA] e 

J A 



-Sir, 
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does not suffer from local gauge invariance, as does a naive integration. 

If the gauge condition (Eq. (1.12)) is ambiguous, it is termed non-ideal 
and an integration beyond perturbation theory may become ill-defined. The 
different solutions to a non-ideal gauge condition belong to the same orbit and 
are called Gribov copies in honor of its discoverer [Gri78]. In the following, 
we assume the gauge condition to be ideal. Note that even popular non-ideal 
gauge conditions, like the Coulomb or Landau gauge, are sufficient within 
the framework of perturbation theory. This is because in perturbation theory 
only small fluctuations of A^(a;) around zero are necessary and with respect 
to infinitesimal gauge transformations 

5^(x;t) = l + iTu;"(x)T« + 0(T2) (r < 1). 

even non-ideal gauge conditions are unique. The problems of Gribov copies 
and nonperturbative quantization will be discussed in Sec. 2.1.^ 



1.1.3 The Faddeev- Popov method 

From ordinary calculus of discrete n-dimensional vectors it is known that 



5"(/) = I 



5"(/M) 



det 



^f^ 



(1.13) 



where the determinant in the last expression is the Jacobian determinant that 
arise due to the substitution rule for integrals with multiple variables. If / 
is invertible near cu then its Jacobian determinant at cu is non-zero (inverse 
function theorem). 

If we assume in the following that Eq. (1.12) represents an ideal gauge 
condition, the identity Eq. (1.13) may be generalized to an identity for func- 
tional integrals 

1 = j [Dcu] 6 [J^i^'A]) Afp[A] (1.14) 

which was first proposed by FADDEEV and PoPOV [FP67]. In this relation, 
the Jacobian determinant'^ 



Afp[A] :=detM[A] 



(1-15) 

^Gauge-fixing is also necessarj' in canonical quantization, but not for stochastic quan- 
tization. Therefore the latter has the advantage to do not suffer from Gribov copies. 
However, it is more complicated than the other two methods. For the standard lattice 
approach to QCD gauge-fixing is also not necessary. See also Chapt. 2 and 3. 

■^In general the absolute value of the FP determinant has to be considered. However, 
the assumption of an ideal gauge condition guarantees the determinant to be nonzero. So 
it cannot change sign which cancels anyway due to normalization. See also the discussion 
in Sec. 2.1.2. 
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is known as the Faddeev-Popov (FP) determinant [FP67] of a matrix 

(1.16) 



u)=0 



that represents the change of JF under local gauge transformation at a; = 0. 
Inserting the identity in Eq. (1.14) now in the naive integration over gluon 
fields we end up with 

J [DA] App[A] S{J^[A]) e-^'-[^l (1.17) 

This represents an integration over the fundamental modular region, but if 
and only if the gauge condition Eq. (1.12) is unique. 

1.1.4 An effective Lagrangian density in covariant 
gauge 

A popular gauge condition for practical calculations is given by the family of 
covariant gauges specified by the condition 

T'^IA] d.A^ix) - B\x) = 0. (1.18) 

Here B"'{x) is an arbitrary fimction^. Since the work of Gribov [Gri78] it 
is well-known that local gauge conditions of this type are ambiguous with 
respect to finite gauge transformations (see also [Sin78; Wil03]), and so, in 
our notation, belongs to the class of non-ideal gauge conditions. With respect 
to infinitesimal gauge transformations, however, they are unique and may be 
treated as ideal ones. 

The family of covariant gauges turns out to be useful in perturbative 
expansions in many applications. It also allows us to represent the delta 
function in the functional integral Eq. (1.17) as a (functional) integral over 
the fields B. Since these fields are arbitrary we can use a Gaussian weight of 
width ^0 to integrate over, i.e. 

J DB" exp J d'^x (5")2(a;)| 5 {d^A^^x) - B^^x)) =: e-^of[^'«ol . 

Here S'cFl^iCo] is defined as the (Euclidean) space-time integral of 



^In Minkowski space -B" would transform as a Lorentz scalar. 
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Cqp is known as the gauge- fixing term which is added to the invariant La- 
grangian density Ci^y. It serves as a substitute for the delta- function that 
specifies the hypersurface in the functional integral Eq. (1.17). The param- 
eter ^0 is the gauge parameter that specifies the particular gauge condition 
in the family of covariant gauges. The special case of = is known as the 
Landau or Lorentz gauge and = 1 as the Feynman gauge. 



Never call a ghost stupid — A few good ghosts can help 

Although the FP determinant (Eq. (1.15)) is not a local function of the 
gauge fields, the functional integration in Eq. (1.17) can be extended such 
that the FP determinant is expressed by an additional (local) term, Cfp, 
added to >Cinv, too. This is done by the familiar device of integrating^ over 
ghost and anti-ghost fields c and c which are independent Grassmann valued 
fields. With the definition of the FP matrix M (Eq. (1.16)) we obtain in 
covariant gauge 



a;=0 



a;=0 



= -9tDl\W\x-y). (1.20) 

Here D"^^^ denotes the covariant derivative in the adjoint representation. Us- 
ing this, the term £fp, known as the ghost term, takes the form 

>Cfp = -{d,^){d'^5'^' + g^r'^ADc" . (1.21) 



Generating functional for QCD in covariant gauge 

In summary, we arrive at an effective Lagrangian density 



(1.22) 



where the individual terms Anv, J^gf and >C|p are defined in Eq. (1.3), (1.19) 
and (1.21), respectively. £eff can be used to define the (Euclidean) generating 



^For any finite N it holds that the determinant of a matrix M can be expressed as a 
functional integral over anti-commuting Grassmann numbers, i.e. 



detM 



N 

JJ / dctdci 



-CjMjiCi 
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functional (see e.g. [AvSOl]) 

Z[f, 77, 77,^, ^] = j [DA] [D^] [D^] [Dc] [Dc] 

• exp |- y d'x {CW - AIjI - j^V' - - ac - ca) I . (1.23) 

Here f), rj, a and a refer to the Grassmannian sources, respectively, for the 
quark, anti-quark, ghost and anti-ghost fields as introduced above. In covari- 
ant perturbation theory this generating functional is used for the calculation 
of Euclidean Green's functions as power series expansions of the interaction 
terms in £eff- Actually, for this the renormalized effective Lagrangian den- 
sity £gg given in Eq. (1.27) must be used instead. Otherwise perturbative 
expansions beyond tree level would be rendered meaningless by divergent 
mathematical expressions. We have indicate this already in Eq. (1.23) by 
giving the suffix r to jCgs. The explicit form of C^^ and the renormalization 
program is discussed in Sec. 1.2. 

Note also that the existence of the generating functional beyond pertur- 
bation theory rather has the status of being postulated than confirmed. So 
far only the continuum limit of a lattice formulation of quantum field theory 
provides a safe definition of the measure in the Euclidean generating func- 
tional, and thus the Euclidean Green's functions as its moments [AvSOl]. 

Under the assumption of its existence, vacuum expectation values of ob- 
servables are obtained from the generating functional as functional deriva- 
tives. For a general observable denoted as O this yields 

(O) oc J [VA][Vi;][V^P][Vc][Vc] O e-^eft[A,vi,v,c,c] ^^ ^4) 

where S^s denotes the effective action, i.e. the space-time integral of the 
efi'ective Lagrangian density. Note that this density also depends on the 
gauge parameter ^q. 



Gauge independence and gauge invariance 

For gauge-invariant observables it can be shown (see e.g. [Wei96]) that func- 
tional integrals of type as given in Eq. (1.24) are independent (within broad 
limits) of the gauge-fixing functional i.e. of the gauge condition. The dif- 
ferent types only result in irrelevant constant factors which are normalized 
away in the ratio of functional integrals. On the contrary, vaciuim expecta- 
tion values for gauge-variant observables dependent on the gauge condition. 
We found it worth to quote in this context a note from Collins's book [Col84]: 
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It is important to distinguish the concepts of gauge invariance and 
gauge independence. Gauge invariance is a property of a classical 
quantity and is invariance under gauge transformations. Gauge inde- 
pendence is a property of a quantum quantity when quantization is 
done by fixing the gauge. It is independence of the method of gauge 
fixing. Gauge invariance implies gauge independence, but only if the 
gauge fixing is done properly [Col84, p.31f]. 

1.1.5 The BRST formalism 

In the last section we ended up with an effective Lagrangian density Eq. (1.22) 
that is no longer local gauge invariant. Local gauge invariance spoils a naive 
integration over gauge fields and thus a gauge has to be fixed before the 
functional-integral formalism is applicable for quantization. However, it is 
a fundamental physical requirement that gauge-fixing is done in such a way 
that matrix elements between physical states are independent of the actual 
choice of gauge condition. The class of effective Lagrangians generated 
by the FP formalism above, can be shown to fulfill this requirement. They 
all yield the same unitary ^-matrix^. 

If phrased in a modern language of quantum field theory, namely the 
BRST formalism, the effective Lagrangian must be BRST invariant in order 
to have a renormalizable theory yielding a unitary S'-matrix. This formalism 
takes BRST invariance as a first principle and can be even used as a substi- 
tute for the FP method, in particular there where the FP method fails. ^ 

The BRST formalism goes back to the discovery of Becchi, Rouet 
and Stora [BRS75; BRS76] who first noted (independent also Tyutin 
[Tyu75; IT76]), that even if Cea is no longer locally gauge-invariant, it is 
invariant under a special type of global symmetry transformation. This sym- 
metry is a supersymmetry that involves ghost fields c"'{x) in an essential way. 
Remember, in the FP approach ghost fields are merely a technical device to 
express the FP determinant in terms of a path integral. In the BRST formal- 
ism, however, they serve as parameters cj'*(,t) = S\c"'[x) of infinitesimal local 
gauge transformations of gauge and fermion fields. Here SX is an (infinitesi- 
mal) ^-independent Grassmann number. In fact, a BRST transformation is 
isomorphic to an infinitesimal gauge transformation (Eq. (1.9)). They are 
written in the form 

5b$ = 5As<I> where $={A^,t/'}. 
''See e.g. [Wei96] for a more details 

^For example, the BRST symmetry is a basis for developing the canonical operator 
formalism. See the book by Nakanishi and Ojima [NO90] for a comprehensive account 
on that. 
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Here s denotes the BRST operator that acts upon the gauge and fermion 
fields according to 

sAl{x) = Dfc\x), (1.25a) 
s^(a;) = -igoTV(x)V^(a;) . (1.25b) 

Upon the ghost field c"(x) the BRST operator s is defined to act as 

sc"(a;) = -^rV(x)c'=(x) . (1.25c) 

This and the Jacobi identity for the structure constants /"^'^ suffice to show 
that the BRST operator s is nilpotent [KU82], i.e. 

s(s$) = 0. 

In addition to the ghost fields, the BRST formalism introduces antighost 
c" and Nakanishi-Lautrup auxiliary fields [Nak66] that transform as 

s^ix) = iB^ix) , (1.25d) 
sB^ix) = 0. (1.25e) 

The introduction of the auxiliary fields linearizes the BRST transforma- 
tions and renders the operator s to be nilpotent also off-shell. 



The BRST charge 

Since the BRST symmetry is a global symmetry^ of £eff there exists a cor- 
responding Noether current Jb that is conserved, i.e. d^JB^ = 0. Its explicit 
form (see e.g. [NO90]) is not of interest for the following discussions. But 
the existence of a corresponding unbroken charge is important (see be- 
low and the discussion concerning the Kugo-Ojima confinement criterion in 
Sec. 2.3.1). 

In general, the charge Q corresponding to a current J^, is defined as the 
spatial integral of Jq. It is a generator of the global symmetry, even if the 
integral is not convergent.^'' If this is the case, however, then the charge is 
ill-defined and is a generator of a spontaneously broken global symmetry. 

^Note that the invariance of Ci„y is a trivial consequence of its gauge invariance. Af- 
ter renormalization (see Sec. 1.2) the BRST symmetry is still a global symmetry of C^g 
(Eq. (1.27)) supposed the BRST transformations are substituted by the renormahzed ones 
given in Eq. (1.30) and (1.31). 

^°If the integral is not convergent, the charge Q is an ill-defined operator, and hence 
neither eigenstates nor expectation values of Q can be considered. However, for any local 
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It has been argued by KuGO and Ojima [K079] that the BRST charge 
Qb is an unbroken charge and so we can consider its eigenstates. In partic- 
ular, states belonging to the physical state space Vphys are assumed to be 
BRST singlet states of Qb, i.e. they are annihilated by 



This assumption plays an important role in the Kugo-Ojima confinement 
scenario to be introduced in Sec. 2.3.1. It is also related to the requirement 
of gauge-independence for physical matrix elements as we discuss now. 

Gauge-fixing and the BRST formalism 

Within the BRST formalism gauge-fixing is neatly performed by considering 
the BRST invariance as a first principle, i.e. a Lagrangian density has to 
be BRST invariant to have a renormalizable theory that yields a unitary 
S'-matrix. Since the BRST operator s is nilpotent, BRST-exact (or BRST- 
coboundary) terms — those are of the form s(*) — can freely be added to 
the gauge invariant Lagrangian density >Cinv, i.e. 



This will not change physics in any order of perturbation theory, but it can 



where was defined for covariant gauges in Eq. (1.18). 

Any change AT in the functional T, for example in JF", must not change 
any matrix element (\I'i|\I'2) of physical states [Wei96], i.e. 




= (*i|sAr|*2) . 



(1.26) 



With Qb being a generator of the BRST symmetry 



is(*) 



quantity $(?/) the (anti-) commutator 




can be considered, because the integrand vanishes for sufRcientljf large x. It follows that Q 
is an generator of infinitesimal symmetry transformations. See [NO90, pl3f.] for details. 
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we obtain that Eq. (1.26) can only hold for arbitrary changes in T if physical 
states are in the kernel of Qb, i-e. 

= (*i|[gB, Ar]|*2) ^ O^iIQb = Qb\iP2) = 

Note that the BRST symmetry of the full quantum Lagrangian is a basis 
for developing the canonical operator formalism. It also is very useful for 
deriving the Slavnov- Taylor-identities (STI). These are used for the proof of 
renormalizability of QCD. 

1.2 Regularization and renormalization 

The quantum field theory of the strong interaction introduced so far is still 
incomplete, because pertiu'bative expansions of Green's functions beyond tree 
level would be rendered meaningless by divergent mathematical expressions. 
In particular, loop integrals produce ultraviolet divergences when the cutoff 
for the internal momentum integral is send to infinity. This is known since 
the early days of QED (see e.g. [Opp30]). 

Fortunately, QCD is a renormalizable theory to any finite order in per- 
turbation theory.^^ That is, all divergences may be absorbed into a change 
of the normalization of the Green's functions and a suitable redefinition of 
all parameters appearing in the Lagrangian. The renormalized theory then 
yields only finite expressions at any order of perturbation theory. 

1.2.1 Regularization 

For this to work, QGD needs to be regularized prior to renormalization, 
using, for example, the Pauli-Villar, the dimensional or the lattice regular- 
ization. Actually, the latter is the only known nonperturbative regularization 
of QCD. In this regularization, the lattice spacing a serves as a regularization 
parameter that renders all momentum loop integrations finite via a gauge in- 
variant ultraviolet cutoff A = ar^. Consequently, arbitrary n-point functions 
calculated within the lattice approach due not suffer from ultraviolet diver- 
gences as long as a > 0. We shall briefly introduce the lattice regularization 
of QCD in Sec. 3.1. A list of references for detailed information can also be 
found there. 

^^The first proof for non-abelian gauge theories to be renormalizable was given by 
'tHoOFT and Veltman [tH71b; tH71a] using Slavnov- Taylor identities (STI) [Sla72; 
TayTl]. Modern proofs take advantage of the BRST formalism. 
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1.2.2 Renormalization 

After a suitable regularization has been carried out, for example by setting a 
cutoff A, the renormalization program introduces so called Z-factors which 
absorb the finite (because of the cutoff), but potentially divergent part in each 
of the fundamental two and three-point functions, i.e. those with tree-level 
counterpart in the Lagrangian density [AvSOl]. These are the inverse giuon, 
ghost and quark propagators as well as the three-gluon, four-gluon, ghost- 
gluon and quark-gluon vertices. The corresponding Z-factors are denoted by 
Zs, Z3, Z2, Zi, Z4, Z\ and Z\p respectively, and are formally introduced by 
writing the renormalized Lagrangian density as (see e.g. [AvSOl]) 

1 , ,2. 



+ Z,^d\' + Z.gJ'^^^^d, (A^c'') - Z.gJ'^'^ {d,Al) A'^At 
+ Z, ^ij'-'^r"' A^AlA^At + Z,^{- ^,d, + Z^m,)^ 

-z.Figri^i.T'^i^A;. (1.27) 

An additional factor, Z^, is necessary to adjust the mass of the quark prop- 
agator to the pole mass and g,. refers to the renormalized coupling constant. 
The latter is related to the bare parameter by considering, for example, 
the ghost-gluon vertex. Using this, 

Sr ^^go ■ (1-28) 

Of course, any other vertex function could be used instead to define g^. If all 
renormalization constants were independent then each vertex would define 
its own renormalized coupling constant. For example, using the three-gluon 
or the quark-gluon vertex this is 

^3 ^3 ^2 



9AAA '— —^go or g^^A '— —7^ go ■ 

But in order to guarantee that they all define the same coupling constant, 
i.e. g^ = qaaa = g^{i)Ai the Z-factors are constrained by the Slavnov- Taylor 
identities (STI) [Sla72; TayTl] giving: 

^ = £i = ^ = ^=: ZX" . (1.29) 

Z3 Zs Z2 Z, ' ^ ^ ^ 

Then the renormalized coupling constant is universal and can be related to 
the bare coupling constant gg by Z^g^ = gp where Zg is defined in Eq. (1.29). 
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Obviously, the renormalization constants are not independent of each other 
and the STIs allow us to express the constants for the vertices by the field 
renormalization constants Z3, Z3, Z2 and an independent one Zg. 

Comparing Eq. (1.27) and (1.22), we see that the renormalized Lagrangian 
(Eq. (1.27)) is related to its bare expression (Eq. (1.22)) by rescaling the fields 

Al ^ Z'J^AI ^ ^ zy\, ^ Zl'\'^ 

and by redefining the parameters appearing in the Lagrangian: 

So = Zgg^, mo = ZmTTlr, Co = Z^Cr ■ 

In an analogous manner, this translates to the BRST transformations 
introduced in Sec. 1.1.5. To be specific, the renormalized Lagrangian density 
is BRST invariant and all considerations made previously remain valid 
if the (infinitesimal) parameter of the BRST transformation is replaced by 

SXr := Z^^^^Z^^^^SX (1.30) 
and the renormalized fields transform as (see e.g. [AvSOl]) 

SrA^ix) = Z^Dfc\x), (1.31a) 
Sri>{x) = -Ziig^T'^c'^^ix) (1.31b) 



s,c«(x) = -Zi|:r'"=cV(x) , (1.31c) 
Sr^{x) - U^^A^ix). (1.31d) 



1.2.3 The MOM scheme 

After this rather formal rescaling and renaming of fields and parameters, any 
unrenormalized, but regularized n-point or Green's function is related 
to their renormalized one (in momentum space) through 

Gr{pi: ■■■,Pn] gr, '^r, Cr) = Zq ■ G,eg{pi: ...,Pn]K go, Co, mo) (1.32) 

where Zq refers to the corresponding product of Z-factors that appear in 
the Green's function. To give two simple examples: for the gluon two-point 
function in momentum space D'^^{p) :~ (A^(p)A^(— p)) this is Zq = Z^, 

~ ~ 1/2 

whereas for the ghost-gluon vertex Zq is given by Zq = Zi = Z^Z^^Zr^ . 

The Z-factors have to be determined such that the renormalized expres- 
sion on the left hand side of Eq. (1.32) is finite. The way this is done is 
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defined by the renormalization scheme. There are different renormalization 
schemes. In each scheme the divergent part is absorbed into the Z-factors, 
but they differ in how much of the finite part is absorbed, too [Mut98]. Com- 
mon renormalization schemes are the subtraction schemes MS, MS or MOM. 
The latter type is considered in this thesis, even though there are infinite 
many different MOM schemes. 

A MOM scheme defines the Z-factors such that the fundamental two-point 
and three-point functions equal their corresponding tree-level expressions at 
some momentum fj,'^, the renormalization point. 

The two-point functions that are relevant in this thesis are the gluon 
propagator (A^(a;)A^(|/)) and the ghost propagator {c"'{x)c^{y)) in Landau 
gauge. Actually, we are interested in the Fourier transform of these two 
expressions. In momentum space the gluon propagator in Landau gauge has 
the following tensor structure: 



/^) = ( - ^) (1-33) 



where Z denotes the form factor or the dressing function of the gluon prop- 
agator. It expresses the deviation of D^^{p) from its tree-level form [Z = 1). 
For the ghost propagator the corresponding tensor structure is given by 

Here J denotes the dressing function of the ghost propagator. 

In a MOM scheme the renormalization constants, for instance of the gluon 
and ghost fields Z^ and Z3, are defined by requiring the renormalized expres- 
sions to equal their tree-level form at some (large) momentum //^. That is, 
Z3 is defined as 

d;1{p; A, go, ^0,^0) ^ ,=--Z3S^'(s^''-^)^ (1.35) 
p-t^^ \ 11/11 

where denotes the unrenormalized gluon propagator. Z3 is given by 



HI/ 

-tab 



G^\p;A,g,,mo,^o) ^ =: Z3 (1-36) 
■=f^^ II 



Therefore, a renormalization constant can be determined by calculating 

the corresponding unrenormalized (regularized) Green's function. Its value 
depends on the renormalization point /i^ and also on the bare parameters 
of the regularized theory. For example: In our lattice simulations we have 
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calculated the bare (quenched) gluon propagator using the bare parameters: 
gQ(A^) = 1, 1/mo = and = 0. Requiring Eq. (1.35) to hold at some 
momentum /x^, we have fixed Z3. The renormalized gluon propagator (at /x^) 
is then obtained via multiplicative renormalization according to Eq. (1.32). 

1.3 The renormalization group 

Obviously, a renormalized Green's function depends on the subtraction point 
H whose choice is not unique. Also the renormalized parameters g^, rrir and 
depend (via the corresponding Z-factors) on /x. Keeping gg, mo, Co and A 
fixed, we could, of course, had chosen another point, say //'. This would yield 
a new renormalized Green's function with the new values g^l/^O) "^ril-^') and 
Cr{l^')- Even though both renormalized Green's functions are different, they 
are related by a finite multiplicative renormalization, i.e. 

Gr{pi; gr(l^'),mr{lJ,'),^r(l^'),l^') = ^(//', //) Gr{pi; g^(//) , m^(//) , Cr(//) , /i) 

(1.37) 

where -2 is a finite number depending on /i and f^'. The finite renormalization 
of Green's functions forms an Abelian group called the renormalization group 
(RG). Physically measurable quantities are invariant under renormalization 
group transformations, i.e. they are independent of the subtraction point jj,. 
Green's functions, are generally not renormalization-group invariant. An 
analytic expression of this property is given by the renormalization group 
equation [Mut98]. 

1.3.1 The renormahzation group equation 

The renormalization group equation is best derived by noting that the un- 
renormalized Green's function does not depend on the renormalization point 
II if all bare parameters (go, mo, Co and A) are fixed, i.e. 

= At^G-bi; So: Co, mo, A) . (1.38) 

On the contrary, the renormalized Green's function depends on n not only 
explicitly, but also implicitly due to the renormalized parameters. By us- 
ing the chain rule for differentiation, Eq. (1.38) yields for the renormalized 
Green's function 
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Here a sum over the different fermion fiavors is implied and the (dimension- 
less) RG functions are defined as (see e.g. [Mut98]) 



P gr: >Cr 

/i 

Ira \ Eri ' 

nir 

7 I gr^ ^Cr 

fj, 



.dgr 



g(),"Xo,Co,A fixed 



rrir dfi 
dlnZc 



djjL 



go ■'"o,^o,A fixed 



gOi^^o.^cA fixed 



djji 



(1.40a) 
(1.40b) 
(1.40c) 
(1.40d) 



gOi'^OiCcAfixed 



The RG equation expresses how the renormalized Green's function, in partic- 
ular their parameters, change under a variation of the renormalization point 
11. In the following we shall assume that ^ rrir always holds. Thus ap- 
proximately, the RG fimctions do not depend on m^.^^ The /3- function has 
also been proven to be gauge independent [MP78], i.e. 

To get rid of the gauge dependence of the other RG functions, we shall 
restrict ourselves in the following to the Landau gauge, because the Landau 
gauge is a fixed point under the renormalization group. To see this note that 
in general covariant gauge we obtain for the RG function (3^ (Eq. (1.40d)) 
depending on the gauge parameter {Z^^r — Co) 



go ,^0, A fixed 



^-zl^- 



dlnZ-, 



go ,^0, A fixed 



Therefore, given the initial condition = the function (3^ vanishes com- 
pletely in Landau gauge. 

We are left with the three RG equations P{gj.), jm{gr) aiid ^{gr) (see 
Eq. (1.40a), (1.40b) and (1.40c)) which in our approximation only depend on 
the renormalized coupling constant g^ilj)- They express how the renormal- 
ized parameters gri^lj), mr{[x) and the renormalized Green's function change 
under a variation of ji. In fact, given the initial values gr{li) and mr{ii) at a 
renormalization point /i, the values g^l/^') mr{ii') at /x', are determined 



^^Siich an approximation corresponds to using a mass-independent renormalization 
scheme, like for example the MS scheme. 
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by the solution of the differential equations (1.40a) and (1.40b), respectively. 
That is 



Similarly, the change of the renormalized Green's function under the RG 
transformation — > /x' is obtained. To see this, consider Eq. (1.37). The RG 
equation (1.39) tells us that [Col84] 



The function 7 is known as the anomalous dimension for reasons that become 
clear in Sec. 1.3.4. A full solution to this RG equation is given by Eq. (1.37) 
where [Col84] 



The values of gr{l^') and rrirdi') are given in Eq. (1.41) and (1-42), respec- 
tively. 

Of course, the explicit form of the RG functions are generally unknown, 
but approximations can be made by taking a finite number of terms in the 
perturbation series for the RG functions 7^ and 7. Since this is an expan- 
sion in gj. it is only valid at sufficiently small g^. We shall see in Sec. 1.3.3 
that for QCD this is realized in the asymptotic region of large Euclidean 
momenta. 

Therefore, the most important application of the RG equation in QCD 
is to study the asymptotic behavior of Green's functions at large Euclidean 
momentum [Col84]. 

1.3.2 Perturbative expansion of the /3-function 

For the /3-function (1.40a) a power expansion in the coupling constant can 
be calculated by choosing one of the four vertex functions which define g^. 
Taking, for example, the three-gluon vertex the renormalized coupling con- 
stant gj. is related to the bare coupling as 




(1.41) 



(1.42) 





(1.43) 
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The renormalization constants, and Zi, are defined at a subtraction point 
— //^ in terms of the bare (transverse) gluon propagator and the bare three- 
point vertex, respectively. Extracting both Z-factors to two-loop order, the 
solutions can be plugged into the RG equation Eq. (1.40a) for g,.{p). After 
some algebra this gives an expansion for the /^-function (see e.g. [Wei96]) 

/5[&(/^)] = -/3o^ - + O [glii,)) (1.44) 

that holds at small g^ifJ^). The first two coefficients are given by 

Po = 11 -^AT^, (l-45a) 
19 

Pi = 51-— Nf. (1.45b) 

Here Nf is the number of quark flavors with masses below the energies of 
interest. 

Since the determination of the Z-factors generally depends on the renor- 
malization scheme used, the explicit form of /5(g;.) depends on the gauge 
and on how the running coupling is precisely defined. However, it can be 
shown that the first two coefficients, Po and Pi, are renormalization-scheme 
independent, whereas those of higher loop-expansions are scheme-depended. 



1.3.3 The running coupling constant 

Using the expansion of the P function we can solve the RG equation Eq. (1.40a) 
for the coupling constant grifJ-) to the given order. The general solution to 
Eq. (1.40a) takes the form given in Eq. (1.41). As mentioned above, it 
describes the variation of g^ under the change /i ^ /i' keeping the bare pa- 
rameters gQ, mo and A fixed. This is usually termed as the running of the 
coupling constant g,. changing the energy or the momentum scale jj,. 

It is common practice to parameterize the running coupling constant^^ 
gril^) by introducing a RG-invariant mass parameter A being the integration 
constant of a solution to the differential equation in Eq. (1.40a). It is defined 

^^Note that in each energy range between any two successive quark masses we have a 
difTerent value of Nf, and also a difTerent A (see Sec. 1.3.3), chosen to make g{fj.) continuous 
at each quark mass.[Wei96, p. 157] 

^"'In the literature sometimes the notation ^[g'^(/x), ln(;[i//x')] is used for the running 
coupHng constant [GW73a; GW73b; Pol73]. It is the same as gril-i-'), because g is defined 
to be the value of the coupling constant renormalized at jj.' (what we call g^(/i')) if it is 
known to have the value g^il^) at [CG79]. 
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by 



A := 



-bl/(2b§) 



• exp 



rilA 



dh 



hi 



p{h) boh^ hlh 



1 1 

+ 



. (1.46) 



Here ho and hi are nothing but the first two coefficients of the /9-function 
given in Eq. (1.44), i.e. 

/3o , . _ /3i 



and hi = 



(1.47) 



167r2 12877^ 
Specifying g^lA*) value of /i is exactly equivalent to fixing A [Col84]. 

Thus if the /3-function were known we could calculate A from the knowledge of 
grifJ') at one value /i and vice versa. Note that renormalization has introduced 
a new parameter A of dimension mass into the theory that is not present in 
the bare Lagrangian density. This is called dimensional transmutation. The 
parameter should be determined by comparing experimental data with QCD 
predictions. 

The definition of A in Eq. (1.46) is scheme dependent, since the coupling 
constant gr(yu) may have even different meanings in different renormalization 
schemes. For example, on the lattice a scale is given by the lattice spacing a 
and the coupling is given by the bare coupling constant go(o) depending on a. 
The scheme dependent parameter A on the lattice is called Alat- In the MS 
scheme the A-parameter is usually called A^, whereas in MOM scheme it is 
called Amom- The knowledge of the ratio allows to relate results obtained in 
different renormalization schemes. For example the ratio of Amqm and A-j^ig 
is given by (see e.g. [CG79; MM94]) 

Amom 



2.895655. 



*MS 



Up to two-loop order the coefficients of the /9-function are scheme in- 
dependent. Using these coefficients an expression for the running coupling 
constant can be given that is valid in any renormalization scheme as long as 
/X ^ A is fulfilled in that scheme. Defining 



47r 



(1.48) 



the solution for the running coupling constant as{n) in two-loop order is (see 
e.g. [Col84; E+04]) 



47r 



/3oln(/.VA') 



2A ln[ln(/iVA^)] 
(31 ln(/xVA^) 



O 



ln^[ln(/i/A)] \ 
ln'(/^/A) )' 



(1.49) 
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This two-loop result for q;s(/x) will be used in Sec. 4.4.1. 

We note in passing that a similar expression can be derived for the running 
mass, but since this will not be a subject of this thesis we refer to standard 
textbooks for it. 

1.3.4 The anomalous dimension 

It is interesting to have a look once more at Eq. (1.37) and (1.43). Assuming 
a renormalized Green's function Gr has been computed at the momentum 
\p using a large value for the scale factor A. Setting /i' — X/i, from Eq. (1.37) 
one knows that 

Under the assumption rrir^Xii) does not get large and using dimensional anal- 
ysis'^ one obtains [Col84] 

= • z{f,,Xfj,) ■ Gr{pi;gr{XiJ,),0,iJ,) (1.50) 

where D is the dimension of G. This makes it evident that the relevant cou- 
pling constant is the effective coupling at the scale of the momenta involved. 
A change in all external momenta using a common scale factor A is equivalent 
to changing the coupling constant. We see further from Eq. (1.50) that the 
overall scale factor is not just A^, as one might expect from naive dimen- 
sional analysis (i.e. 3{g,.) = 0), but that it includes an extra factor z defined 
in Eq. (1.43). This is the reason why 7 is called the anomalous dimension. 
Note that this dimension arises from the fact that a scale changes the renor- 
malization point, and G is not necessarily invariant under this operation. 

To lowest order in perturbation theory the anomalous dimension ^(gr) is 
given by the expansion j{g,.) = CQgf. + O (gf) [Col84] where Cq is the zeroth- 
order coefficient of the anomalous dimension. Using cq and the corresponding 

^^If one scales the momenta Pi — > Xpi then using dimensional arguments, a Green's 
function behaves as: 

Gr{Xpi,g^, n) = M^/(A^pi • Pj/n'^) 

where D is the dimension of G and / is dimensionless. This is because G is Lorentz 
invariant, and hence can only be a function of the various dot products Pi ■ pj [Kak93]. 
Hence 

Gr{Xpi,g^,Xlx) = A^/X^/(ft -Pj/li^) = X^Gr{pi,gr,H)- 
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coefRcient bo (Eq. (1.47)) of the /3-function, we obtain from the definition of 
z{X) — z{Xii, /i) (Eq. (1.43)) to lowest order in perturbation theory [Col84] 

oc [In A] ~ A — > oo 

where 5 := Co/(26o)- For the gluon and ghost propagators in the quenched 
case {Nf — 0) these exponents are Sd — 13/22 and Sq — Q/4:A, respectively. 
Therefore, the corresponding dressing functions, Z and J, behave in the far 
ultraviolet momentum region like 

Z{p') ~ (in^) and J{p') ~ (in ^) . 



" gr(V) 



-| co/bo 



[1 + 



(1.51) 
(1.52) 
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Chapter 2 



Infrared QCD and criteria 

for confinement 



e discuss briefly some problems and recent developments concerning non- 



V V perturbative quantization with a particular focus on the infrared region of 
QCD. We give a summary of results obtained by using the Dyson-Schwinger (DS) 
approach. Thereby we restrict ourselves to those results which are of relevance 
for the discussion of our lattice results presented in subsequent sections. In the 
second part of this chapter we introduce different criteria for confinement that can 
be formulated for QCD in covariant gauges. Note that in Chapt. 5 we shall try to 
confirm these criteria for lattice QCD in Landau gauge. 



2.1 Nonperturbative approaches to QCD 

It should be clear from the discussion of the renormalization group that it is 
possible to define an effective (running) coupling constant gj. which is a func- 
tion of the momentum. This functional dependence is intimately related to 
the momentum dependence of the vertex functions of QCD. For sufficiently 
large Euclidean momenta q, the coupling constant becomes small and thus 
perturbation theory seems an appropriate calculation tool in this asymptotic 
regime. On the other hand, gj. grows large with decreasing Euclidean mo- 
mentum and diverges at a certain momentum, the so called Landau pole. Its 
position signals the breakdown of perturbation theory. Even though signifi- 
cant progress has been achieved in recent years in the perturbative calculation 
of higher order corrections to renormalization group functions, like, for ex- 
ample, for the P function or the anomalous dimension^, perturbation theory 
only applies to the region of large Euclidean momenta. Note that there is 
no imique definition of the running coupling constant, since its definition 
depends on the renormalization scheme employed [CG79]. 

^For details on the progress made and also for recent results see, for example, the work 
by Chetyrkin [CROO; Che97; Che05] and references therein. 
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2.1.1 Brief remarks on nonperturbative methods 

In any case, the infrared region corresponds to strong coupling rather than 
weak coupling, and hence is of interest for studying confinement of QCD. 
Since perturbation theory is of no avail in studying QCD at low momentum, 
this region can be explored only in genuinely nonperturbative approaches. 
One such approach is the vast field of Monte Carlo simulations of lattice 
QCD. It is a first principle approach that contains the full nonperturbative 
structure of QCD and has the striking feature of being manifestly gauge 
invariant. (For some details and a list of references see Chapt. 3). 

However, lattice simulations are limited by the enormous computational 
effort they require and thus its application to a wider range of unresolved 
problems connected with QCD can only be extended when computer technol- 
ogy continues to improve. Moreover, lattice calculations are always afflicted 
with uncertainties in extrapolating to the infinite volume and continuum 
limit which is necessary to connect with the physical world. Therefore, it is 
worthwhile to pursue other approaches that preserve features of QCD that 
the lattice formulation lacks. For example, lattice simulations cannot made 
definite statements about the far infrared region of QCD due to finite lattice 
volumes. Although, we did our best in this thesis to do so. 

In this context another complementary nonperturbative approach based 
on the infinite tower of Dyson-Schwinger equations (DSE) has gained much 
attention in recent years. In Sec. 2.2 we give a brief overview about the recent 
developments that are relevant for this thesis, but before we would like to 
stress a particularly important point concerning the problem of Gribov copies 
in the nonperturbative regime. 

2.1.2 The problem of Gribov copies 

The Dyson-Schwinger equations of QCD are derived from a generating func- 
tional. In the continuum, this requires gauge-fixing for the definition of an 
integration measure. For several reasons the gauge condition is mostly cho- 
sen within the family of covariant gauges which are known to be not unique 
beyond perturbation theory due to the existence of Gribov copies. Remem- 
ber that in the asymptotic regime, where perturbative QCD is relevant, the 
problem of Gribov copies could be safely ignored. This is because all Gribov 
copies {^A^} of a particular (gauge-fixed) gluon field carry an additional 
dependence — they are related by a particular local gauge transfor- 
mation (Eq. (1.7)) — and hence can be neglected within the framework of 
perturbation theory. Therefore, we can use the FP method or, even more 
elegantly, the BRST formalism for the class of covariant gauges to obtain a 
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quantized gauge theory which is manifestly Lorentz covariant and for gauge- 
invariant observables gauge-independent. Moreover, there are elegant BRST 
proofs of multiplicative renormalizability and unitarity to any order of per- 
turbation theory [Bau85]. 

Beyond perturbation theory we have to face the problem of Gribov copies. 
Actually, it represents the main impediment to nonperturbative gauge-fixing. 
As pointed out in [Wil03], there is no known Gribov-copy-free gauge fixing 
which is a local function of gluon fields A^. Therefore, one has to adopt 
either a nonlocal Gribov-copy free gauge or attempt to maintain local BRST 
invariance at the expense of admitting Gribov copies [Wil03]. Concerning the 
latter point, there is, however, the well-known Neuberger problem of pairs 
of Gribov copies with opposite sign giving expectation values 0/0 [Neu87; 
Tes98]. It is still unknown whether local BRST invariance for QCD can be 
maintained in the nonperturbative regime [Wil03]. 

Of course, it is desirable to overcome the Xeuberger problem and el- 
evate the BRST formalism to the nonperturbative level. A first success- 
ful step in that direction was done recently for the massive Curci-Ferrari 
model [KvSWOo]. Also in a recent work [GKW05], a generalization of the 
FP method was given that is valid beyond perturbation theory and, most 
notably, circumvents the Xeuberger problem. In fact, in [GKW05] a path in- 
tegral representation of the absolute value of the FP determinant in terms of 
auxiliary bosonic and Grassmann fields has been presented. Remember that 
usually the absolute value of the FP determinant is dropped, but this can- 
not be done beyond perturbation theory using a gauge condition that suffers 
from the Gribov ambiguity. The resulting gauge-fixing Lagrangian density is 
local and enjoys a larger extended BRST and anti-BRST symmetry, though 
it cannot be represented as a BRST exact object [GKW05]. 

2.1.3 Nonperturbative quantization in Landau gauge 

Moreover, progress has been made concerning the nonperturbative quanti- 
zation in Landau gauge. It fact, it has been argued by ZwANZlGER that an 
exact nonperturbative quantization of continuum gauge theory is provided 
by the FP method in the Landau gauge if the integration is restricted to the 
Gribov region fl [Zwa04]. That is, the vacuum expectation value of a gluonic 
observable O is given by 



^ JQ 

Here we have indicated the restriction to the Gribov region by giving a suffix 
il to the expectation value. 




(2.1) 
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To remind the reader, the Gribov region Q is defined as the region in the 
space of transverse gluon fields where the FP operator, M[A] = — 
is positive, i.e. all its (nontrivial) eigenvalues are positive. 



It can be shown that the Gribov region is convex, it is bounded in every 
direction and it contains the origin = (see e.g. [Zwa04]). The boundary 
dD, of the Gribov region is called the Gribov horizon. It consists of transverse 
gluon fields for which the lowest nontrivial eigenvalue of the FP operator 
vanishes. Thus the FP determinant, being the product of all eigenvalues 
detM = rin'^"' positive inside the Gribov horizon and vanishes on it. 

Due to the fact that the Gribov region is not free of Gribov copies [Zwa04, 
and references therein], Eq. (2.1) was generally abandoned to be used for an 
exact quantization in favor of an integration over the fundamental modular 
region A [Zwa94] 



(See Ref. [Zwa04] for a discussion). The fundamental modular region is 
the set of unique representatives of every gauge orbit, i.e. its interior is free 
of Gribov copies. The boundary of the fundamental modular region, dA, 
contains different points that are Gribov copies of each other, but which 
have to be identified topologically [vB92]. Like Q, the fundamental modular 
region is convex, bounded in every direction and included in the Gribov 
region [Zwa82; Zwa92]. The boundary of the Gribov region, i.e. the Gribov 
horizon dQ, touches dA at so called singular boundary points. 

It is diflficult to give an explicit description how to restrict the functional 
integration to the fundamental modular region. One usually introduces the 
gauge functional 



to characterize the different regions. Here denotes a gluon field lo- 
cally gauge-transformed by g (see Eq. (1.7a)). The gauge transformation is 
parameterized by uj as given in Eq. (1.8). The Gribov region corresponds 
to the set of all relative minima of FA[g] with respect to local gauge transfor- 
mations g. Those minima which are absolute characterize the fundamental 
modular region [Zwa04]. 



Q := {A^ : d^A^ = 0; -d^D^[A] > 0} . 




(2.2) 




FA[l]-i / d'^xTrujd^A^+ / d'^x Tr cud /.D^lu + O (uj^) 



(2.3) 
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Returning to Eq. (2.1) and (2.2), it is of advantage if the integration 
needs not to be restricted to the fundamental modular region. That this 
could be the case, has been argued recently by Zwanziger in [Zwa04]. In 
fact, in this reference it is stated that functional integrals are dominated by 
the common boundary of A and Thus, in the continuum, expectation 
values of correlation functions over the fundamental modular region A are 
equal to those over the Gribov region Jl, i.e. 

(0)a = {0)a , (2.4) 

even though the Gribov region is not free of Gribov copies. The Gribov 
copies inside fl are argued to not affect expectation values in the continuum. 

This is fortunate for many reasons. For example, in Monte Carlo sim- 
ulations of lattice QCD gauge-fixing is usually implemented by an iterative 
maximization (or minimization depending on conventions) of a lattice ana- 
logue to the gauge functional^ in Eq. (2.3). Therefore, lattice gauge-fixing 
automatically restricts to the Gribov region. Note, however, that it has also 
been pointed out in [Zwa04] that on a finite lattice the distinction between 
the fundamental modular region and the Gribov region cannot be ignored, 
but hopefully in the thermodynamic limit the relation Eq. (2.4) holds. We 
have found indications in our data [SIMPS05d] that support this assumption 
to be true. See Sec. 4.2.3 for more details. 

For our discussion in the next section it is worthwhile to mention that if 
the functional integral is cut off at the (first) Gribov horizon, both the (Eu- 
clidean) gluon and propagators have to be positive [Zwa03b]. This is fulfilled 
by the solutions obtained for truncated systems of Dyson-Schwinger equa- 
tions for the gluon and ghost propagators summarized in the next section. 
Therefore, they automatically restrict to the Gribov region, even though no 
direct restriction to dfl has been done [Zwa04]. 

2.2 The Dyson Schwinger equations of QCD 

Before we start to introduce the lattice approach to QCD (Chapt. 3) and 
report on the results we have obtained in studying some infrared properties 
of Landau gauge gluodynamics (Chapt. 4 - 6), we summarize results achieved 
in recent years within the DSE approach to QCD. Those results presented 
here will then be compared to our lattice data in subsequent chapters. 

The Dyson-Schwinger equations (DSEs) of QCD are infinite towers of 
coupled nonlinear integral equations relating different Green's functions of 

^See Eq. (3.14) for a definition of the lattice gauge functional as used in this study. 
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QCD to each other (see below for two examples). Generally speaking, the 
DSEs follow from the observation that an integral of a total derivate van- 
ishes. The DSEs can be directly derived using the generating functional 
given in Eq. (1.23) where the existence of a well-defined integration measure 
is assumed. 

The relevant DSEs for the subsequent discussion are those for the ghost 
and gluon propagators. Therefore, these are briefiy recalled, but for an ex- 
plicit derivation of all of the DSEs of QCD we refer to Ref. [AvSOl] (see also 
[RW94]). Following [AvSOl], the DSE of the (inverse) ghost propagator takes 
the form 

iG''r\k) = -6^%e + dr'Z, J ^zfc^G-(g)P/(g, k)D%{k - q) 

in momentum space. G denotes the full ghost propagator, D the full gluon 
propagator and F^-^* the full ghost-gluon vertex. For other notations used 
herein we refer to Chapt. 1. In Fig. 2.1 we also give a pictorial representation 
of the ghost DSE. This DSE contains the inverse of the tree-level propagator 
(dashed line), the tree-level ghost-gluon vertex (small filled circle) and the 
full ghost-gluon vertex (open circle). The latter is coupled to a fully dressed 
ghost and gluon propagator (dashed and wiggled lines carrying a filled circle, 
respectively). 

The gluon DSE is obtained ^S^a^ 

in a similar manner as the ghost ___^__r^ ^ _ 

DSE, but in comparison to that, 
it is much more complex. There- 

„ 1 • • . • 1 Figure 2.1: The ghost DSE [AvSOll. 

tore, we only give a pictorial 

representation of the gluon DSE in Fig. 2.2. Please refer again to Ref. [AvSOl] 
for an explicit expression. As in Fig. 2.1, wiggled lines refer to the gluon prop- 
agator in Fig. 2.2, whereas a line and a dashed line refer to the quark and 
ghost propagator, respectively. Lines carrying a full circle correspond to fully 
dressed propagators. Open circles denote full vertex functions. 

For the DSE of the quark propagator a similar diagram as for the ghost 
DSE can be given. See Ref. [AvSOl] for both an explicit and a pictorial 
representation. 

2.2.1 Infrared behavior of ghost and gluon propagators 
in Landau gauge 

If the full solution of QCD's Dyson-Schwinger equations were available it 
would provide us with a solution to QCD. However, this has not been ob- 
tained so far. The main impediment is that the infinite tower of coupled 
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Figure 2.2: Diagrammatic representation of the gluon DSE. Adapted from Ref. [AvSOl]. 



nonlinear integral DSEs has to be truncated in order to be manageable. 
Thereby, the challenge is to use suitable truncations that respect as much as 
possible the symmetries of the theory. Doing so, solutions of the truncated 
systems might provide deep insights into many phenomena in hadron physics, 
because they may serve as input into bound state calculations based on the 
Bethe-Salpeter equations for mesons or the Faddeev equations for baryons. 

Studying DSEs has become a subject of growing interest in recent years. 
See Refs. [RW94; RSOO; AvSOl; MR03] for a comprehensive overview. In par- 
ticular, in Landau gauge or general covariant gauges considerable progress 
has been made in studying the low-momentum region for the coupled sys- 
tem of quark, gluon and ghost propagators. In contrast to earlier attempts 
[Man79; ADJS81; AJS82; BP89] where contributions of ghost fields were ne- 
glected, the more recent attempts, initiated by VON Smekal et al. [vSAH97; 
vSHA98], have shown that the inclusion of ghost fields is important for the 
generation of a consistent infrared behavior of QCD [BloOl]. In fact, it has 
been shown in [vSAH97; vSHA98] that in the infrared momentum region a 
diverging ghost propagator is intimately related to a vanishing gluon propa- 
gator, both following a power law at low momentum. 

We have given explicit expressions for the gluon and ghost propagators in 
Eq. (1.33) and (1.34), respectively. The corresponding renormalized dressing 
functions have been denoted by Z and J. According to [vSAH97; vSHA98] 
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these dressing functions are predicted to follow the power laws 

Z{q\^^') (X {q'/^^'Y- (2.5a) 
J^q^li-") « (gV/^')-'^" (2.5b) 

in the limit — > with infrared exponents, and kq^ satisfying 

KD — '^K'G ■ (2.6) 

This has been confirmed later by Atkinson and Bloch [AB98b; AB98a] 
at the level in which the vertices in the DSEs are taken to be bare. Also 
investigations of the DSEs in flat Euclidean space-time performed by Fis- 
cher et al. [FAR02; Fis03; FA02] without angular approximations, as used 
in the former studies, support such an infrared behavior. The value of kq 
depends on the truncation used, but in Landau gauge it has been argued that 
Kg ~ 0.595 should be expected [LvS02; Zwa02]. Thus the ghost propagator 
is supposed to diverge stronger than and the gluon propagator to be 
vanishing in the infrared regime. 

These findings are, as we shall see in Sec. 2.3.1 and 2.3.2, in agree- 
ment with the Kugo-Ojima confinement criterion [K079] as well as with 
the Zwanziger-Gribov horizon condition [Zwa04; Zwa94; Gri78]. According 
to the latter condition, a diverging ghost and a vanishing gluon propagator 
result from restricting the gluon fields to the Gribov region Q [Zwa04]. 

Note, however, that quite recently Fischer et al. have investigated DSEs 
on a torus. They have found quantitative differences with the infinite volume 
results at small momenta [FGA06; FP06]. Moreover, Bloch has developed 
truncation schemes for the DSEs of the ghost, gluon and quark propagators 
[BloOl; Blo02] which preserve multiplicative renormalizability, in contrast 
to those schemes used in the references cited above. Furthermore, Bloch 
has demonstrated that a definite conclusion about the existence of infrared 
power-behaved gluon and ghost propagators cannot be reached yet. 

We shall show in Sec. 4.3 that our lattice data and also those by others do 
not confirm Eq. (2.6), at least for the lattice momenta available at present. 
The values of the infrared exponents have been found to be different. There- 
fore, we think it is still an open question whether the infrared behavior as 
favored by current DSE studies is realized for lattice QCD in Landau gauge. 
However, if it is realized then it has interesting consequences for the running 
coupling constant discussed next. 
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2.2.2 A nonperturbative running coupling constant 

From the renormalization of the ghost-gluon vertex, the renormalized cou- 
pling constant at the renormaUzation point ji is defined by 

z3(A^/.^)z|(A^/.^) 

as[^i ) = «(A ) ~ . (2.7) 

Here a{A?) := ^{A?)/At: denotes the bare coupling constant that depends on 
the ultraviolet cutoff A and the multiplicative constants, Z3, and Zi, refer 
to the renormalization constants of the gluon and ghost propagators and 
of the ghost-gluon vertex, respectively. Long time ago, it has been shown 
by Taylor [Tay71] that to any order in perturbation theory the ghost- 
gluon vertex in Landau gauge is finite (see Sec. 2.2.3 for some discussions). 
Therefore, one may chose _ 

Zl(A^//^) = l. (2.8) 

Based on this and by using the relation between the bare and renormalized 
ghost and gluon dressing functions 

J(g^/x2) = Z^\^^^^)JB{^\q^) (2.9a) 
Z{q^^?) = Z^\A',i,')ZBiA',q') (2.9b) 

where J{ii^,fj,'^) — Z^ji^^jj?) — 1, one can show that the product [vSAH97; 
vSHA98] 

as{q^) ■■= as{^^)Z{q\^^)j\q\i^^), (2.10) 

is renormalization group invariant and defines a running coupling constant 
within the context of DSEs [AvSOl; AFLE05]. In fact, with Eq. (2.7), (2.8) 
and (2.9) it holds that 

a,{^x')Z{q\^x'')J\q\^x'') = «(A2)Zs(A^g2)J|(A^g2). 

Obviously, the right hand side does not depend on the renormalization point 
//^ and thus the combination of (renormalized) ghost and gluon dressing 
functions on the left hand side is renormalization- group-invariant. Evaluating 
the left hand side once at an arbitrarily chosen ji^ and once at /x^ = q^ one 
obtains the product given in Eq. (2.10) [AFLE05]. 

The definition of a running coupling constant by Eq. (2.10) has been 
first derived in [vSAH97; vSHA98]. Later it has been shown by Bloch 
[BloOl; Blo02] that after a reformulation of the DSEs for the ghost, gluon 
and quark propagators this coupling constant enters directly the kernel of 
the DSEs. 
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Assuming multiplicative renormalizability to hold beyond perturbation 
theory, and assuming also the power laws (Eq. (2.5)) for the dressing func- 
tion being realized at infrared momenta according to the relation (2.6), then 
as{q'^) has finite infrared fixed point [vSAH97; vSHA98]. The precise value 
of as{Q) depends on kg, but under certain assumption it has been shown to 
be [LvS02] 

a,(0) ^ Smh/Nc for SU{N^). 

However, we would like to stress already here that our lattice data for 
as{q^ > 0) do not indicate the existence of a finite value at zero momenta 
(see Sec. 4.4.1). We have also found no indications for deviations from Zi 
being constant (Sec. 4.4.2). Our lattice data presented in this thesis are in 
agreement with other lattice studies and, most notably, also agree qualita- 
tively with recent investigation of DSEs on a torus [FGA06; FP06]. This, 
once more, puts the proposed form for infrared power laws into question. 

2.2.3 The finiteness of the ghost-gluon vertex 

Arguments were given first by TAYLOR [TayTl] that in Landau gauge the 
renormalization constant of the ghost-gluon vertex Zi = 1. He showed that 
in this gauge the ghost-gluon vertex is finite in the ultraviolet momentum 
region and stays bare for a vanishing incoming ghost momentum. To motivate 
this we recall arguments given in Refs. [MP78; AFLE05] using the DSE for 
the full ghost-gluon vertex. A pictorial representation of this DSE is shown 
in Fig. 2.3. 

Considering this figure, the argument goes at follows [MP78; AFLE05]: 
The bare ghost-gluon vertex in the interaction diagram (rightmost diagram) 
is proportional to the internal loop momentum Z^. Since in Landau gauge 
the gluon propagator D^^{1 — q) is transverse, it holds that l^D^^{l — q) = 
qixDfj,u{l — q) and thus the interaction diagram vanishes in the limit q^ ^ 
[AFLE05]. This implies (see Fig. 2.3) 

Tf^-(k, q) = ^A,.(^"''g' + B^'^e + C"^V') + D'^'^kxh^ + etc. 

as A;2, and so that r^Tf^''{k,q) vanishes [MP78]. Therefore, the 

full vertex equals the bare vertex at the subtraction point = = = 
and there is no renormalization. This, however, is not true if one renormal- 
izes at another point //^ > 0, even though it remains true to lowest order in 
perturbation theory [MP78]. Note also that this argument would be invali- 
dated if the two-ghost-two-gluon scattering kernel T^^{k, q) had an infrared 
divergence [AFLE05]. 
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Figure 2.3: The Dyson-Schwinger equation of the ghost-gluon vertex. 



We shall show in Sec. 4.4.2 that our lattice results for the renormalization 
constant Zi do not support the existence of such a divergence in the particular 
MOM scheme = r 2 and = 0. This was also seen by others [CMM04] 
in a similar lattice study of quenched SU{2) gauge theory. Furthermore, 
semiperturbative calculations of Zi using either = r"^ = k"^ = or = 
r'^;k'^ = show that such a divergence is absent. The analytical study 
presented in [AFLE05] also shows that the dressing of the ghost-gluon vertex 
remains finite if all external momenta vanish. 

We also refer to the recent study [B"'"05b] where arguments in favor of a 
singular infrared behavior of the ghost-gluon vertex function is given. 

2.3 Criteria for confinement in finear 
covariant gauges 

Despite the success of QCD in describing strong interaction processes at high 
energies, there is still the unsolved and theoretically demanding problem 
what kind of mechanism confines QCD degrees of freedom, the quark and 
gluon states, and keeps them from being observed in the physical particle 
spectrum. In contrast to QED, where a one-to-one correspondence between 
basic fields and stable particles can be assumed, in QCD it cannot. Only 
hadronic amplitudes are physical and only color singlets can contribute to 
the physical state space of QCD. 

For QCD in covariant gauges there is the additional subtlety that it re- 
quires a state space with indefinite metric. In order to allow for a quantum 
mechanical interpretation there has to be some mechanism which generates a 
subspace of (colorless) physical states that has a positive semi-definite metric. 

In the following three criteria for confinement are briefiy introduced, 
namely the Kugo-Ojima confinement scenario, the Gribov-Zwanziger hori- 
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zon condition and the violation of reflection positivity of propagators corre- 
sponding to confined particles. These criteria are proposed to be sufficient to 
indicate confinement; and we will see in succeeding chapters that our lattice 
simulations show evidence for those criteria to be satisfied in lattice QCD in 
Landau gauge. 

2.3.1 The Kugo— Ojima confinement scenario 

In this section we temporarily switch from the functional formalism to the 
covariant operator formalism which is more suitable in the present context. 
This formalism takes fields as operators, rather than as c-numbers, and they 
satisfy (anti-) commutation relations. Also the notion of BRST symmetry is 
important for a formulation of a covariant operator formalism. For a compre- 
hensive account of this subject we refer to the book [XO90] by Nakanishi 
and Ojima where most of the material summarized below is exhaustively 
discussed. 

Requirements for a physical S-matrix to exist 

Covariant quantum gauge theories require state spaces V with indefinite met- 
ric. If one can show, however, that (a) the Hamiltonian operator of a theory 
is hermitian so that a ^'-matrix S exists in V satisfying (5'\I'i|5'\E'2) = (^i|^2) 
and (b) there is a subspace Vphys Q V which not only is invariant under time 
evolution, but also (c) has a positive semi-definite inner product (-I-), then a 
physical ^'-matrix .Sphys can be defined in the completed quotient space^ 

"Wphys = Vphys /Vo • 

T^phys denotes a Hilbert space with positive definite inner product, and S'phys is 
unitary with respect to this Hilbert space structure. The subspace Vq C Vphys 
contains the zero-norm states of the positive semi-definite subspace Vphys and 
is orthogonal to it, i.e. Vo -L Vphys [NO90, A. 2]. 

It can be shown that the first two requirements, (a) and (b), are satisfied 
automatically for QCD in covariant gauges if the Lagrangian density /^gff 
(Eq. (1.27)) is hermitian'^ and the physical subspace Vphys is defined as the 
kernel of the BRST-charge Qb, i.e. 

Vphys = kerQe = {^r e V : Qs'if = 0} (2.11) 

^The overline denotes the completion of this space, i.e. all the limiting states of Cauchy 
sequences are incorporated in this, too. For a proof and further details see [NO90, A. 2]. 
''For this the ghost fields must satisfy c°'^ = and c"''' = c" [NO90]. 
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where is assumed to be an unbroken generator of BRST symmetry. 

The proof of condition (c), the positive semi-definiteness of the inner 
product in Vphys, however, is a nontrivial problem and requires a detailed 
analysis of the inner product structure of both the total state space V and 
the physical subspace Vphys [NO90]. 

Representations of the BRST algebra 

Such an analysis has been done by KuGO and Ojima [K079] long ago from 
the viewpoint of the BRST-algebra given by 



Here Qc denotes the FP-ghost charge which is assumed to be unbroken^ as 
it is assumed for Qb- 

The BRST-algebra has two types of irreducible representations, namely 
singlet and doublet representations. While a doublet consists of a so-called 
parent state \-k) and a daughter state \5) = QbItt) ^ (i.e. the BRST trans- 
form of Itt)), a BRST-singlet is a state that is annihilated by without 
having a corresponding parent state in V. 

KuGO and Ojima have shown that to each BRST-doublet there exists al- 
ways another being FP-conjugate to it, i.e. with opposite eigenvalues of iQc- 
This FP-conjugate pair of BRST-doublets is called a BRST-quartet. Any 
state in V can be classified to be either a BRST-singlet or to belong to a 
BRST-quartet and it can be shown that this exhausts all possible represen- 
tations of the BRST-algebra in spaces with indefinite inner product. We 
will see below that this classification is important, in that under certain con- 
ditions, colorless asymptotic states are BRST-singlets which then can be 
identified with physical particle states. On the other hand colored asymp- 
totic states are members of a BRST-quartet and, therefore, do not appear 
in S'-matrix elements. 

To see this note first that daughter states \6) and BRST-singlet states 
(with no parents) belong to the physical space Vphys as defined in Eq. (2.11). 
However, daughter states are orthogonal to all states ^ e Vphys 



and so cannot contribute to any element of the physical ^'-matrix -S'phys- All 
the physical content is contained in the physical Hilbert space 



{Qb,Qb} = 0, 



[iQcQb] = Qb, 



[QcQcl^o. 



= (^igslTT) = 



^phys = Vphys/ Vo ^ Vs. 



(2.12) 



^See also our discussion concerning the BRST-charge in Sec. 1.1.5. 
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which is isomorphic to the space of BRST-singlet states. Here 

Vo = img^ = m e V : = QbItt), {tt) e V} 
is the set of all zero-norm (daughter) states. 

Confinement of colored asymptotic fields 

Under the assumption of and Qc to be unbroken, the classification of 
representations of the BRST-algebra can be translated into the properties 
of creation and annihilation operators of asymptotic fields. Intuitively, these 
are understood to create asymptotic particle states which are observable in 
scattering experiments long before and long after collisions. 

Asymptotic states constitute two Fock spaces, Vin and Vout- By the postu- 
late of asymptotic completeness they are equal to the whole Hilbert space V, 
i.e. 

Vin = V = Vout. 

This is an important postulate and in particular it guarantees that any op- 
erator is expressible in terms of asymptotic fields [NO90]. For any two states 
\f),\g) &V an asymptotic field (p^ is defined as the weak limit 

lim {f\4>^^^-\x)-r'{x)\g) = 

of the corresponding renormalized operator 0^^^ = Z^^^'^cj). 

Due to confinement there must be some mechanism that causes colored 
asymptotic fields, if any. not to contribute to any physical iS-matrix element. 
Such a mechanism, the quartet mechanism, has been proposed by KUGO and 
Ojima [K079] for covariant gauge theories. 

They have analyzed the total state space V as the Fock space of asymp- 
totic fields and showed that any asymptotic field is either a BRST-singlet 
or a quartet member, because no other irreducible representation exists as 
they have shown (see above). If combined with the fact that under certain 
conditions (specified below) the charge of global gauge transformation is 
BRST-exact^ and can be written as (see next subsection) 

Q" = {gB,C} where C '.^ J d^x{Docy{x). (2.13) 

then for any physical states \ '^2) ^ ^phys it holds that 

^An BRST-exact operator ^ is a BRST variation Sb of another operator B, i.e. A is 
of the form: A = 5bB = {iQe, B}. 
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The charge vanishes in the physical Hilbert space Hphys defined in Eq. (2.12) 
Consequently, if there were colored asymptotic fields, they would belong to 
the BRST-quartet representations. BRST-singlets, on the other hand, are 
necessarily colorless and thus can be identified with physical particles. This 
is known as confinement by the quartet mechanism. 

One example is the so called elementary quartet that can be easily de- 
duced from the BRST transformation of and c" (Eq. (1.25) or (1.31)). 
Obviously, it consists of the parent states and |c") and of the daughter 
states \D'^^c'^) and IB"). As shown in [Kug9o] the corresponding (massless) 
asymptotic states — they describe longitudinally polarized giuons, ghost and 
antighosts — also form a BRST-quartet representation and are therefore not 
observable in the physical spectrum. 

It is also expected that the quartet mechanism applies to transverse gluon 
and quark states, as far as they exist asymptotically. A violation of (refiec- 
tion) positivity (see Sec. 2.3.3) for such states entails these are not observable 
either [NO90; AvSOl]. In Sec. 5.3 it is shown that the lattice gluon prop- 
agator in Landau gauge violates refiection positivity explicitly which thus 
supports this expectation. 

The Kugo-Ojima confinement criteria 

Color confinement by the quartet mechanism can only take place if there is 
an unbroken and BRST exact color charge Q". 

In general, Q"^ is a generator of the global gauge symmetry that, in ad- 
dition to the BRST symmetry, is left in the gauge-fixed effective Lagrangian 
density C^s- The corresponding symmetry transformations are given in 
Eq. (1.7) if there the parameters u)"' are taken to be space-time indepen- 
dent parameters. Being a global symmetry there exist Noether currents 
which are conserved , i.e. 

d,j;^0. (2.14) 

As pointed out first by Ojima [Oji78], these currents enter the equation 
of motion for the gauge fields in the form 

gr^ + d,F:^ = {QB,D,c'^}. (2.15) 

This equation is usually referred to as the quantum Maxwell equation in 
the non-abelian case, because for any physical states \ipi), \iJ2) G "^'phys the 
classical Maxwell-type equation {ipiKdi^F^^^ + gJ^)\'i,lj2) = holds. 

Since the conservation law in Eq. (2.14) allows adding arbitrary terms of 
the form d^f^^^^ to with /[^^] being a local antisymmetric tensor, one could 
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be tempted to define the global charge operators Q"" as the spatial integral 
of the current 

fa TO _| o rpa 

such that the BRST-exact expression in Eq. (2.13) is retrieved. But this 
naive definition of Q" is ill-defined due to massless one-particle contributions 
to J^^ix), dvF^v ^i^d {Qb, -Dy^iC"} which cause the integral to not converge 
[Kug95]. 

If however these contributions are consistently incorporated in the defi- 
nition of Q", a well-defined expression is obtained (see [Kug95] for details). 
This is important, since with the Goldstone theorem^ this automatically im- 
plies that condition (B) formulated in the original work [K079] of KUGO 
and Ojima, namely 

(B) is not spontaneously broken, 

is satisfied in the indefinite metric space V. If furthermore a certain pa- 
rameter u"**, the Kugo- Ojima confinement parameter, turns out to satisfy 
condition 

(A) = 

then the color charge takes the BRST-exact form as given in Eq. (2.13) 
and color confinement by the quartet mechanism 

(^i|g''|^2) = o ^i,^2GVphys 

takes place. But this holds only if condition (A) is fulfilled. Otherwise we 
cannot identify BRST-singlet states in the physical Hilbert space TCphys with 
color singlets. 

To investigate whether condition (A) is realized, the Kugo-Ojima confine- 
ment parameter u"** can be obtained as the zero-momentum limit 

u"^ := lim u^^p^) (2.16) 

of a function ?//''^(p^) which itself may be defined through the correlation 
function (see, for instance, [Kug95] or [AvSOl]) 

J d^xe'P^^'-y^ {Dl^c\x)gX'"^A'i{y)c\y)) =: (^5^^^ - u''\p^). (2.17) 

'''In Ref. [K079] various versions of the Goldstone theorem are given which altogether 
state that the following conditions concerning a conserved current and its global 
charge Q are equivalent: (1) Q = / d^xJo is a well-defined charge; (2) Q does not suf- 
fer from spontaneous symmetry breaking; (3) contains no discrete massless spectrum: 
(0|J^vl/(p2 = 0))=0. 
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To our knowledge a direct determination of u in terms of a lattice cal- 
culation of u"'^{p'^) has never been done. There are a few explorative lattice 
studies [NFOOb; NFOOa; FX04a:^FN04b], but these are based on data the 
ghost renormalization function ^3 (see below for the relation between Z3 and 
u.). The major problem is that in a lattice simulation can be calcu- 

lated only at finite momenta p and the data then have to be extrapolated to 
p = for which a suitable ansatz has to be chosen. 

The ghost propagator in the infrared is related to u{p) 

In Landau gauge the calculation of the correlation function in Eq. (2.17) can 
be even circumvented, because it has been shown [Kug95] that in this gauge 
the ghost dressing function J is related to m"^(p^) := 5"''^u{p^) according to 

Jip") = . ^ ( 2^^ 2 ( (2-18) 

and v{p'^) is an arbitrary function (see [Kug95] for a definition). In the zero- 
momentum limit this yields 

Therefore, if condition (A) is realized in QCD then the ghost dressing 
function in Landau gauge should diverge for vanishing momenta. Turning 
the argument around, if the ghost propagator is found to be more singular 
than a simple pole this may serve as a sufficient criterion for the Kugo-Ojima 
confinement scenario to be realized. 

Based on this arguments the infrared behavior of the ghost propagator as 
extracted from the corresponding truncated systems of DSE, has been stated 
to realize the Kugo-Ojima confinement scenario for QCD in Landau gauge. 
In support of this, in Sec. 4.3 we will show that also the ghost propagator as 
calculated in lattice simulations diverges stronger than a simple pole, albeit 
with different exponent compared to the one found in DSE studies. 

More importantly, however, in Sec. 5.2 we present data for the Kugo- 
Ojima confinement parameter u°'^{p) at finite momentum p and compare 
these data to those of the ghost propagator. The data indicate u°'^{p) to 
become proportional to —5"^ in the limit of vanishing momenta. So the 
Kugo-Ojima confinement scenario seems to be realized for QCD in Landau 
gauge. 



(2.19) 
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2.3.2 The Gribov-Zwanziger horizon condition 

In the previous section we have seen that the infrared behavior of the ghost 
propagator is intimately connected to one condition of the Kugo-Ojima con- 
finement criteria. In fact, if condition (A) is fullfilled then the ghost dressing 
function must be infrared divergent. 

In coincidence with this, there is another condition, namely the Gribov- 
Zwanziger horizon condition [Gri78; Zwa93], that also requires the ghost 
dressing function to diverge at vanishing momentum. In fact, the horizon con- 
dition states that the ghost propagator G in Landau gauge diverges stronger 
than l/p^ in the zero-momentum limit^, i.e. [Zwa02] 

lim [p^G{p)r^ = 0. (2.20) 

Basically, this limit is a consequence of the expectation that the infrared 
modes of the gauge fields are very close to the Gribov horizon^ dO, and hence 
give rise to an accumulation of small non-zero eigenvalues of the FP operator 
[Zwa94; Zwa93; ZwaQlb]. Since the ghost propagator, essentially, is the 
inverse of this operator, it must diverge in the infrared. 

Remarkably, the horizon condition as given in Eq. (2.20) allows us to 
renormalize the ghost propagator at p = in the form of a nonperturbative 
formula for the corresponding renormalization constant (see [Zwa04] for 
details). This formula disagrees with the usual perturbative expression for 
Z3, but it satisfies the perturbative renormalization-group flow equation. If 
used together with the DSE for the ghost propagator it even yields an infrared 
anomalous dimension kq for it, such that it behaves like [Zwa04] 



Gip) -^ K 



2\ «G 



p2 \^p2 

in the infrared. This is in agreement with the infrared behavior of the 
ghost and gluon propagators extracted from their Dyson-Schwinger equa- 
tions [vSAH97; AvSOl; LvS02]. 

Furthermore, it has been argued by ZwAXZlGER [ZwaQlb; Zwa92] that 
the gluon propagator in Landau gauge vanishes in the infrared, i.e. 

lim D{p) = 0, 



^But only, if the restriction to the Gribov region is done properly. 

^By definition, the Gribov horizon occurs where the lowest nontrivial eigenvalue of the 
FP operator vanishes. For typical configurations on large Euclidean volumes this operator 
is expected to have a high density of eigenvalues near zero [Zwa04]. 
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because the infrared components A{k) of the gluon field are suppressed by 
the proximity of the Gribov horizon in infrared directions [Zwa02]. 

In our study we will not only check for the infrared limits of both the 
gluon and ghost propagators as given above, but we will also show that the 
accumulation of near-to-zero eigenvalues of the FP operator increases with 
enlarging the physical volume (see Sec. 5.1, Sec. 5.2 and Sec. 6.2). Therefore, 
gauge configurations at the Gribov horizon seem to dominate the infrared 
properties of lattice Landau gauge theory in the thermodynamic limit. 



2.3.3 Violation of reflection positivity as a criterion for 
conflnement 

The mechanism for confinement by Kugo and Ojima introduced in Sec. 2.3.1 
relies on the existence of an unbroken BRST symmetry beyond perturbation 
theory. This, however, has not been proven yet and thus the Kugo-Ojima 
scenario should not necessarily apply to QCD. Nevertheless, the numerical 
results that we will discuss in Sec. 5.2, but also recent studies of truncated 
systems of Dyson-Schwinger equations for the ghost and gluon propagators 
(see next section) favor the Kugo-Ojima confinement scenario to be realized 
for QCD in Landau gauge. 

In any case, there is another particular criterion for confinement that has 
been proposed in recent years (see e.g. [AvSOl; ADFM04]) and is focused on 
in this study too, namely the violation of reflection positivity. In fact, refiec- 
tion positivity is an essential part of the famous Osterwalder-Schrader axioms 
[OS73; OS75] for Euclidean quantum field theory^". Arbitrary partial sums 
of Euclidean n-point functions have to fulfill those axioms to ensure their 
analytic continuation to the physically interesting functions in Minkowski 
space. They thus guarantee the reconstruction of a Carding- Wightman rel- 
ativistic quantum field theory. This is important in order to arrive at a 
physical interpretation. 

For the purpose of our study there is no need to go into detail about all the 
Osterwalder-Schrader axioms^^. We will rather focus on the notion of reflec- 
tion positivity which states that an Euclidean Green's function (Schwinger 



^"Reflection positivity of lattice gauge theory assures that gauge-invariant excitations 
have a physical spectrum [Zwa92]. 

^^The justification of this axioms is far from being trivial. For a comprehensive account 
on this the reader is referred to the books by Haag [Haa92] or Glimm and Jaffe [GJ87] 
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function) A has to satisfy 



i=l 



Xi 



r{@Xi,...,QXn) 



A{exi,...,ex„,yi,...,yn)f{yi,...,yn) > 0. 



(2.21) 



Here / refers to a complex valued test function^^ with support for positive 
(Euclidean) times, i.e. f{xi, . . . ,x„) = for any xf < 0; and © is the reflec- 
tion operator that acts on Xi according to: Qxi — (xj, —xf). 

In particular for a generic (Euclidean) 2-point fimction A{x, y) = A{x — y) , 
reflection positivity is a necessary and sufficient condition (see [AvSOl]) for 
the existence of a Kallen-Lehmann representation [Kal52; Leh54]. This is 
a spectral representation of 2-point functions^^ with positive, but generally 
unknown, spectral density p(m^) which in momentum space takes the form 
[AvSOl; A+97; CMT05] 



A(p) 



^0 



dm 



with p{m^) > 0. 



(2.22) 



The absence of a Kallen-Lehmann representation for a particular 2-point 
function is a sufficient condition for confinement of the corresponding particle, 
because then it cannot be interpreted in terms of stable particle states. 

Considering the temporal correlator C(t, p^), the absence of such a rep- 
resentation can be even formulated more straightforward. In fact, C{t, p^) is 
the Fourier transform of Eq. (2.22) and takes the form^^ [AvSOl] 



C{t, p') 



dp4 cos 



Jo 



dm 



2 I 9 



TT 



dm p{m ) — e 
2a; 



-Ult 



where uj"^ — m^ + p^. After substitution this gives 

/CO 
duj p{uj^ - p') e- 



■LOt 



(2.23) 



One clearly sees that if the spectral density p is a positive function then 
the temporal correlator C(t, p^) is positive as well, but it does not hold in 

^^Those functions belong to the Laurent Schwartz space of infinitely often differentiable 
functions, decreasing together with their derivatives faster than any power as x moves to 
infinity in any direction [Haa92]. 

^^This representation is very helpful for describing the analytic structure of propagators. 
Combined with the positivity requirements it yields bounds on their asymptotic behavior 
and the magnitude of rcnormalization constants. See e.g. [Wei95; PS95] for more details. 

^''Note, A(j3) is an even function of p4 which simplifies the Fourier transform. 
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general the other way. If on the contrary C{t, p^) is found to be negative for 
a certain range in t, i.e. 

C(i,p2)<0 (2.24) 

there cannot be a positive spectral density and thus reflection positivity is 
violated. This is an indication for confinement [AvSOl; CMT05]. 

For the gluon two-point function reflection-positivity violation, of course, 
is expected to happen. However, it has never be shown explicitly in lattice 
simulations for the case of SU{3). In three-dimensional pure SU{2) gauge 
theory numerical evidence for reflection-positivity violation of the lattice Lan- 
dau gluon propagators has been given [CMT05]. This study shows that also 
the SU{3) gluon propagator in Landau gauge violates reflection positivity 
for the quenched and unquenched case. 

We note in passing, that if the gluon propagator D{p) in momentum space 
would be infrared vanishing, as expected from the proximity of the Gribov 
horizon (see the previous section), then it would violate reflection positivity 
maximally. This can be seen from 

This can only happen if the gluon propagator in coordinate space, D{x), 
contains positive as well as negative contributions of equal integrated strength 
[ADFM04]. 

The ghost propagator violates reflection positivity trivially which can 
be deduced already from its bare expression. Hence ghosts are explicitly un- 
physical indeed. The unphysical spin-statistic relation of ghost fields already 
suggests this. 
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Chapter 3 



QCD Green's functions in 
LATTICE Landau gauge 



This chapter very briefly introduces the lattice regularization of QCD. We 
concentrate on the Wilson formulation of QCD with and without clover 
improved Wilson fermions as employed for this study. After specifying some general 
aspects we focus on lattice QCD in Landau gauge and define all relevant gauge- 
variant observables which are analyzed in the following chapters. 

3.1 Basics of lattice QCD 

We know that the very definition of a quantum field theory, as QCD, requires 
a regularization which often breaks some of the underlying symmetries of 
the classical theory and introduces a new scale into the theory [Has98]. At 
present the only known regularization of QCD beyond perturbation theory is 
the lattice regularization which discretizes Euclidean space-time into a lattice 
of points. Thereby, the lattice spacing a serves as a regulator of the theory 
that renders all ultraviolet divergences, usually encountered in QCD, finite. 
Unlike field theory with a naive ultraviolet cutoff, the lattice formulation 
maintains exact gauge invariance [Mut98; BHL"'"05]. Furthermore, the lattice 
regularization offers the possibility to investigate nonperturbative aspects of 
QCD in terms of numerical Monte Carlo (MC) simulations. Therefore, it is a 
valuable tool for cross-checking results obtained using other nonperturbative 
methods, for example the DS approach to QCD. 

It is beyond the scope of this thesis to adequately describe all aspects of 
lattice QCD. Therefore, we shall recall only some basic points important for 
the following. For a comprehensive introduction on the subject we refer to 
standard textbooks [Cre83; MM94; Rot97; Smi02] as well as to the lectures 
[Kog83; LepOO; Gup97; Dav02; Liis03; Dav05], to name but a few. Most 
of the material summarized below can be found there. To find out about 
the most recent developments in the field the proceedings [AHI+04; B"'"05a; 
IMM05] of the yearly Lattice Conferences provide a good starting point. 
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3.1.1 The lattice and its fields 

Very briefly, lattice QCD is a discretization of QCD in ^——^ ^ 
Euclidean space. It replaces the four-dimensional space- 
time continuum through a hypercubic lattice and restricts ^ z> ^ 
fermion and antifermion flelds, ipx and ipx respectively, to 
dwell on the lattice sites x. Rather than specifying the ^ ^ ^ A 
gauge flelds by gluon flelds on the lattice gauge U^^n 
flelds are associated with links joining adjacent lattice 
sites, X and x + fi. Here and in the following /i is a unit vector in the 
Xfj, direction of space-time. The lattice spacing a is the distance between ad- 
jacent lattice sites. The gauge flelds — also known as link variables or even 
just as links — are usually denoted by Ux^n and take values in a compact 
Lie group, here SU{3). They are the lattice version of the parallel transport 
matrix between adjacent sites and therefore are related to the continuum 
gluon flelds A^{x) by the line integral (see e.g. [LSWP99]) 

Ux,^, = P exp jigo ^ Af,{x + atfl) dt^ ~ ^^ag^A^{x+M2) _^ q (^^3^ _ ^3_;^^ 

Here P denotes path ordering of the gluon flelds along the integration path 
such that gluon flelds v4^(x + atfi) with larger t stand to the left of those 
with smaller t. The gauge flelds U are taken in place of the gluon flelds 
for the purpose of maintaining explicit gauge invariance on the flnite lattice 
[BHL"'"05]. A gauge invariant formulation of lattice QCD directly in terms 
of the gluon flelds is not possible [LepOO]. 

Quantization of lattice QCD is done in the functional integral formalism, 
i.e. expectation values of different observables are given in terms of path 
integrals over the gauge and quark held variables [MM94] 

(C) = ^ j [VU, Vi>, Vi)] 0[i), i), U] e-^ocD[c/,^,^] _ ^3_2) 

Here the partition function Z is chosen such that (1) = 1 and C[V^,V^,f/] 
denotes an arbitrary function of the held variables, for instance, n pairs of 
quark flelds and a product of link variables 

For simplicity, the indices Oj and hi denote the set of all internal symmetry 
indices, like color and spinor degrees of freedom. The color indices of link 
variables are hidden. Lattice sites are denoted by Xj, Ui or Zi and /ij refers to 
one particular direction on the lattice. 
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In contrast to any formulation in the continuum, the lattice action Sqcd 
maintains exact gauge invariance. Since the gauge fields take values only in 
the group S'C/(3), i.e. they are restricted to a compact manifold, the func- 
tional integration measure is well-defined and gauge-fixing is not necessary if 
gauge-invariant observables '?/^, L'^] are considered. However, in this the- 
sis gauge-variant quantities are investigated, in particular in Landau gauge. 
Therefore, gauge-fixing becomes necessary here as well. In Sec. 3.2 this is 
discussed in more detail. 

3.1.2 The Wilson action with clover-improved fermions 

Apart from the requirement of local gauge invariance the definition of Sqcd 
is not unique. There are infinite ways to define a lattice action of QCD, but 
any definition has to be such that it takes the classical continuum form for 
vanishing lattice spacing a. This freedom allows for a clever reduction of 
systematic errors caused by finite lattice spacings. Any reasonable formula- 
tion will give the same continuum theory up to finite renormalizations of the 
gauge coupling and the quark masses [Lus03]. 

Common actions used in the literature consist of two parts, a gauge Sq 
and a fermionic part Sp, i-e. they are of the general form 

SQCD = SG[U]+SF[U,i,,^P] . (3.3) 

Both parts depend on the gauge fields U = {U^^fj} and the fermionic part 
contains in addition bilinear expressions in the fermion fields ijj and tjj con- 
structed such that the whole action is manifestly gauge invariant. 

For the gauge part we have employed the standard Wilson gauge action 
[Wil74] throughout this study. It is given by the sum 



(3.4) 



X l<fj,<i/<4 

over traces of plaquettes denoted here by 

U,,,U,+/,,mU,^^UI, . (3.5) 

It represents a square of four links on the lattice. The parameter /3 is defined 
such that Sg takes its classical continuum form in the limit a — > 0, i.e. it is 
defined as 

27V 

P--=-^ (3.6) 
where = 3 for SU{3) and go is the bare coupling constant. 
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For the fermionic part Sp there are various definitions in use which ei- 
ther belong to the family of Wilson or of Staggered fermion actions or start 
therefrom. For the purpose of this study fermions of the Wilson type have 
been employed. In fact, our choice are clover-improved Wilson fermions 
[SW85; LSSW96]. The corresponding part of the action can be written as 

Sp[U,ij,ij]^a'Y,i^lQxyi^i (3.7) 

f,x,y 

where for each flavor / = 1, . . . ,Nf the fermion matrix Q is defined to act 
upon the (Grassmann valued) fermion fields, ipx and ipy, according to 



\ CI Qj 
X K 11 



ixw ) 

(see e.g. [GHP+05]). The fermion fields are normalized such that they cor- 
respond to the continuum fields by rescaling if) — > \I\/2ki\). Due to the last 
term in Eq. (3.8), Sp is suflicient to remove all O (a) errors for on-shell quan- 
tities^ The value of the parameter Csw depends on and has to be tuned 
appropriately. The clover field-strength tensor is given by 

^x,^T o- 2 E i^^'t^" ~ ^l,f^u) ' (3-9) 

where the definition of the plaquette (see Eq. (3.5)) has been extended such 
that the 11, v directions can be negative [GHP+05]. The hopping parameter 
K is related to the (subtracted) bare quark mass via 



2a \K 

where Kc is defined as the value of k at which the pion mass vanishes. In this 
study we consider only the cases of either infinite heavy quarks, i.e. no quarks 
{N f = 0) or Nf = 2 mass degenerate quarks. The former case is known as 
the quenched approximation of QCD, whereas the latter is an approximation 
to the real, unquenched world of two light and a number of heavier quarks 
that do respond to the gauge field. 

^For gauge dependent quantities it is an open question whether further, gauge non- 
invariant (but BRST invariant) terms must be added [SWOl]. However, one usually as- 
sumes any such term to be small. 
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3.1.3 Vacuum expectation values from MC simulations 

After having defined the lattice action we now come back to the definition 
of vacuum expectation values in Eq. (3.2) and outline the way Monte Carlo 
(MC) simulations are employed to calculate them. First note that the integral 
over the fermionic variables can be done instantly due to their Grassmann 
nature. We obtain [MM94] 

= 1 1 [DU] e-^e«[c/] gp^^ ^^^Z Q^'vM ■ ■ ■ QZvM (3.11) 

where the efi^ective action 5'eff[^7] (see Eq. (3.12)) depends only on gauge 
fields. Q\U] denotes an arbitrary function depending on link variables only. 
The tensor e^^^' ;^" := 1 (e^jy;;.'^" := — 1) if 2^1, . . . , is an even (odd) permu- 
tation of the lattice sites Xi, . . . , ,x„; or zero else. This tensor multiplied with 
elements of the inverse fermion matrix Q^^ appears due to the integration 
over (Grassmann valued) fermion fields present in the observable. It is absent 
if vacuum expectation values of pure gluonic observables ^[C/] are considered. 

Since the fermionic part Sp is usually of the form as given in Eq. (3.7), 
the effective action can be written as 

^eff[^] = Sg[U] - logdetg[t/] (3.12) 

where in our case the fermion matrix Q is defined in Eq. (3.8). The deter- 
minant in Eq. (3.12) is known as the fermion determinant. In the quenched 
approach to QCD this determinant is set to equal one. This simplifies numer- 
ical calculations of expectation values enormously, however, at the expense 
neglecting quark loops. 

To evaluate the remaining integral in Eq. (3.11) we can perform numerical 
MC simulations which allow us to estimate vacuum expectation values as 
statistical averages. Note that after discretization the generating functional 
Z of our lattice theory corresponds to a partition function of a statistical 
system. 

In a typical MC simulation of lattice QCD, sets of gauge field configu- 
rations U^^^ , C/^^^ . . . are successively generated by an appropriately chosen 
Markov process. Ideally, this process samples a large number of, say N , 
independent configurations being a realization of the Boltzmann weight 



-exp{-5eff[^]} 
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This is known as important sampling. On each configuration t/W the observ- 
able of interest O is measured such that the sample average 

1 ^ 

i=l 

is an estimator of the ensemble average {{0)u)- The latter is equal to the 
expectation value (O). The difference between the ensemble average and its 
estimator {0)u is important, because the latter is just an average over a finite 
sample of N gauge fields and therefore is naturally afflicted with a statistical 
error of aa/^/N. Only in the limit N ^ oo they would match. Here 
refers to the (usually unknown) standard deviation of the observable O. It 
can be estimated for example by the jackknife (e.g. [Wu86]) or the bootstrap 
method [ET93]. 

In addition to the statistical error an estimate is also afflicted by sys- 
tematic errors due to the finite volume V and the finite lattice spacing a, 
respectively. To estimate those, MC simulations are usually performed on 
different lattice sizes and at different setups of the parameter /3 (and eventu- 
ally k). The latter two are related to the lattice spacing and the quark mass. 
With respect to all these effect, ideally the continuum value is obtained by 
taking the multiple limit 

lim lim lim {0)u ■ 

V-^oo a^O JV— »oo 
y=const. 

This limit, in this order, corresponds to the prescription of an axiomatic field 
theory. 

3.2 The Landau gauge on the lattice 

The major focus of this thesis is to investigate the infrared behavior of the 
gluon and ghost propagators and other observables in Landau gauge. For 
their calculation gauge-fixing is necessary as well. In the following we discuss 
how the Landau gauge condition is imposed in our lattice simulations and 
then we define all observables relevant for this study. 

Gauge-fixing in the continuum usually includes a parameter, the gauge 
parameter ^, that causes the measure of the functional integral to peak 
around a particular gauge field on the gauge orbit. Lattice gauge-fixing 
achieves the same by a two-step process. First an ensemble of lattice gauge 
field configurations is generated using standard MC methods. Since the lat- 
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tice action S^s as defined in Eq. (3.12), is invariant under gauge transforma- 
tion 

^ '^U^,n = 9x C/x.M gl+fi, (3.13) 

tlie ensemble generated does not satisfy any gauge condition. Then in a sec- 
ond step, for each such generated configuration U = {Ux,^} a gauge trans- 
formation g — {g^} is chosen such that ^C4,^ satisfies the (lattice version of 
the) gauge condition. In this way, a particular configuration on the gauge 
orbit of Ux^fj, is chosen. 



3.2.1 The gauge functional 

For the particular case of Landau gauge, one usually searches for a gauge 
transformation g = {qx}, keeping U fixed, that maximizes a certain func- 
tional, the gauge functional Fu[g]. This typically reads 

X //=1 

where V denotes the lattice volume. Obviously, in the trivial case of U — 1 
the largest value Fu[g] = 3 is obtained for g — 1. For any other U, choosing 
a maximum of Fu[g] makes all on average as close to unity as possible. 
A continuum analog of the gauge functional has been given in Eq. (2.3). 

The functional Fu[g] has many different local maxima which can be 
reached by inequivalent gauge transformations g, the number of which in- 
creases with the lattice size. As the inverse coupling constant /3 is decreased, 
increasingly more of those maxima become accessible by an iterative gauge 
fixing process starting from a given (random) gauge transformation g. In 
this study we have employed two popular algorithms for gauge-fixing: over- 
relaxation [MO90a] and Fourier-accelerated gauge-fixing [D"'"88].^ The dif- 
ferent gauge copies corresponding to the maxima reached are called Gribov 
copies, due to their resemblance to the Gribov ambiguity in the continuum 
[Gri78]. All Gribov copies {^U} belong to the same gauge orbit spanned by 
the Monte Carlo configuration U. They all satisfy the differential Landau 
gauge condition (lattice transversality condition) (V^M^)(a;) = where 

4 

{V,A,){x) = (V • A){x) := J2 [M^ + A/2) - A,{x - /i/2)] . (3.15) 



^For a comparison of both algorithms see App. A. 2.1 
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Here An{x+fx/2) is the non-Abelian (hermitian) lattice gauge potential which 
may be defined at the midpoint of a link 

A,{x + /i/2) := 1 - C/t J - 1 IV - Ul^) . (3.16) 

In this way it is accurate to O {a?). Note that this association between lattice 
and continuum gauge fields is not unique, but the one chosen here represents 
the maximally local choice for such an assignment [M087]. The bare gauge 
coupling gQ is related to the inverse lattice coupling via /3 = 6/g^ in the case 
of 5?7(3) (see Eq. (3.6)). 

In the following, we will drop the label g for convenience, i.e. we assume 
U to satisfy the Landau gauge condition such that g = 1 maximizes the 
functional in Eq. (3.14) relative to the neighborhood of the identity. To 
simplify notation we will also use a more compact notation 

— + A/2) 

for the lattice gluon fields, but it is always understood that they dwell at the 
midpoint of a link.^ Additionally, we give the adjoint expression of a lattice 
gluon field 

Al^ := Alix + /i/2) = 2 • 3mT,{T^U,^} (3.17) 



3.2.2 The Faddeev- Popov operator 

Before we can go further and introduce the observables relevant for this study, 
it is necessary to give first a lattice expression for the Faddeev-Popov (FP) 
operator in Landau gauge. This operator can be easily derived by considering 
a one-parameter subgroup of the local SU{3) gauge group defined by (see 
e.g. [Zwa91a]) 

g^ (r, x) = exp {iTU^T"} r, a;^ e R . 

The generators T"^ of the SU{3) group have been defined in Sec. 1.1.1. In 

fact, if we assume U to represent a local maximum of the gauge functional 
then for any r it holds that Fu[l] > Fu[g^^{T)] for all u. Consequently, at 
r = the first derivative of the one-parameter function /u;(t) := Fu[gu){T)] 
with respect to r should vanish. One can easily show that 



- §^Ur) 



r-0 2 

T—IJ x,C fj, 



^This one should keep in mind if the Fourier transform of the gluon field has to be 
calculated. 
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and so any maximum of the gauge functional automatically satisfies the lat- 
tice Landau gauge condition. Then the second derivative of /(^(r) at r = 
defines a symmetric quadratic form 



9t2 



x,y,c,d 



whose kernel 



Mfy = S,,y - J2 5,+^,, + (3.18) 



with 



= E^^eIV[{r«,r^}(C/.,^ + C/._^,,)], (3.19a) 
= 2-lHcTr[r^r"C/,,J , (3.19b) 
C^l = 2-9leIV[r"r''C/,_^,^] . (3.19c) 

is the Hessian of Fi/[g]. M defines a real symmetric matrix that in the case 
of U satisfying V ■ A = equals the FP operator 

M[U] = -V • D[U] = -D[U] • V <^ V • A = 

where D[U] refers to the covariant derivative [Zwa94]. A lattice definition for 
D[U] can be found, for instance, in the same reference. After some algebra 
one can show that in the adjoint representation this can be written in the 
form 

(D^[U]y' = 2 Tr [T'T'^U,,^] S,+f,,y - 2 IHe Tr [T'^T'U,,^] 5,,y . (3.20) 

\ /XV 



3.2.3 Defining F, ft and A on the lattice 

On the lattice the Gribov ambiguity of the Landau gauge condition finds its 
expression in the ambiguity to find a local maxima of the gauge functional 
in Eq. (3.14). Therefore, terms like transversal plane, Gribov region and 
fundamental modular region also translate to the lattice formulation. 

The transversal plane is constituted by all (gauge transformed) configu- 
rations U that satisfy the lattice Landau gauge condition V • A{U) — using 
the definition (3.15), i.e. 



r:= {U:V-A{U) = 0} 
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The subset of F whose elements in addition give rise to a semipositive definite 
FP operator M is called the Gribov region 

n:^{U -.U e T,M[U] > 0}. 

Of course any element in H is a local maximum of the gauge functional, 
but only those which are global maxima constitute the fundamental modular 
region 

A:={U:Fuil)>Fu[g] for all . 

For a finite lattice it has been proven that the interior of A consists of non- 
degenerate absolute maxima (or minima depending on the definition). Gri- 
bov copies may only occur on the boundary dA [Zwa94]. 



3.3 Lattice definition of our observables 

3.3.1 The (inverse) FP operator in momentum space 

In subsequent sections we shall derive expressions for the numerical calcu- 
lation of the ghost propagator and the ghost-gluon-vertex renormalization 
constant in momentum space. For theses purposes the following Fourier 
transform 

{M-r'ik) = l,J2^-"'-iM-Xy (3-21) 

x,y 

of the inverse FP operator is of interest. Here and in the following the 
scalar product 

k.x = J2^7T^ (3.22) 

of lattice momentum k and lattice site x is understood. denotes the lattice 
extension in direction /x. 

Due to its eight trivial zero eigenvalues the inverse needs to be 

defined with care. If we are just interested in non-zero momenta k we are 
automatically in a subspace orthogonal to the space spanned by the (space- 
time constant) zero modes. Hence, for non-zero momenta we can apply, for 
example, the conjugate gradient method to solve the sparse linear system 

[M^,r = E ^-,a.v^r = er(^) (3-23) 

a,x 

using a fixed source ^5 with 8V complex components ^^^{k) := d'^^e^^'^ . Here 
c and z label the vector components oi^i,, while index h specifies the different 
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sources on the right hand side of Eq. (3.23), i.e. which of the color components 
are non-zero. The solution iph to this linear system can then be used to write 
the Fourier transform in Eq. (3.21) as the following scalar product in space- 
time^ 

{M-'f{k) = ^^e-'=--V^r(^) ■ (3.24) 

X 

With Eq. (3.23) it is clear that i/)^'"' represents the 8V" vector components 
with respect to the matrix multiplication of M"^ with ^b, i.e. 

V'r(^) = E^-!^.^^''' (^>o)- 

y 

For the numerical calculation of we rather solve the two independent 
linear systems 

[Mcbik)Y' = (3.25) 
[Msb{k)Y' = 6'*sm{k-z), (3.26) 

than that given in Eq. (3.23), because ■06^ = c^^ + is^^. With this notation 
we can write the Fourier transform in Eq. (3.24) as 

{M-^y'' {k) = ^J2 • ^)cr(^) + sin(A; • a;)s^^(A;) 

+i [cos{k ■ x)sr(fc) - sin(A; ■ x)cr(A;)] . (3.27) 

We shall see subsequently that the calculation of both, and s^, is even 
not always necessary depending on the observable considered. For example, 
for the calculation the ghost-gluon-vertex renormalization constant we only 
have to solve Eq. (3.26). This, as we shall see later, relies on the fact that 
the FP operator is symmetric and thus 



X x,y 



Ecos(A;-y)M„^;^^sin(A;-x) 

y,x 



{M=M'^) 



J]sin(A;-x)M^^;„^cos(A;-y) 

y,x 

J2cos{k-x) ■st^'ik) = J]]sin(A;-x) •c^^(A;). (3.28) 



^For our convenience, ipb carries the index b as well, in order to trace back afterwards 
the color index of non-zero components of ^(,. 
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3.3.2 The gluon and ghost propagator 
The gluon propagator 

Studying nonperturbatively gauge-dependent quantities on the lattice the 
gluon propagator is perhaps the simplest object to start with. Given the 
definition for the lattice gluon fields = + (1/2) in Eq. (3.16), the 
gluon propagator is estimated in lattice simulations by the MC average of 
the corresponding two-point function 

In this study we are in particular interested in the Fourier transform of this 
two-point function which on the lattice is given by 

Note that the term A'^^^Ay ^^ in general is not translational invariant. How- 
ever, we have imposed this invariance for the vacuum expectation value by 
summing over all differences {x—y) available on the lattice. If we assume that 
also on the lattice the gluon propagator in Landau gauge has the continuum 
tensor structure 

DfM = 5''' {^5''' - ^) D{q') (3.30) 
then all the physical information is contained in the scalar function 

D{<f) = J^^Y.D7M) (3-31) 

where A^c = 3 denotes the number of colors of SU{Z). 

Connecting a lattice momentum to its continuum counterpart 

Before we proceed with the ghost propagator we can define already here how 
the lattice momentum k is related to its continuum counterpart q. In fact, it 
is well-know that due to lattice artifacts the tree-level gluon propagator 
does not simply reproduce its continuum expression, but rather has the form 

- - '-^) (3.32) 
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where q^ik) is defined as 



2 

= — sin 

a 




) 



(3.33) 



Therefore, one usually employes Eq. (3.33) for translating a lattice momen- 
tum to the corresponding continuum momentum q^. Even more impor- 
tantly, in Ref. [LSWP99] it has been verified that the lattice gluon propaga- 
tor in Landau gauge shows scaling at (3 — 6.0 and 6.2 in the entire range of 
studied momenta if these are defined according to Eq. (3.33). Also sys- 
tematic effects with respect to the tensor structure of the gluon propagator 
(see Eq. (3.30)) are reduced. So the definition in Eq. (3.33) is a reasonable 
definition of q^ though it becomes worse at larger k^. Since we are interested 
mainly in the infrared properties of the gluon propagator we will use this 
definition of in our study. 

The ghost propagator 

Beside the gluon propagator we are also interested in the Landau gauge ghost 
propagator. Given the lattice definition of the FP operator in Eq. (3.18) this 
propagator (in momentum space) can be estimated in lattice simulations by 



Although the inverse of the FP operator itself is not translational invariant, 
the vacuum expectation value (here the ensemble average) has to be. There- 
fore, a sum over all possible differences {x — y) at the same momentum k is 
taken here again. 

Since the continuum ghost propagator is of the form G°''^{q) = 5"''^G{(f) 
we are interested in the scalar function 



^(^'(^)) = 1^E^""(^'(^)) = l^<TrA1-^(^)>^ (3.34) 



where in the last step the expression {M.~^Y"-{k) as defined Eq. (3.21) has 
been used. With respect to the discussion in Sec. 3.3.1, the trace can be 
calculated by using Eq. (3.27) and (3.28) which finally yields 



where d^^ and s'^^ are solutions to the two independent linear systems given 
in Eq. (3.25) and (3.26). 




TrA^-^(A;) = ^ cos(A; • x) ■ c^^ + sin(A; • x) ■ s^^ 



.ax 
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In our lattice simulations we solved the two linear systems by applying 
the pre-conditioned CG algorithm (PCG) where as pre-conditioning matrix 
we used the inverse Laplacian operator with diagonal color substructure. 
This significantly has reduced the amount of computing time as it is discussed 
in more detail in App. A. 3. 

Eigenmode expansion of the ghost propagator 

For a better understanding of the infrared behavior it is interesting to analyze 
the ghost propagator also by exploiting the spectral representation of the 
inverse FP operator for a given gauge field U in terms of its real (ascendent) 
eigenvalues Aj and its (normalized) eigenvectors 4>i{x) in coordinate space 

[M-\U)]:l^J2^t{x)UKy)- (3.35) 

i=i 

Here 0"(x) are the components of <Pi{x). Taking the Fourier transformed 

vectors ^i{k) at lattice momentum A;^j and averaging over a Monte Carlo 
(MC) generated ensemble of gauge field configurations one can compute the 
ghost propagator from truncated mode expansions 

Gr,{q'{k)) = {G{k\n))u (3.36) 

where 

1 " 1 

G{k\n) = - 5]-*,(fc)-*,(-fc) (3.37) 

denotes the contribution of the eigenvalues and eigenmodes on a given gauge 
field configuration. Here the vector and scalar product notation refers to the 
color indices. The Fourier momenta k^ are related to the physical momenta 
q^{k^) by Eq. (3.33). 

In fact, if the whole eigenvalue spectrum and all eigenvectors were known 
the ghost propagator would be determined completely. However, this is nu- 
merical too demanding. Nevertheless, restricting the sum in Eq. (3.37) to 
the n lowest eigenvalues and eigenvectors {n N = 8V — 8), we can figure 
out to what extent this sum saturates the full value G'(g^) determined using 
Eq. (3.34). We will see in Chapt. 6 that at lowest momentum the low-lying 
part of the FP spectrum gives the major contribution to G. 
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3.3.3 The Kugo-Ojima confinement parameter 

The next observable we are interested in is the Kugo-Ojima confinement 
parameter. According to Sec. 2.3.1 this parameter = ^"^(O) is defined 
as the zero-momentum Umit of some function u°'^{p'^) which itself has been 
introduced in Eq. (2.17) using a particular correlation function. On the 
lattice this requisite correlation function is of the form 

x,y c,d,e 

x,y c,d,e / ^ 

Here the (adjoint) lattice covariant derivative Z?^ acts upon a color-space 
vector as given in Eq. (3.20) and the lattice gluon fields are defined 
according to Eq. (3.17). With respect to the Lorentz-structure as given in 
Eq. (2.17) the function u"'^{q'^) is obtained by the sum 



over Lorentz indices where gg denotes the bare coupling constant and = 4: 
is the number of dimensions. For the actual lattice calculation of ii"^(g^) we 
note that the correlation function W^^ik) can be written as a scalar product 
of a plane wave and a vector -i/^b^j, multiplied by the covariant derivative 



Here ■06^,^ is the solution of the linear system Miph,^{k) — (f)b,v{k) with M 
being the FP matrix and a source (t>b,v{k) defined by the components 



ik-y 



As for the ghost propagator the pre-conditioned conjugate-gradient algorithm 
(see appendix A. 3) has been employed to extract ^^^^ for each non-zero mo- 
mentum k, separately. Then by using Eq. (3.38) the function u"''^ {(f {k)) is 
determined. 
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3.3.4 Renormalization of propagators 

Quantities like giuon or ghost propagators as they come out from typical 
lattice simulations have to be renormalized yet. Assuming multiplicative 
renormalization the renormalized continuum propagators Dr and Gr are 
related to the bare, dimensionless lattice propagators via 

a^D{a^q^) = Z^{ii^ ,a^)DR{q^] ij^) and 

The renormalization constant and Z^ are determined by imposing a renor- 
malization condition at some chosen renormalization scale [i^ . In this study 
we apply the MOM scheme for renormalization according to which the renor- 
malized propagators equal there tree-level form at some momentum jj?. Here 
for the scalar functions of gluon and ghost propagators it holds that 

Dn{q,iJi) =\ (3.39) 

GR{q,fi) , , = A- (3.40) 

In the subsequent chapter we consider in the majority of cases not the 
propagators themselves but their dressing functions Z and J (see Eq. (1.33) 
and (1.34)). They describe the deviation of the propagators from their tree- 
level forms. With Eq. (3.39) and (3.40) it is clear that at the renormalization 
point /x^ the dressing functions equal one. 



3.3.5 The ghost-gluon- vertex renormalization constant 

Apart from the renormalization constants of the gluon and ghost propagators 
we are also interested in Zi, the renormalization constant of the ghost-gluon 
vertex. In Landau gauge, the most general tensor structure of this ver- 
tex with gluon momentum s and ghost momenta q and t is given by (see 
e.g. [SMWA05]) 

T'^J'^is; q, t) = zgo [Qu ir'^ + A-'^ q\ t')) + s^B^"^ (s^; q\ f)] . 

Here A"^" and B"^" are scalar functions which describe the deviation from 
the tree-level from. They are assumed to have the same color structure as in 
perturbation theory, i.e. A""^^ =: f'^^Ai-s^; q^,t^) and B^^^ =: f'^^Bis^; q^,t^). 

In a MOM scheme the renormalized vertex Vr = ZiV equals its tree-level 
expression f at a renormalization point /x^. Therefore, if we consider the 
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particular renormalization scheme MOM which is defined by subtracting the 
vertex function at the asymmetric point — (f — /i^ and — 0, then the 
renormalization constant Zi is defined by 



where the tree-level expression is given by 

rf'^(0;g,t)=^gort.. 

Since in this renormalization scheme the tensor structure of the ghost- 
gluon vertex boils down to 

rf'=(0;g) = zgort.r(0;a 

where r(0;g^) = 1 + A(0;g^,g^) we arrive at the (inverse) vertex renormal- 
ization constant given by [CMM04] 

On the lattice, the definition for Z^^ can be derived in a similar manner 
as in the continuum (see e.g. [CMM04]). Only the tree-level expression of 
the vertex is different. Using the tree-level expression known from lattice 
perturbation theory the definition of the ghost-gluon- vertex renormalization 
constant is given by [CMM04] 

- ^E^-- (^) ^^"'^ :Jtnr»;^(0,,(^)) (3.41) 

where the constant c := 2/ {gQNc{N^ — 1)) and quiK) is defined in Eq. (3.33). 
Here L^, refers to the number of lattice points in direction u and ki, takes 
values in the interval (L,^/2, Lj,/2] as usual. 

The vertex Vf" can be obtained by amputating the external ghost and 
gluon legs from the three-point function of gluon, ghost and anti-ghost fields. 
In the MOM scheme considered here this yields [CMM04] 

fabc/r, \ 

where D and G refer to gluon and ghost propagators, respectively, and G"^*^ 
is given on the lattice by the MC average 

(0,g(fc)) = (^^(0) {M-'Y' (k)) ^ . (3.43) 
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In this expression .4"(0) denotes the Fourier transform of the gluon fields 
^4^ at zero momentum, i.e. 

X 

with A'^j^ defined in Eq. (3.16). Ai~^ is the inverse of the FP operator in 
momentum space considered at non-zero momenta (see Eq. (3.21)). If com- 
bined with the definition in Eq. (3.41) the renormalization constant Zj~^(g^) 
can be estimated by the MC average 

^- '('^) - T^mwmV'^" if) ■ 

where c := 2/{gQNc{N^ — 1)) and Gj^ represents the average 

GMk)) = E/"'' ^^G:'"^i^,q) = /j2A:{o)-r{k)) . (3.45) 

a,b,c \ a / u 

^"(0) is the same as defined above and 4>"'{k) is given by 

b,c x,y 

= E {cos{k ■ x) ■ s^^ - sin(A; • x) ■ c^^) 

6,c x,y 

= E 2/"^^ cos(A; • x) ■ s^^. (3.46) 

In the derivation we have used the antisymmetry of Z"^"^, Eq. (3.27) and 

(3.28). We see that if we calculate the ghost propagator anyway we can use 
the set of solutions Sc (c = 1,...,8), obtained at an intermediate step, to 
extract the ghost-gluon- vertex renormalization constant Z~^{q^) at the same 
momenta as the ghost propagator. 



Chapter 4 



Results for lattice QCD 
Green's functions 



In this chapter results for lattice QCD in Landau gauge, both in the quenched 
and unquenched cases, are presented. We start with a discussion of systematic 
effects due to finite volumes, discretization and the problem of Gribov copies. For 
this we restrict our attention to the quenched approximation. The infrared behavior 
of the gluon and ghost dressing functions is analyzed then, along with a discussion 
of unquenching effects. Subsequently, we report on results for the running coupling 
constant and the renormalization constant of the ghost-gluon vertex. 



4.1 General prerequisites 

4.1.1 Specification of our lattice samples 

In this study we have analyzed pure SU{3) gauge configurations, all ther- 
malized with the standard Wilson gauge action at three values of the inverse 
coupling constant P — 5.8, 6.0 and 6.2. The different lattice sizes studied are 
given in Table 4.1. For thermalization an update cycle of one heatbath and 
four micro-canonical over-relaxation steps was used. 

In addition, we have analyzed dynamical SU (3) gauge configurations pro- 
vided to us by the QCDSF collaboration^ Those configurations were gener- 
ated using the same gauge action, supplemented with the interaction with 
two flavors of clover-improved Wilson fermions. A deflnition of that action 
was given in Eq. (3.12). The different pairs of couplings are specifled 

in Table 4.2 together with the lattice sizes that have been used. 

^We thank the QCDSF collaboration, in particular Gerrit Schierholz and Dirk 
P LETTER, for giving us access to their configurations via the International Lattice Data 
Grid (ILDG). 
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4.1.2 Gauge-fixing 

Each gauge configuration U was transformed into Landau gauge by search- 
ing for a local gauge transformation g = {qx} that maximizes the functional 
Fu[g] defined in Eq. (3.14). For this purpose either one of the two algo- 
rithm commonly used, namely the over-relaxation (RLX) [MO90a] method 
or Fourier-accelerated gauge-fi,xing (FAG) [D+88] were employed. Both al- 
gorithms are iterative in nature and each iteration-cycle increases the func- 
tional Fu[g] until a (local) maximum is reached^. As stopping criterion not 
the functional itself, but the violation of transversality (see Eq. (3.15)) was 
used. In fact, the iteration process stopped as soon as 

max^leTr [(V/^,,^)(V/^,,^)T < e := 10"^" (4.1) 

was fulfilled at each lattice site, i.e. the lattice average was even lower. For 
some configuration we used e — 10~^^ which is also appropriate. 

4.1.3 The fc-bc strategy 

A subset of our quenched gauge configurations (see Table 4.1) were gauge- 
fixed more than once in order to investigate the infiuence of the Gribov 
ambiguity of gauge-dependent observables. In [SIMPSOSd] we have investi- 
gated such a dependence, following the strategy of choosing for each gauge 
configuration U the first (f c) and the best (be) gauge copy among N^p copies. 
Each copy has been fixed to Landau gauge always starting from a new ran- 
dom gauge copy of U. As best we have considered that copy with largest 
functional value among the Ncp gauge copies. Of course, the first gauge copy 
is as good as any other arbitrarily selected gauge copy. 

For each set of f c and be copies we then have measured the ghost and 
gluon propagators as well as the eigenvalue spectrum of the FP operator (see 
Chapt. 6) to determine the systematic effect caused by the Gribov ambiguity. 
In the following we call this particular way of studying the dependence on 
Gribov copies as the fc-bc strategy. Below we shall also motivate why we 
think that our values for A^^^p given in Table 4.1 are sufficient for this purpose. 
Of course, the more gauge copies one gets to inspect, the bigger the likeliness 
that the copy labeled as be actually represents the absolute maximum of the 
functional in Eq. (3.14). But as it is discussed in more detail below, the ex- 
pectation value of gauge variant quantities, evaluated on be representatives, 
is converging more or less rapidly with increasing number N^p. 

^In App. A. 2.1 wc compare both algorithms and analyze how the iteration numbers 
scale with the lattice size. 



4.1 General prerequisites 



73 



no. 




lattice 


[GeVl 

I J 


a [fml 

L J 


#conf 


cp 


S-1 


5.8 


16^ 


1.446 


0.1364 


40 


30 


S-2 




24^ 






25 


40 


F-1 




24^ 






40 


30 


S-3 


5.8 


32^ 


1.446 


0.1364 


34 


1 


S-4 


6.0 


16^ 


2.118 


0.0932 


40 


30 


S-5 




24^ 






30 


40 


S-6 




32^ 






40 


1 


S-7 




48^ 






20 


1 


A-1 




24=^ X 48 






30 


1 


A- 2 




32^ X 64 






40 


1 


A-3 




163x 128 






30 


1 


A-4 


6.0 


24=^ X 128 


2.118 


0.0932 


30 


1 


F-2 


6.2 


12^ 


2.914 


0.0677 


150 


20 


F-3 




16^ 






100 


30 


S-8 




16^ 






40 


30 


F-4 




24^ 






35 


30 


S-9 




24^ 


2.914 


0.0677 


30 


40 



Table 4.1: The /3 values and lattice sizes used in simulations of the quenched case. Also 

numbers of configurations used arc given. iVcp specifies the number of different random 
gauge copies considered for each configuration. The 5th (and 4th) row lists the (inverse) 
lattice spacings corresponding to (3. Labels in the first row are used in the text. 



no. 


P 


K 




ma 


a [fm] 


[GeV] 


#conf 


D-1 


5.29 


0.13500 


0.13641(9) 


0.03828 


0.0957 


2.063 


90 


D-2 


5.29 


0.13550 


0.13641(9) 


0.02462 


0.0898 


2.196 


60 


D-3 


5.29 


0.13590 


0.13641(9) 


0.01376 


0.0850 


2.320 


55 


D-4 


5.25 


0.13575 


0.13625(7) 


0.01352 


0.0904 


2.183 


60 



Table 4.2: The f3 and k values of all dynamical gauge configurations used in this study. 
The Kc values are taken from Ref. [G"'"06] where also the values for the Sommer scale in 
lattice units ro/a are specified. The latter were used to assign physical units to a. We also 
give values for ma = l/2(l//c — 1/kc)- The lattice size is 16"^ x 32 for the first (D-1) and 
24^ X 48 for the other three sets (D-2, D-3, D-4). The numbers of configurations used are 
given in the last row. For all sets A^cp = 1- Labels in the first row are used in the text. 



74 Chapter 4 Results for lattice QCD Green's functions 



4.1.4 Selection of momenta 

Our implementation of the lattice gluon propagator first constructs the lat- 
tice gauge fields A^^^ according to Eq. (3.17) and then Fourier-transforms 
these using a Fast-Fourier transformation (FFT) algorithm^. Because a FFT 
provides us with all lattice momenta at once, the data for the gluon prop- 
agator D{q^{k)) have been determined for all momenta available. 

It is obvious from Eq. (3.33) that different lattice momenta give rise to 
the same value q^(k). Naively, one would average over all data of D{q^{k)) 
at different /c^ but same (fik). This however leads to systematic errors due 
to finite volume and discretization effects. It has been shown [LSWP99] that 
both of these systematic errors can be reduced by applying two cuts on the 
data, i.e. only a subset of momenta is used. But within this subset, all data 
with same q^{k) are then averaged over all different realizations of k^ and 
configurations. 

One of these cuts, known as the cylinder cut [LSWP99], reduces errors 
due to finite lattice spacings. Conceived in general terms, it selects data of 
D{k) = D(q^{k)) with k lying in a cylinder with radius of one momentum 
unit along one of the (lattice) diagonals n — 1/2(±1, ±1, ±1, ±1). To be 
specific, 



./i=0 



L,, 



4 I. 



E 



./i=0 1^ 



< ^ (4.2) 



where Ls is the extension in spatial direction. In the special case of Lt — Lg, 
Eq. (4.2) reduces to A;^ - {J^^k^n^f < 1. In agreement with [LSWP99] 
this recipe has drastically reduced lattice artifacts for the gluon propagator, 
in particular for larger momenta. 

The other of the two cuts is known as the cone cut [LSWP99]. It addresses 
finite volume errors by removing all data D(k) with one or more vanishing 
momentum components A;^. We have applied this cut to our data only for 
smaller lattice volumes. In the next section we shall show in more detail 
how finite volume efi'ects have influenced our data at lower momenta. In 
particular, if asymmetric lattice geometries are used the cone cut is necessary. 

For the ghost propagator and some other quantities considered in this 
thesis, the numerical calculation involves an individual inversion of the FP 
operator for each vector k. Due to limited computing time we did not had 
the chance to obtain data for all momenta allowed by the cuts. However, our 
selection of different k were guided by either cuts. 



^For all FFTs we have employed the FFTW-hbrary [FJ98], see also online: 
http : / / www . f f tw . org. 
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4.1.5 Mapping to physical units 

The P values chosen for the quenched case allow us to apply the results of 
Ref. [NS02], giving a parameterization of the functional dependence of the 
lattice spacing a, or better of ln(a/ro), on {3. Using such parameterization and 
the Sommer scale tq = 0.5 fm [Som94] we obtain for the three values (3 — 5.8, 
6.0 and 6.2 a'^ = 1.446 GeV, 2.118 GeV and 2.914 GeV or a = 0.1364 fm, 
0.09315 fm and 0.0677 fm, respectively. In our opinion, this mapping between 
a and /? is more appropriate as another one formerly used by us and others 
(see [SIMPSOSa; SIMPS05b; SO04; LSWP99]). 

For our sets of unquenched gauge configurations the values of ro/a were 
provided to us by the QCDSF collaboration (see Table II in Ref. [G"'"06]). 
Using again ro = 0.5 fm we can assign to each pair of (5 and k a lattice 
spacing in physical units. These are given in Table 4.2 together with other 
specifications. 

Finally, with the help of Eq. (3.33) we can then map lattice momenta A;, 
or better a?q^{k), to physical momenta (f. If not otherwise stated, (f is 
always given in GeV^. 

4.2 Systematic effects on gluon and ghost 
propagators at low momentum 

In this section different systematic effects on the gluon and ghost propagators 
are discussed. We start with effects due to finite volume and finite lattice 
spacings, followed by a warning to refrain from using quite asymmetric lattice 
geometries. Finally, attention is paid to the dependence on Gribov copies. 



4.2.1 Finite volume and discretization effects 

Precisely because we have applied the above-mentioned cuts to our data, 
it is quite natural to analyze here the different systematic effects on the 
gluon and ghost propagators of changing either the lattice spacing a or the 
physical volume V . However, due to the preselected set of momenta for the 
ghost propagator and the three chosen (3 values, our study is partial and 
limited to a region of intermediate momenta. For the gluon propagator this 
has been done in more detail by other authors (see e.g. [BBL+01]). 
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Finite volume eflfects 

Keeping first the lattice spacing fixed we have found that both the ghost 
and gluon dressing functions calculated at the same physical momentum 
(f decrease as the lattice size is increased. This is illustrated for various 
momenta in Fig. 4.1. There both dressings functions versus the physical 
momentum are shown for different symmetric lattice sizes at /? = 5.8. 6.0 
and 6.2. In contrast to our study [SIMPS05d] here we show data obtained 
on f c gauge copies, because the larger lattice sizes were manageable only 
without repeating gauge-fixing several times. 

In this figure we have not dropped data with vanishing momentum com- 
ponents (i.e. those excluded by the cone cut) to emphasize the influence 
of a finite volume on those (low) momenta. We also show data from simula- 
tions on a 8^ and 12^ lattice. One clearly sees that the lower the momenta 
the larger the effect due to the finite volume. In comparison with (3 — 5.8 
and 6.0 this is even more drastic aX (5 — 6.2. At this (5 the lattice spacing is 
about a = 0.06 fm. Thus the largest volume considered at /? = 6.2 is about 
(1.4 fm)^, which is even smaller than the physical volume of a 16^ lattice at 
13 = 5.8. 

In summary we can state that for both dressing functions finite volume 
effects are clearly visible at volumes smaller than (2.2 fm)^, which corresponds 
to a 16^ lattice at (3 — 5.8. The effect grows with decreasing momentum or 
decreasing lattice size (see the right panels in Fig. 4.1). At larger volumes, 
however, the data for g > 1 GeV coincide within errors for the different 
lattice sizes (left and middle panels). For q < 1 GeV we find only small finite 
volume effects for both dressing functions at the lowest momentum if data 
obtained on a 24^, 32^ and a 48"^ lattice at /3 = 5.8 and 6.0 are considered. 

Discretization errors 

Based on our chosen /3 values and lattice sizes we can pick up equal physi- 
cal volumes (with different coarseness) only approximately. Hence also the 
physical momenta are only approximately the same if the ghost and gluon 
dressing functions are compared at different /?, i.e. at different lattice spac- 
ings. Therefore, it is difficult to analyze the systematic effect of changing a if 
for both dressing functions this leads to small hidden variations in q^. Conse- 
quently, in Fig. 4.2 we show data for the ghost and gluon dressing functions 
obtained at approximately the same physical volume V ~ (2.2 fm)^ for two 
different a as functions of g^. This allows us to disentangle by inspection 
a change of data due to varying a, given the physical dependence of the 
propagators on q^. Looking at Fig. 4.2 one concludes that the gluon dressing 
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Figure 4.1: The ghost (upper panels) and gluon (lower panels) dressing functions for 
different lattice sizes as functions of momentum . Prom left to right the panels show 
data at /3 = 5.8, 6.0 and 6.2. Only data on f c gauge copies are shown here. Lines are 
drawn to guide the eye. 



function at the same physical momentum and volume increases with decreas- 
ing the lattice spacing. A similar effect (beyond error bars) is not observable 
for the ghost dressing function. 



4.2.2 Asymmetric lattices cause strong systematic 
errors 

Besides discretization and finite volume errors, there are also systematic ef- 
fects involved using asymmetric lattice geometries. It could be tempting 
to use an asymmetric lattice such that in one direction, for example in 
time directions, the lattice extension is much longer than in all other di- 
rections. Going this way, one would tend to believe that much lower mo- 
menta can be studied compared to using symmetric lattice geometries at 
same computational costs. This approach has been pursued for example in 
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Figure 4.2: The ghost (left) and gluon (right) dressing functions are shown as functions 
of momentum at (approximately) fixed physical volume V k, (2.2 fm)**. The data at 



(3 = 5.8 (6.0) correspond to a lattice spacing of about a 
obtained on f c gauge copies are shown here. 



0.136 fm (0.093 fm). Only data 



[SO05b; OS05a; SO05a; OS05b] aiming at the infrared behavior of the gluon 
propagator, in particular, at the determination of an infrared exponent for 
this propagator. 

However, as a comparison of our data obtained either for symmetric or 
asymmetric lattice geometries suggests: "There is no free lunch!". The more 
asymmetrically a lattice has been chosen, the larger are the systematic errors 
encountered by that. To make this point more clear, we have calculated both 
the ghost and gluon propagators not only on symmetric lattices as discussed 
above, but also on lattices of sizes 32^ x 64, 24^ x 48, 16^ x 128 and 24^ x 128 
all at the same P = 6.0 (see runs A-1 to A-4 in Table 4.1). The corresponding 
data are plotted in Fig. 4.3 together with those obtained on a 32"^ and a 48^ 
lattice at the same value of /? (S-6, S-7 in Table 4.1). The lower panels of 
this figure show the corresponding dressing functions. 

In Fig. 4.3 we clearly see that data obtained on a 16'^ x 128 lattice suffer 
under the largest systematic errors in the low momentum region compared 
to the other data. In fact, the data on a 32'^ x 64 lattice, shown using open 
squares in this figure, exhibit the smallest asymmetry effects. There only 
the data point at the lowest on-axis momentum for the ghost propagator 
is lower in tendency than an extrapolation including the lowest momentum 
data point referring to a 48^ lattice would give. For the gluon propagator 
this deviation is even smaller. Also for data on a 24^ x 48 lattice, shown 
using open diamonds, the asymmetry effect is small. There only the data 
point at the lowest on-axis momentum k = (0, 0, 0, 1) is lower (larger) than 
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Figure 4.3: The ghost G and gluon propagators D (upper panels) and the corresponding 
dressing functions (lower panels) are shown sX (3 — 6.0 as functions of the momentum in 
lattice units. The data were obtained using different lattice geometries. Lines connecting 
data points for symmetric lattice geometries are drawn to guide the eye. 



the corresponding result for the ghost (gluon) propagator on a 48"^ lattice 
at the same momentum a?q^{k). At larger momenta this effect is negligible 
within errors for these particular sets of data (24^ x 48 and 32'^ x 64 lattice). 

Considering instead data generated on the lattices sizes 16'^ x 128 and 
24^ X 128, the results at the six or seven lowest momenta are significantly 
different to what is expected from the data obtained on symmetric lattices. 
A data inspection yields that the distorted momenta are only the on-axis 
momenta k = (0, 0, 0, kt) allowed by a cylinder cut, but forbidden by a cone 
cut (see Sec. 4.1.4 for a definition of these cuts). 

Similar systematic effects due to asymmetric lattice geometries have been 
reported recently for the gluon propagator in three-dimensional pure SU (2) 
gauge theory [CM06], but in fact already in the well-known reference [LSWP99] 
for pure SU{3) gauge theory. Therefore, extractions of an infrared exponent 
for the gluon propagator using only on-axis momenta on asymmetric lattices 
(see e.g. the references mentioned above) should be taken with reservations. 

We note in passing that we also see in Fig. 4.3 points from asymmetric 
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lattice geometries at not so small momenta which are slightly larger (lower) 
compared to those for the ghost (gluon) propagator on symmetric lattices. 
With respect to our discussion of finite volume effects above, we can, however, 
interpret this observation more generally as a finite volume effect caused by 
the lower spatial volumes of these asymmetric lattices. 

4.2.3 Facing the problem of Gribov copies 

Apart from the systematic effects inherent in any lattice study, it is known 
that the presence of Gribov copies may also causes a systematic error on the 
data of gauge-variant quantities. Below it is shown that this is indeed the 
case for the ghost propagator, while the infiuence on the gluon propagator 
seems to be negligible (i.e. hidden in the statistical noise). In [SIMPS05d] we 
have investigated this following the f c-bc strategy introduced in Sec. 4.1.3. 
For each set of f c and be copies we have measured the ghost and gluon 
propagator to determine the systematic error caused through the Gribov 
ambiguity. 

Estimating the number of gauge copies 

In pursuing the f c-bc strategy, at some stage one has to figure out how large 
Ncp has to be approximately. This does not mean that we can be confi- 
dent that the set of best copies then in any respect represents the real ones, 
i.e. those which give rise to the absolute maximum of the gauge functional 
Fu[g] (Eq. (3.14)) for each U. It turns out, however, that for each gauge- 
variant observable there exists a certain Ncp that warrants the convergence 
of that particular observable. 

Numerically, it turns out that the dependence of the ghost propagator 
on the choice of the best copy is most severe for the smallest momentum. 
In addition, this sensitivity depends on the lattice size and (3. Therefore we 
studied first the dependence of the ghost and gluon propagators at lowest 
momentum (ignoring the cone cut) on the (same) best copies as function of 
the number of gauge copies N^p under inspection. This was done at /3 = 6.2 
where we used 12'\ IG'' and 24^ lattices. The number of thermalized config- 
urations used for these three lattice sizes are given in [Table 4.1: F-2 to F-4]. 
To check the dependence on (3 also a simulation at (3 — 5.8 on a 24^ lattice 
was performed [Table 4.1: F-1]. 

The results of this investigation are shown in Fig. 4.4 and 4.5. While 
there the ghost propagator is shown as an average over the two realizations 
k = (1, 0, 0, 0) and k = (0, 1, 0, 0) of the smallest lattice momentum a^(f{k), 
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Figure 4.4: The upper panels show the ghost propagator G{k) as average over two 
reaHzations k — (1, 0, 0, 0) and k = (0, 1, 0, 0) of the smallest lattice momentum, measured 
always on the best gauge copy among N^p copies. In the middle panel the same dependence 
is shown for the gluon D{k) propagator, however, as average over all four permutations 
of fc = (1,0,0,0). The lower panels show the relative difference 5F = \ — F'^^'^/F^'^ of 
the corresponding current best functional values F'^^'^ to the value F^"^ of the overall best 
copy. The data are obtained at (i — 6.2 using the lattice sizes 12^ (left) and 16^ (right). 
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Figure 4.5: The same as in Fig. 4.4, however, the data refer to /3 = 5.8 and 6.2 on a 24'' 
lattice. 
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Figure 4.6: The ratios J^^/J*"" and Z^'^/Z^'' of ghost (l.h.s.) and gluon (r.h.s.) dressing 
functions, respectively, vs. momentum . J*"^ and Z*"^ refer to values determined on first 
(f c) gauge copies, whereas J*"^ and Z^'^ to those on best (be) copies. The ratios were 
calculated using jackknife. 



the gluon propagator has been averaged over all four non-equivalent realiza- 
tions. Note that D{k) = D{-k). 

In Fig. 4.4 and 4.5 it is clearly visible that the expectation value of 
the gluon propagator does not change within errors as N^p increases, in- 
dependently from the lattice size and /3. Contrarily, the ghost propagator 
at /3 = 5.8 on a 24^ lattice saturates (on average) if calculated on the best 
among N^p = 15 gauge copies. At /3 = 6.2 the number of necessary gauge 
fixings attempts reduces to 5 < Ncp < 10 on a 16^ and 24^ lattice. On the 
12"^ lattice a small impact of Gribov copies is visible, namely 1 < A'cp < 5 
is sufficient for convergence. The lower panels of Fig. 4.4 and 4.5 show the 
average of the relative difference 6F = 1 — F<=^'=/F^<^ of the corresponding 
(current best) functional value F'^^'^ to the value F^'^ of the overall best copy 
after A^'cp = 20, respectively N^p = 30, attempts. This may serve as an indi- 
cator how large A^cp has to be, on average, for the chosen algorithm to find a 
maximum of F close to the global one. In Table 4.1 the respective numbers 
Ncp for the different simulations are given. 

Gribov copies may cause systematic errors 

Having specified the necessary amount of gauge copies we focus now on the 
influence the Gribov ambiguity might have on the gluon and ghost propaga- 
tors. To investigate this a combined study of the gluon and ghost propagators 
on the same sets of f c and be representatives of our thermalized gauge field 
configurations has been performed (see [SIMPS05d]). This has allowed us to 
assess the importance of the Gribov copy problem for the ghost propagators 
in the low-momentum region. 
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To be specific we have generated different sets of quenched gauge con- 
figuration using the inverse coupling constants (5 — 5.8, 6.0 and 6.2 and the 
lattice sizes 16^ and 24^. The corresponding entries in Table 4.1 are S-1, 
S-2, S-4, S-5, S-8 and S-9. Following our f c-bc strategy, each gauge con- 
figuration U has been fixed to Landau gauge a number of N^.^ times using 
over-relaxation and starting always from an random gauge copy of U. Then 
for each U the first (f c) and the best (be) gauge copy (among A^cp copies) 
form the two ensembles for measurements of the ghost and gluon propagator. 

In Fig. 4.6 we illustrate the effect of the different Gribov copies by plotting 
the ratios J^" j 3"°" and Z^" jZ^"^ of the ghost and gluon dressing functions, 
respectively. J^"^ and Z'^" refer to values determined on first (f c) gauge copies, 
whereas J*"" and Z^"^ to those on best (be) copies. Obviously, in Fig. 4.6 there 
is no influence visible for the gluon propagator within the statistical noise. On 
the contrary, for the ghost propagator the Gribov problem can cause O (5%) 
deviations in the low-momentum region ((^ < 1 GeV) . For better gauge copies 
the ghost dressing function becomes less singular in the infrared. This can 
also be seen in Fig. 4.7 and 4.8 where both f c and be data for the ghost 
propagator are shown. 

A closer inspection of the data in Fig. 4.6 and 4.8 indicates that the influ- 
ence of Gribov copies on the ghost propagator becomes weaker for increasing 
the lattice size. In fact, comparing in Fig. 4.6 ratios for the ghost dressing 
function at g < 1 GeV, the rise at /? = 6.0 is obviously larger than that at 
(5 = 5.8. In both cases the data are from simulations on a 24"^ lattice. Thus, 
it seems that by increasing the physical volume (lower /?) the effect of the 
Gribov ambiguity gets smaller if the same physical momentum is considered. 
This we can also deduce from Fig. 4.7 and 4.8. There obviously the differ- 
ence between J^'' and J^" at lower momentum reduces when /3 decreases. We 
think this observation is not biased by a too small number A^cp of inspected 
gauge copies since, judging from Fig. 4.4 and 4.5, A^cp = 40 as specified in 
Table 4.1 seems to be on the safe side. However, in a future project one 
should definitely check the ratios on a 32^ lattice at /3 = 6.2 to eliminate the 
last doubts. 

We conclude: the ghost propagator systematically depends on the choice 
of Gribov copies, while the impact on the gluon propagator is not resolvable 
within our statistics. However, there are indications that the dependence 
on Gribov copies decreases with increasing physical volume. This is also in 
agreement with the data listed in the two lattice studies [BIMMP04; Cuc97] 
of the SU{2) ghost propagator G, while it is not explicitly stated there. In 
fact, in Ref. [BIMMP04] the ratio G^^/G°^ at /? = 2.2 on a 8^ lattice is larger 
than that on a 16"^ lattice at the same physical momentum. In [BBMPM05] 
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Figure 4.7: The ghost dressing function determined on first (f c) and best (be) gauge 
copies vs. momentum . Only data on a 24" lattice are shown here at three different 
values of (3. Lines connecting data points belonging to the same (3 are drawn to guide the 
eye. 
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Figure 4.8: The same as in Fig. 4.7, but focussing on a lower range of momenta. In 
addition, f c data on a 32" lattice at /3 = 5.8 and on a 48" lattice at /? = 6.0 are shown. 
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similar indications have been found for SU (2) taking non-periodic Z{2) trans- 
formations into account. 

Note also that this observation is in agreement with a recent claim by 
Zwanziger according to which in the infinite volume limit averaging over 
gauge configurations in the Gribov region should lead to the same result as 
averaging over configurations in the fundamental modular region [Zwa04]. 



4.3 The infrared behavior of gluon and ghost 
propagator in the quenched and 
unquenched case 

After discussing different systematic effects on the gluon and ghost prop- 
agators or on the corresponding dressing functions we concentrate now on 
the infrared behavior. Thereby, however, we shall study a further system- 
atic effect, namely the change of data at lower momenta due to unquenching 
the theory. As we shall see below, this change is negligible for the ghost 
propagator, but conspicuously large for the gluon propagator at not so small 
momenta. 



4.3.1 The gluon propagator 

From the material presented above it is clear that the gluon propagator in 
Landau gauge does not suffer from the Gribov ambiguity, at least in the mo- 
mentum range considered. Therefore, if we restrict ourselves in the following 
to consider only f c gauge copies, the systematic effect involved will stay most 
probably within statistical errors. This has the advantage of allowing us to 
use data obtained on larger lattice sizes even though they have been collected 
only from the ensemble of f c gauge copies. 

Just as a reminder, in our simulations the gluon propagator was deter- 
mined for all momenta at once. However, only a subset of data was used for 
the final analysis. For details see Sec. 4.1.4 where our selection of momenta 
is specified. The way how lattice momenta are connected to physical ones 
is discussed in Sec. 4.1.5. Beside this, it is also clear that all data obtained 
in lattice simulations have to be subject to a renormalization, discussed in 
Sec. 3.3.4. In the following the data presented for the gluon and ghost dress- 
ing functions, Z and J, are renormalized such that at the renormalization 
point = 4 GeV they are fixed to Z {/j,"^ , /j,"^) = J{ii^,ii^) = 1. 
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The quenched gluon propagator 

In Fig. 4.9 the dressing function of the gluon propagator is collected from 
data at different (3 values as a function of the (physical) momentum q^. Based 
on the experiences gained in previous sections only data obtained on larger 
lattice sizes are presented there. Beside of a cylinder cut applied to all data, 
a cone cut has been imposed to data related only to the lattice sizes 16'^ 
and 24'^. Since finite volume errors on the larger lattices are small, it is not 
necessary to impose this cut there. 

Obviously, the data surviving all cuts lie along a smooth curve which not 
only is nonperturbatively enhanced around q — 1 GeV as expected, but also 
decreases in a sufficiently large range of low momenta. Therefore, it appears 
natural trying to fit the ansatz 

fD{q')^An-{qY'' (4.3) 

to the data at lower momenta. As mentioned above, from studies of truncated 
DSE for the gluon propagator an infrared exponent Kd = 2k, with k ~ 0.595 
is expected to describe the infrared behavior [LvS02; Zwa02]. 

We have tried to fit this ansatz to the data by imposing a upper-momentum 
cut < qf that has been varied. The parameters Ad and kd extracted and 
the values of qf used are listed in Table 4.3. As can be seen in this table 
shifting qf to larger momenta a fit with ansatz (4.3) becomes worse. The 
parameters of our best fit to the data have been given in bold letters and in 
Fig. 4.9 we have plotted the corresponding fitting function foiff)- Although 
this fit supports the conjecture of an infrared vanishing gluon propagator our 
value of kd = 0.83(2) is lower than expected from the DSE studies. However, 
there is a tendency for hd to rise as the interval of momenta, contributing to 
the fit, is shrinking towards the infrared. Therefore, we cannot exclude the 
expected infrared behavior to become realized at much lower momenta. 

The unquenched gluon propagator 

Beside of our measurements in the quenched approximation we have also 
determined the gluon propagator on gauge-fixed copies (gauge-fixed only 
once) of dynamical gauge configurations specified in Table 4.2. The results 
of the corresponding dressing functions are shown in Fig. 4.10. With respect 
to finite volume effects reported above for the quenched case we have plotted 
in this figure only data obtained on a 24^ x 48 lattice, nevertheless for three 
different pairs of (3 and k (see D-2, D-3 and D-4 in Table 4.2). Note, the data 
related to the 16"^ x 32 lattice (D-1) turn out to exhibit similar finite volume 
effects as reported for quenched configurations. 
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Figure 4.9: The dressing function of the gluon propagator renormalized at ^ = 4 GeV 
is shown as a function of momentum using various lattice sizes. All data have been 
obtained from first gauge copies in quenched simulations (see Table 4.1). The line refers to 
an infrared fit on the data as explained in the text. The momentum qf marks the largest 
momentum used for this fit. 




Figure 4.10: The same as in Fig. 4.9, including now data obtained on first gauge copies of 
unquenched configurations (see D-2 to D-4 in Table 4.2). For comparison some quenched 
data at /3 = 6.0 [Table 4.1: S-6, S-7] have been included in this figure as well. A dotted 
line connecting quenched data is drawn to guide the eye. 
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In order to perform a comparison with our quenched data with the great- 
est of ease some of the quenched data have been included in this figure as 
well, namely the results at /3 = 6.0 associated with a 32^ and a 48'* lattice. 
The lattice spacing associated with this is comparable to those of the dy- 
namical configurations. As above the dressing functions were renormalized 
at // = 4 GeV. 

As can be seen from Fig. 4.10, the unquenching effect becomes clearly 
visible for the gluon dressing function, in particular around ~ IGeV^. 
There the non-perturbative enhancement, characteristic to this function, is 
drastically reduced compared to the quenched data {ma = oo). It also be- 
comes softer as the quark mass is decreasing, even though this is a small ef- 
fect. This has been observed also in recent lattice computations of the gluon 
propagator using configurations generated with dynamical AsqTad-improved 
staggered quarks [BHL"'"04] as well as from unquenching studies for the ghost 
and gluon propagators within the DSE approach [FAC+06; FWC05; FA03]. 
We refer also to recent lattice studies with dynamical Kogut-Susskind and 
Wilson fermions [FN05a; FN05b]. The difference between quenched and un- 
quenched data in the ultraviolet asymptotic tail is consistent with what is 
expected from perturbation theory (see Sec. 1.3.4). 

Unfortunately, the amount of data collected for momenta in the range 
< IGeV^ does not allow for a fit to the ansatz given in Eq. (4.3). In- 
tuitively, one may guess, looking at the low-momentum tendency of the 
unquenched data points in Fig. 4.10, that in the infrared the unquenched 
gluon dressing function might match that of the quenched case. 

4.3.2 The ghost propagator 

In the discussions of different systematic effect we have found that the Gribov 
ambiguity has an influence on the ghost propagator. Nevertheless, in Fig. 4.8 
we have also seen that even though the effect is visible in the data, the 
systematic error one would encounter by neglecting the Gribov copy problem 
is not that large. Since at present our data collected on f c gauge copies cover 
much lower momenta, we analyze in the following only those. Incidentally, 
this also allows for a fair comparison of our quenched and unquenched data. 

The quenched ghost propagator 

The ghost dressing function has been determined on the same gauge con- 
figurations as that of the gluon propagator discussed above. However, the 
list of momenta studied is much shorter than that of the gluon propagator. 
Anyway, we have tried to cover as much as possible the whole range of mo- 
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menta with special focus on the infrared region. The selection of momenta 
was guided by the cylinder cut and for smaller lattices sizes sharpened by 
the cone cut. 

Our results for the quenched ghost dressing function are shown in Fig. 4.11. 
As for the gluon dressing function we used = 4 GeV as renormalization 
point such that J(/i^,//.^) = 1. The data all lie along a smooth curve. Thus 
we may conclude that both cylinder and cone cut work well also for the ghost 
propagator. Note that in contrast to the gluon dressing function here a cone 
cut has not been imposed on a 24"^ lattice. 

As expected the ghost dressing function seems to diverge with decreasing 
momenta and we have tried to fit the infrared behavior again by an ansatz 
motivated through the mentioned studies of the ghost propagator within the 
DSE approach. To be specific, we fitted the ansatz 

(4.4) 

to the data at lower momenta. From the mentioned studies one usually 
expects to get kq ~ 0.595 [LvS02; Zwa02]. Even more importantly, one 
expects Kg = 2^^ to hold. See below for a discussion about that. 

The ansatz (4.4) has been fitted to our data by imposing an upper mo- 
mentum cut < qf. However, in contrast to the gluon dressing function, 
here the parameter Aq and kq turn out to be more robust against shifting 
gf. This can be seen in Table 4.3 where the parameters obtained are given. 
In Fig. 4.11 we show the corresponding fitting function. Even though data 
at the lowest three momenta available on a 24'^ lattice are shown in Fig. 4.11 
they do not contribute to the fit. We also have used data only at (3 = 6.0, 
because including data at /3 — 5.8 in the fit makes it even worse. 

Our fits to the data suggest that kq — 0.20(1) which is far away from 
values larger than 0.5. Therefore, we cannot confirm the expected infrared 
exponent n ^ 0.595 on the basis of our data. Also the not so small values 
found for x^/ndf indicate that the power ansatz (4.4) seems to be not reason- 
able for the range of momenta considered. Note also that the linear rise of 
the data in Fig. 4.11 (momentum log-axis) for decreasing momenta suggest 
the ghost dressing function might diverge logarithmically at least in the con- 
sidered range of momenta. However, we do not know arguments that would 
support this. 

Together with our results for the gluon dressing function we can neither 
confirm the expected values of the infrared exponents nor that there the 
relation 

kd — 2k,g = A 

holds with A = 0, at least from the data available to us. Our fits yield 
A f» 0.43. Furthermore, one should remark here that even though we have 
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Table 4.3: Parameter extracted corresponding to the fitting Ansatze (4.4) and (4.3) for 
the ghost and gluon dressing function, respectively. The first row specifies the upper- 
momentum cutoff rj^, i.e. < qf . Bold letters indicate parameters used for the fitting 
function shown in Fig. 4.11 and 4.9, respectively. 



fitted the Ansatze (4.4) and (4.3) to the corresponding dressing functions, the 
fits are quite unstable. For the ghost dressing function the fitting function 
does even not really describe the tendency of data towards the infrared (see 
Fig. 4.11). 

The unquenched ghost propagator 

In Fig. 4.12 we present our full QCD result for the ghost dressing functions. 
There we have not discard data at the lowest on-axis momentum in order to 
illustrate the systematic effect encountered through the asymmetric lattice 
geometry. For comparison, selected data of the quenched case (i.e. for infinite 
quark mass) are also shown, namely those on 32^ and 48^ lattice at (3 — 6.0. 
As in the quenched case the dressing function has been renormalized such 
that Jin'^.iJi'^) = I Hi 11 = A GeV. 

In contrast to the gluon dressing function, in Fig. 4.12 the unquenching 
effect is very small. In the ultraviolet the unquenched data are slightly lower 
and at lower momenta they are slightly larger than the quenched data, but 
altogether the effect stays within error bars and thus can be neglected in 
the region of momenta considered here. Negligible unquenching effect are in 
agreement with what is expected from unquenching studies for the ghost and 
gluon propagators within the DSE approach [FAC+06; FA03]. We refer also 
to recent lattice studies with dynamical Kogut-Susskind and Wilson fermions 
[FN05a; FNOob]. Thus, even though we cannot confirm the predicted infrared 
exponent (from the DSE studies), we affirm tliat the influence of fermions on 
the infrared behavior of the ghost propagator is negligible. 
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Figure 4.11: The dressing function of the ghost propagator renormaUzed at /i = 4 GeV 
is shown as a function of momentum using data from various lattice sizes. All data 
have been obtained from first gauge copies in quenched simulations (see Table 4.1). The 
line refers to an infrared fit to the data as explained in the text. The momentum gf marks 
the largest momentum used for this fit. 



3.0 



2.5 



— I — I I I I I I \ \ \ — I I I I I I 

ma) = (5.29, 0.0246) V-^ 
(5.29,0.0138) I — ^ 
(5.25,0.0135) r -e- 
quenched I 1— 



1.5 



2.0 



1.0 



0.1 



10 



100 



g2 [GcV2 



Figure 4.12: The same as in Fig. 4.11, however, most of the data are obtained on first 
gauge copies of unquenched configurations using a 24^ x 48 lattice (see D-2 to D-4 in 
Table 4.2). For comparison, some quenched data (3 — 6.0 [Table 4.1: S-6, S-7] have 
been included into this figure as well. A dotted line connecting quenched data is drawn to 
guide the eye. The unquenched data at lowest momentum refer to the lattice momentum 
k — (0,0,0, 1) on a 24^^ x 48 lattice and are most probably affected by lattice-asymmetry 
effects. 
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4.4 The running coupling and the ghost-gluon 
vertex 

After extensively reporting on results obtained for the gluon and ghost dress- 
ing functions, now we focus on the running coupling constant as{(f). As 
explained in Sec. 2.2.2, as(g^) may be defined by a renormalization-group- 
invariant combination of the gluon and ghost dressing functions (Eq. (2.10)) 



Remember that this definition relies on the assumption that the ghost-gluon 
vertex in Landau gauge stays bare also beyond perturbation theory. Below 
we give numerical evidence confirming Zi 1 in a MOM scheme where the 
gluon momentum equals zero. This we show for both the quenched and un- 
quenched case of SU{?)). Similar results indicating this, directly [CMM04] 
and indirectly [BCLM04], were presented in studies of the quenched SU{2) 
gauge theory, but also in semiperturbative calculations within the DSE ap- 
proach [SMWA05]. 

4.4.1 Results for the running coupling constant 

Based on our data for the renormalized gluon and ghost dressing functions, 
we have estimated the product in Eq. (4.5) using the bootstrap method with 
drawing 500 random samples. Since the ghost-gluon-vertex renormalization 
constant Zi has been set to one, there is an overall normalization factor which 
has been fixed by fitting the data for > to the well-known perturbative 
results of the running coupling a2-ioop at 2-loop order (see also [BCLM04]). 
Defining x = Q'^/Aj.ioop, the 2-loop running coupling is given by 



The /3-function coefficients /3o and Pi have been defined in Eq. (1.45a) and 
(1.45b). They are independent of the renormalization prescription. The 
value of A2-100P has been fixed by the same fit. The lower bound has been 
chosen such that an optimal value for x^/ndf has been achieved. 

The results are shown in Fig. 4.13. There also the 1-loop contribution 
is shown where we used the lower bound — 50 GeV^. The best fit of the 
2-loop expression to the data gives A2.100P = 1.15(15) GeV (x^/ndf = 7.5), 
while Ai_ioop = 0.75(30) GeV is obtained (x^/ndf = 5.2) using just the 1-loop 
part. For the 2-loop expression we used q^ = 30 GeV^. The value for A2_ioop 




(4.5) 




(4.6) 
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Figure 4.13: The running coupling as{q'^) as a function of the momentum . The data 
refer to first (f c) gauge copies of quenched configurations. 
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Figure 4.14: The running coupling as{q'^) as a function of momentum determined on 
gauge-fixed configurations in the unquenched case. For comparison, we also show quenched 
data (crosses) obtained at /3 = 6.0 on a 32'' and 48* lattice. 
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is similar within errors to the SU(2) result given in Ref. [BCLM04]. Note 
that we have imposed again a cone cut to data obtained on a 16^ and 24^ 
lattice. 

Approaching the infrared limit in Fig. 4.13 one clearly sees a.si'f) increas- 
ing for > 0.4 GeV^. However, after passing a maximum at ~ 0.4 GeV^ 
as{q^) decreases again. The same behavior is found on our sets of dynam- 
ical gauge configurations. These results are presented in Fig. 4.14 where 
quenched data {f3 — 6.0: 32"' and 48'') are also shown for comparison. In 
fact, looking at Fig. 4.14 we clearly see that the same turnover, as found in 
the quenched case, can be assumed for the unquenched case, even though 
the data at the lowest (on-axis) momenta have to be taken with special care. 
(See our discussion concerning systematic effects due to asymmetric lattices.) 
Therefore, the data points at ~ O.lGeV^ have to be checked carefully on 
larger symmetric lattices, in order to eliminate the last doubts. 

We observe a quite clear unquenching effect for Q!s(g^) extending from the 
perturbative range down to the infrared region. We think that this is caused 
in the majority due to unquenching effects as found for the gluon propagator 
in Sec. 4.3.1. 

Note that a similar infrared behavior of Q;s(g^) as presented here has 
been observed in other lattice studies [FN04a; FN04b]. But opposed to 
Ref. [FN04b] we argue that the existence of a turnover is independent on 
the choice of Gribov copies. In fact, qualitatively we have found the same 
behavior if Oigiq^) is calculated on be gauge copies. To illustrate this, in 
Fig. 4.15 we show the corresponding data extracted in be copies. There we 
have not impose the cone cut on data associated with the 24"^ lattice. Ob- 
viously, cxsiq'^) decreases also in this case and a decreasing running coupling 
constant is not due to restricting to f c gauge copies. The values for Ai.ioop 
and A2.100P are the same within error bars as obtained on f c gauge copies. 

We should remark here that our observation of an infrared decreasing 
running coupling is in agreement with recent studies of DSE results ob- 
tained on a torus [FAR02; FA02; FGA06; FP06]. In those studies a,(g2) 
was shown to tend to zero for — > in one-to-one correspondence with 
what one finds on the lattice. This would indicate very strong finite-size 
effects and a slow convergence to the infinite- volume limit. However, on the 
lattice we do not find any indication for such a strong finite-size effect, ex- 
cept the convergence to the infinite-volume limit would be extremely slow. 
An alternative resolution of this problem has been proposed by Boucaud 
et al. [B+05b]. These authors have argued^ that the ghost-gluon vertex in 
the infrared might contain g^-dependent contributions which could modify 



We thank A. Lokhov for bringing us the arguments in [B"'"05b] to our attention. 
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Figure 4.15: The running coupling as{q^) as a function of momentum determined on 
best (be) gauge copies of quenched configurations. 

the DSE results for the mentioned propagators. Note, however, that recent 
DSE studies of the ghost-gluon vertex did not provide hints for such a modi- 
fication [SMWA05; AFLE05]. Thus, at present there seems to be no solution 
of the puzzle. 

For completeness we mention that running coupling constants decreas- 
ing in the infrared have also been found in lattice studies of the 3-gluon 
vertex [BLM+98a; B+03] and the quark-gluon vertex [SK02]. 



4.4.2 The vertex renormalization constant 

For discussing the infrared behavior of the running coupling it is interesting 
to have an independent check of whether the renormalization constant Zi 
really stays constant for all renormalization points. Remember, this is always 
assumed if renormalization-group invariance is demonstrated for the product 
of gluon and ghost dressing functions in Eq. (4.5) and hence is important for 
this particular definition of a nonperturbative running coupling constant. 

In principle the finiteness of Zi has to be checked in each renormalization 
scheme separately. A first attempt to calculate Zi on the lattice was made 
in [CMM04] for quenched SU (2) gauge theory using the particular renormal- 
ization scheme of zero gluon momentum (see also [BCLM04] for an indirect 
determination). To confirm these findings also for the gauge group SU{?>) in 
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Figure 4.16: The (inverse of the) ghost- 
gluon- vertex renormahzation constant Z\ is 
shown versus the momentum . The data 
refer to quenched simulations at (5 = 5.8, 
and 6.0 using a 32^^ and 48"^ lattice. Only 
data on f c gauge copies are available. 



Figure 4.17: The same as in Fig. 4.16, but 
the data are obtained on unquenched gauge 
configurations using a 24^ x 48 lattice. The 
points at the lowest momentum refer to the 
on-axis momentum k = (0,0,0,1) on that 
lattice. Labels refer to the parameters k 
and P given in Table 4.2. 



the quenched and unquenched case, we have used the same renormalization 
scheme as in [CMM04] for a calculation of Zi. The necessary observable to 
be estimated on the lattice has been derived in Sec. 3.3.5 (see Eq. (3.44)). 
Since this is a combination of MC averages for the ghost and gluon propa- 
gators and for the 3-point function G^'"^ of gluon, ghost and anti-ghost fields 
(see Eq. (3.43) and (3.45)) we used again the Bootstrap method to estimate 
errors. 

Our data for the inverse of Zi from quenched simulations are shown in 
Fig. 4.16. To simplify matters, we have restricted ourselves to present only 
data obtained on the largest lattice available at P = 5.8 and 6.0. For data 
on smaller lattices we refer to our recent conference proceeding [SIMPS06]. 
The data there agree within errors with those presented here. 

In Fig. 4.16 we clearly see that Zi stays constant in the region of momenta 
considered. This holds for the data at both /? = 5.8 and /3 = 6.0. Only 
a slight variation is visible in the interval 0.3 GeV^ < Q'^ < 2 GeV^, but 
this remains within error bars. The same can be concluded from our data 
obtained on unquenched configurations using a 24^ x 48 lattice. These are 
shown in Fig. 4.17 for three different settings of n and P (see the data set 
entries D-2, D-3, D-4 in Table 4.2). For all three cases, the renormalization 
constant Zi does not differ beyond error bars from being constant, but there 
is a certain trend of deviation from unity which systematically depends on 
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the parameter setting {(3, k). 

For each tuple there is also one data point at the lowest (on-axis) mo- 
mentum available on a 24^ x 48 lattice (i.e. k = (0,0,0,1)) which is much 
lower than unity. This particular deviation we also find at the lowest on-axis 
momentum in our data obtained on quenched configurations at /3 = 6.0 when 
using a 24^ x 48 and 32^ x 64 lattice (not shown). Since the data taken in 
simulations on a 48^ lattice at the same (3 do not show such a deviation, 
see Fig. 4.16, this efi^ect is most probably caused by the asymmetric lattice 
geometry as it was reported also for the propagators in Sec. 4.2.1. Data 
at this momentum on an asymmetric lattice would have to be ignored. To 
eliminate the last doubts we have also inspected our data for obtained 
at the lowest on-axis momenta on a 24^ x 128 lattice aX (3 — 6.0. We find 
for these data that this particular deviation is even more dramatic. Hence, 
asymmetric lattices are not appropriate for studying the ghost-gluon-vertex 
renormalization constant at very low momentum. 

In summary, our results for the quenched and unquenched case of SU (3) 
are in full agreement with those presented in [CMM04] for the case oiSU{2). 
Even though there is a weak deviation of Zi from being constant this will not 
have a dramatic influence on the running coupling. Together with our data 
for the running coupling we can thus confirm that in the special MOM scheme 
considered here (gluon momentum equals zero) the product in Eq. (4.5) is 
indeed renormalization-group invariant and thus defines a nonperturbative 
running coupling which monotonously decreases with decreasing momentum 
in the range < 0.3 GeV^. _ 

It is worthwhile to continue the lattice calculations of Zi using other 
renormalization schemes, for example the scheme with a symmetric subtrac- 
tion point. This has been undertaken in a recent study [SMWA05] where a 
semiperturbative calculation of Zi within the DSE approach has been pre- 
sented. There are also data shown for the same renormalization scheme as 
used by us which qualitatively agree with our results and those in [CMM04]. 
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Chapter 5 



Confinement criteria under 
the lattice microscope 



Here we discuss data for the SU{3) Landau gauge gluon and ghost propagators 
in the Hght of the Gribov-Zwanziger horizon condition and the Kugo-Ojima 
confinement scenario. For the latter we also present recent data for the function 
u{q'^) whose zero-momentum limit u{q^ = 0), the Kugo-Ojima confinement param- 
eter, is expected to be minus one. Note that our lattice estimate of the function 
n(g^) is (to our knowledge) the first presented in the literature. The data have been 
renormalized using a minimization process, developed here for the first time. We 
show that, under the assumption the ghost dressing function being divergent in the 
infrared, the function n(g^) will reach minus one at zero momentum. Finally, we 
discuss numerical evidence for the gluon propagator violating reflection positivity 
explicitly. 



5.1 Is the Gribov-Zwanziger horizon condition 
satisfied? 

The Gribov-Zwanziger horizon condition has been introduced in Sec. 2.3.2. 
It states that on the one hand the ghost propagator in Landau gauge diverges 
in the limit of vanishing momenta more rapidly than whereas on the 

other hand the gluon propagator vanishes in the same limit. 

It should be clear from the discussion in the previous chapter that our 
data support the picture of a diverging ghost propagator, even though we 
cannot confirm an infrared exponent > 0.5, a value expected from studies 
of the ghost DS equation. Therefore, the Gribov-Zwanziger horizon condition 
appears to be satisfied with respect to our data of the ghost propagator in 
Landau gauge. 

Concerning the gluon propagator, how-ever, we cannot give a conclusive 
statement of whether it vanishes or stays finite in the infrared. We have 
tried to fit the ansatz given in Eq. (4.3) to our data for the gluon dressing 
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function. As mentioned in Sec. 4.3.1 this ansatz is reasonable only if applied 
to data at the lowest momenta available to us. We also could not confirm 
the corresponding value of the infrared exponent as expected from studies of 
the gluon DS equation. In any case, with respects to the range of momenta 
available to us we cannot state without doubt whether a power law as that 
given in Eq. (4.3) really describes the behavior of the gluon dressing function 
at much lower momenta. Therefore, we cannot judge on the existence of an 
infrared finite or vanishing gluon propagator, even though our fit with ansatz 
Eq. (4.3) supports the latter option. 

There have been attempts in the recent literature [BBL+01; B+06c] to 
argue for a finite gluon propagator in the infrared based on lattice data for 
the gluon propagator at zero four-momentum, i.e. 

r)(0) ^ \^ / /la Aa \ 

( ]\[2 _ l\ \ x,iJ. x+y,fj,/u- 

^ ' IJ,,a x,y 

Here ^4°^ refers to the lattice gluon field as defined in Eq. (3.17) and V 
represent the volume in lattice units. For notations we refer to Chapt. 3. In 
order to estimate the corresponding value -Doo(O) in the infinite volume limit, 
in the Refs. [BBL+01; B+06c] the ansatz [BBL+01] 

D{Q) = ^ + D^{0) (5.1) 

has been applied to the data of D{0) for different volumes^ To set the 
scale for -D(O), the gluon propagator D{(f,y?') has been renormalized in 
MOM scheme choosing as renormalization point ji = A GeV. In doing so, 
both studies [BBL+01; B"'"06c] independently confirm a manifest finite value, 
namely 

An extrapolation of our data yields the estimate 

5oo(0;Ai = 4 GeV) = 8.27(10) GeV'^ 

that approximately agrees with those of these studies. See also Fig. 5.1. Note 
that for c we obtain c = 102(18) fm^GeV^^ which agrees within errors with 
that in [B+06c], but not with [BBL+01]. In [BBL+01] the tree-level, mean- 
field improved gauge action of Liischer and Weisz [Wei83; WW84; LW85] 
has been employed, while in [B+06c] and in this thesis the standard Wilson 
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Figure 5.1: The gluon propagator at zero four-momentum, £'(0), plotted as a function 
of the inverse lattice volume. For comparison, data of the Adelaide group given in 
[BBL+01] have been included into this figure, too. The Hne refers to a fit of the linear 
ansatz given in Eq. (5.1) to the data at /3 = 5.8 and 6.0. 

gauge action has been used. This might be the reason for the deviation in c. 

It is important to note that estimates of D^o have to be taken with cau- 
tion, because these actually do not represent the true infinite volume limit 
i?oo(0). It is also not clear whether a linear ansatz is correct [BBL"'"01]. In 
this context, it is worthwhile to recall a comment already made in [BBL"'"01] 
concerning the infinite volume limit. This cannot be taken such simple as 
outlined above. In fact, the infinite volume limit is given by 

L>oo(0) = lim lim D{0). (5.2) 

V'=const. 

That is, the continuum limit (at fixed physical volume) has to be taken before 
the infinite volume limit. This was not done for the given estimates of Dqo and 
thus a complete systematic extrapolation to the infinite volume limit remains 
to be carried out yet [BBL+01]. For this one would first need to calculate 
D{0) using different lattice spacings a, while keeping the physical volume 
fixed. The results would then have to be extrapolated to the continuum limit. 
Repeating this for a variety of lattice volumes the different continuum limits 

^Note that we distinguish between -Doo(O) and DodO)- See below for a discussion. 
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obtained one could try to extrapolate to the infinite volume limit D^{0). 
This, of course, is the essence of the limit in Eq. (5.2). 

Therefore, we think a conclusive statement of whether in Landau gauge 
the ginon propagator at zero momentum is finite or not cannot be made 
neither from our data at zero and non-zero momentum nor from the data 
available in the literature at present. 



5.2 The Kugo-Ojima confinement parameter 

In section Sec. 2.3.1 we have introduced the Kugo-Ojima confinement sce- 
nario. According to this scenario, colored asymptotic states, if any, are con- 
fined from the physical state space of QCD in covariant gauges by the quar- 
tet mechanism, if for the function u{q^)5'^^ := with defined in 
Eq. (2.17), the zero-momentum limit 

u :— lim u{q^) — —1 (5.3) 

is realized. As pointed out by KuGO [Kug95], in Landau gauge this limit is 
connected to an infrared diverging ghost dressing function J. This can be 
easily seen from Eq. (2.18) which yields 

-^^^l + u{q') + qMQ') '-^ l + u{Q). 

We have made an attempt to confirm the realization of the limit in 
Eq. (5.3) — here for lattice QCD in Landau gauge — not only by giving 
numerical evidence for a diverging ghost dressing function (see Sec. 4.3.2), 
but also by estimating the function u{q'^) for different momenta q^ from our 
lattice simulations. 



5.2.1 Expected infrared behavior 

Before discussing numerical results it is interesting to figure out first what 
can be expected for the momentum dependence of u{q'^) in the infrared. In 
fact, due to Eq. (2.18) we know that in the limit of vanishing momenta the 
function u{q^) approaches the asymptote 

If we assume that, for example, the power law J(g^,/Li^) oc {q^ / jJ^'^Y'^ de- 
scribes the ghost dressing function at very low momenta, then the asymptote 
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Figure 5.2: The asymptote —u{q^,^^) at low momentum (see Eq. (5.4)) is illustrated 
in logarithmic and linear momentum scale. As ansatz for the ghost dressing function we 
used J{q^, fP) = (g^/Ai^)"" with /i = 4 GeV. The k values are 0.2, 0.4 and 0.6. 



/i^), and so the function /i^), should expose an infrared behavior as 
illustrated in Fig. 5.2 for different infrared exponents k. Even if such a power 
law will turn out not to be appropriate for J{q^, ^J?) — the results presented 
in Sec. 4.3.2 at least give rise to some doubt — it is commonly believed (and 
in agreement with our results) that the ghost dressing function diverges in 
the infrared; and this of course independent of the renormalization point /i 
chosen. Therefore, based on Eq. (2.18) the function u{q^,iJ?) is expected to 
reach minus one and to join u{q'^,fj,'^) at vanishing momentum, irrespective 
of the chosen /i^. 

To confirm this point with numerical data, in Fig. 5.3 we show tt(g^,/i^) 
as obtained from our data for the ghost dressing functions renormalized at 
different momenta /i. In this figure we clearly see that the different curves 
referring to different n approach each other slowly with decreasing momen- 
tum. With respect to Fig. 5.2 we expect that all these curves run slowly 
towards minus one in the zero momentum limit by definition (see Eq. (5.4)) 
if the ghost dressing function diverges in the infrared no matter how. 

5.2.2 Explicit lattice data for the function u(q^) 

Apart from the asymptotes u relying solely on data for the ghost dressing 
functions we have also calculated on the lattice the function 

explicitly in terms of MC averages for u'^^i^q^) at different q^ (see Eq. (3.38)). 
These estimates were obtained on 32^ lattice from our (quenched) gauge-fixed 
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Figure 5.3: The asymptote —u{q^,^'^) as defined in Eq. (5.4) is shown as a function 
of momentum . For the ghost dressing function we used our data at /3 = 5.8 and 6.0 
renormaUzed either at = 3, 4 or 7 GeV. The lattice size is 32^. Lines are drawn to guide 
the eye. 



configurations thermalized at /3 = 5.8 and 6.0 (see runs labeled as S-3 and 
S-6 in Table 4.1). 

Renormalization 

As for other observables the lattice data of uii^q^) = UL{q^,a^) have to be 
renormalized yet, i.e. assuming multiplicative renormalizability we have to 
define a factor Z„ that relates the bare estimate to a renormalized one 

For this we have made use again of Eq. (2.18). It relates the renormalized 
ghost dressing function J(g^,/i^) to the renormalized function ti(g^,/i^) at 
finite . 

J J = 1 + ^(^'' /^') + ^Mq') (5.5) 

As in Sec. 4.3.2, we have renormalized the ghost dressing function at the 
renormalization point /i = 4 GeV and so with Eq. (5.5) we can renormalize 
u(g^,/i^) at the same point. However, the function f(g^) (see [Kug95]) is not 
available to us, causing some difficulties for the determination of Zu{fi'^,a^). 
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A 


B 


C xVndf 


1.146(7) 


0.22(3) 


-0.21(5) 


0.06(3) 0.59 



Table 5.1: The parameters obtained by MINUIT after minimization of the function 
defined in Eq. (5.6). For this we used our data of and J(g^,/Lt^ = (4 GeV)^) for 

<2 GeV^ at /3 = 5.8 and 6.0. The value Z„ has been used to renormalize u&i fi = A GeV 
in Fig. 5.4. 

Nevertheless, we know that for small momenta the term q^v{q^) be- 
comes less dominant in Eq. (5.5). With the ansatz of a Taylor expansion 

v{q^)^A + Bq^ + Cq^ + ... 

we have made a fit using ^ 

and our lattice data for uiiq^) and for J{q^,ii^) renormalized at n — 4 GeV. 
By applying the momentum cut qf < 2 GeV^ we have gained in this way Zu 
and the other parameters A, B and C. The parameters are given in Table 5.1 
together with a value for x^/ndf. 

Of course, the number of terms in the Taylor expansion depends on the re- 
gion of momenta considered. We have found that the given order is necessary, 
but also sufficient in our case. Note that due to the number of parameters 
and the amount of data available for q^ < 2 GeV^ in the minimization pro- 
cess we have not distinguished between data at /? = 5.8 and 6.0. In general, 
the factor Zu{n'^, a^) differs for the two P values, but due to multiplicative 
renormalization there is an additional factor that relates both. From data 
inspection we found that we can approximate this to be one in our case. 

Discussion of numerical results 

In Fig. 5.4 we show data for the renormalized function u^q^^iJ?) aX {3 — 5.8 
and 6.0 as a function of momentum q^. Since the ghost dressing function 
was renormalized at ;U = 4 GeV so was the fimction u{q^, ji^) mediated by 
the minimization process mentioned above. In this figure we also show data 
for the asymptote u{q'^,n'^) renormalized at the same /i. Note that only Z^ 
has been of relevance for providing this figure. 

^For minimization we used the new C-I--I- implementation of the library MINUIT [JW04] 
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Figure 5.4: Data for the function u(<7^,/i^) sX P — 5.8 and 6.0 are shown using full and 
open squares. Additionally, data of the asymptote are shown at the same /3 

values (circles). All data refer to the same quenched configurations on a 32^ lattice and 
are renormalized at /i = 4 GeV as described in the text. Lines are drawn to guide the eye. 



We clearly see that u{q'^,fi'^) and u{q'^,fi'^) approach each other with de- 
creasing momenta. Therefore, as the ghost dressing function is expected to 
diverge in the infrared, we expect to approach minus one as does 

the asymptote u{q^,iJ?) by definition. In Fig. 5.4 the growth of —u{q'^,^'^) 
becomes slower as the momentum decreases, but this behavior we anticipate 
from our discussion above (see Fig. 5.2) concerning the infrared behavior 
expected for u. 

Although we cannot give a reasonable extrapolation of towards 
the zero momentum limit, our results for both and con- 

firm that in the infrared the realization of the limit given in Eq. (5.3) can 
be studied solely in terms of the ghost propagator. As this propagator we 
have found to diverge stronger than 1/g^, the function M(g^,/i^) will reach 
minus one at zero momentum. This is because in the infrared the difference 
|M(g^,/i^) — £t(g^,/i^)| seems to vanish as we have confirmed in this study for 
the first time (see Fig. 5.4). 

Here a remark is in order. The studies of FuRUl and Nakajima [NFOOb; 
NFOOa; FN04a; FN04b]) yield as a limit for M(g^,/i^) values ranging between 
—0.7 and —0.83. We cannot confirm these results in our study, because 
with respect to our data such limits would correspond to linear or quadratic 
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extrapolations towards the zero momentum limit in a figure using a linear 
momentum scale. However, linear or quadratic extrapolations would not be 
valid in this context. 



5.3 The gluon propagator explicitly violates 
reflection positivity 

In this section we will show that our data for the gluon propagator presented 
in Sec. 4.3 show evidence for a violation of reflection positivity. To substanti- 
ate this, a subset of data obtained on our larger lattices has been selected for 
both the quenched and the unquenched case and then the temporal correla- 
tor C{t,p^ — 0) of the gluon propagator has been calculated. The continuum 
expression of this correlator is defined in Eq. (2.23); and the lattice equivalent 
is given by 

C{t) ^ C{ty = 0) := ^|^'D(fc4,0)exp|^| (5.7) 

fc4=0 ^ ^ 

where Lt denotes the number of lattice points in = 4 (time) direction and 
D refers to the gluon propagator in momentum space. 

Starting with the quenched gluon propagator aX j3 — 6.0, the results of the 
correlator are shown in Fig. 5.5 for the lattice sizes 32"^ and 48"^. It is obvious 
from this figure that the gluon propagator violates reflection positivity in a 
finite range of t. The same holds for the gluon propagator in the unquenched 
case as can be seen in Fig. 5.6. There we made use of our measurements 
of the gluon propagator on a 24^ x 48 lattice at different values of (5 and k. 
These values are given in Table 4.2 and the labels D-2, D-3 and D-4 refer to 
these. 

Following Refs. [AvSOl; A+97; CMT05], the statement of reflection-posi- 
tivity violation can be made even more clear by considering the quantity 

G(t,a):=^^^m±2^^}_S:lMA. (5.8) 
This is a discretized expression of the continuum quantity G{t)[C{i)Y where 

Here we adopted the notation [A+97; CMT05] 

(•) / {■)e--'p{u')du . 
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Figure 5.5: The upper panel shows the 
real space propagator C{t) of the gluon 
fields in the quenched approach for differ- 
ent lattice sizes at l3 — 6.0. In the lower 
panel the same data are shown, however, 
as logC(i) for C{t) > 0. 



Figure 5.6: The same as in Fig. 5.5, how- 
ever, for the unquenched gluon propagator 
on a 24^ x 48 lattice. The labels refer to 
the corresponding set of parameters in Ta- 
ble 4.2. 



If the spectral density p{uj^) is positive definite, so is G(t). But if instead 
G(t) is found to be negative then there cannot be a positive definite spectral 
density, which entails violation of refiection positivity. This can easily be seen 
in Fig. 5.7 where values of —G(t,a = 1) are shown for the data considered 
above. There G{t, 1) is clearly negative for all t and therefore refiection 
positivity must be violated. 

Beside of G{t), one usually also discusses in this context the rise of an 
effective gluon mass. On the lattice this can be defined as 

meff(t):=-log|^^|^| (5.9) 

where C{t) again refers to the real space propagator of the gluon fields 
(Eq. (5.7)). Already in the first numerical study of the gluon propagator 
[M087] this effective mass was observed to rise with increasing distance t 
(see also [BPS93; MMST93]). In fact, the definition of meff in Eq. (5.9) 
yields that [CMT05] 

If refiection positivity were satisfied, then G{t, a) would be positive semi- 
definite and therefore meff(t) > meff(t + a). However, for the gluon propaga- 
tor in Landau gauge the opposite was found. This has suggested from the 
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t [l/GeV] i [1/GeV] 

Figure 5.7: The quantity —G{t, a) as defined in Eq. (5.8) is shown as a function of t at 
a—1. Here we used the same quenched (left) and unquenched (right) data as in Fig. 5.5 
and 5.6, respectively. 
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Figure 5.8: The effective gluon mass rries of the quenched (left) and unquenched (right) 
gluon propagator. 



beginning [M087] that the SU{3) lattice gluon propagator violates reflection 
positivity. With our data in Fig. 5.5 and 5.6 we show this explicitly for the 
quenched and unquenched SU{3) gauge theory. 

For the sake of completeness we have calculated the effective gluon mass 
also for our data. The results are plotted in Fig. 5.8 where rries is shown 
as a function of t for the quenched and unquenched case, respectively. The 
effective gluon mass is clearly rising within error bars for moderate values 
oft. 

Summarizing, our study has found a clear signal for the violation of reflec- 
tion positivity for both the quenched and the unquenched gluon propagator. 
If we consider this statement in conjunction with the results presented in the 
former section which support the Kugo-Ojima conflnement criterion to be 
satisfled for lattice QCD in Landau gauge, then this supports the assump- 
tion that the transverse gluon states are confined by the quartet mechanism. 
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Chapter 6 



Spectral properties of the 
Faddeev-Popov operator 



Spectral properties of the Landau gauge FP operator are important for under- 
standing many aspects of QCD in Landau gauge. In this chapter we report on 
a study of some of these properties, restricting ourselves to the quenched approxi- 
mations of lattice QCD. We shall start with a discussion of the low-lying eigenvalue 
distribution where the impact of the Gribov ambiguity is shown in particular. Con- 
cerning the infrared behavior of the ghost propagator discussed above here we will 
analyze the contribution of the different eigenvalues and eigenmodes of the FP op- 
erator to the ghost propagator at low momentum. Localization properties of the 
eigenmodes analyzed in addition. 

6.1 Specification of lattice samples 

Our investigation of the infrared behavior of ghost and giuon propagators (see 
Sec. 4.3) has revealed that unquenching effects are negligible within errors for 
the ghost propagator. We expect the same to hold for the FP operator itself. 
Therefore, we have analyzed its spectral properties solely in the quenched 
approximation of QCD. 

For this analysis we have used a subset of our pure SU{3) gauge config- 
urations thermalized with the standard Wilson action at /3 = 5.8 and 6.2 
using the lattice sizes 12^, 16^ and 24^. To study the influence of the Gribov 
ambiguity we followed again our f c-bc strategy introduced in Sec. 4.1.3 and 
have generated two ensembles of first (f c) and best (be) gauge-fixed config- 
urations. On those ensembles the low-lying eigenvalues A of the FP operator 
and the corresponding eigenmodes have been separately extracted. For this 
we used the parallelized version of the ARPACK package [LMSY], PARPACK. 

To be specific, the 200 lowest (non-trivial) eigenvalues and their corre- 
sponding eigenfunctions have been calculated at (3 = 6.2 using the lattice 
sizes 12'^ and 16'' (see Table 6.1). Due to restricted amount of computing 
time only 50 eigenvalues and eigenmodes have been extracted on the 24^ 
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No. 


/? 


lattice 


# conf 


# copies 


# eigenvalues 


F-1 


5.8 


24^ 


25 


40 


90 


F-2 


6.2 


12^ 


150 


20 


200 


F-3 


6.2 


16^ 


100 


30 


200 


F-4 


6.2 


24^ 


35 


30 


50 



Table 6.1: Statistics of data used in our analysis. The last column lists the number of 
eigenvalues extracted separately on f c and be copies of U . At /3 = 6.2 the corresponding 
eigenmodes were calculated, too. Labels given in the first row refer to the corresponding 
entries in Table 4.1. 

lattice at the same (5. In addition, 90 eigenvalues have been calculated on a 
24^ lattice at /3 = 5.8 providing us with an even larger physical volume. This 
allows us to check whether low-lying eigenvalues are shifted towards A ^ 
as the physical volume is increased. The eight (trivial) zero eigenvalues with 
the corresponding constant zero modes have always been discarded. 

6.2 The low- lying eigenvalue spectrum 

We start the discussion with the two lowest (nontrivial) eigenvalues of the FP 
operator and then give an estimate for the density of low-lying eigenvalues. 

6.2.1 The lowest and second lowest eigenvalues 

The distributions of the two lowest-lying eigenvalues, Ai and A2, of the FP 
operator are shown for different volumes in Fig. 6.1. There h{X,X + AA) 
represents the average number (per configuration) of eigenvalues found in 
the interval [A, A + AA]. To disentangle the distributions for the two different 
sets of gauge copies, open (full) bars refer to the distribution on f c (be) 
gauge copies. 

It is obvious from this figure that both eigenvalues, Ai and A2, are shifted 
to lower values as the physical volume is increased. In conjunction the spread 
of A values is shrinking. This would be even more obvious, if we had shown 
both distributions as functions of A in physical units. 

It is also visible that the two low-lying eigenvalues A*"^ {i = 1, 2) on 
f c gauge copies tend to be lower than those on be copies. However, this 
holds only on average as can be seen from Fig. 6.2. There the differences 
A^'' — A*'' of the lowest eigenvalues on f c and be gauge copies are shown for 
different lattice sizes at (3 = 6.2 and 5.8. It is evident that there are only few 
cases where A^'^ < A^"^, even though F^'^ > F^"^ always holds for the gauge 
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Figure 6.1: The frequency h{X) per configuration of the lowest (left panels) and second 
lowest (right panels) eigenvalue A of the Faddeev-Popov operator is shown. Filled boxes 
represent the distribution obtained on best (be) gauge copies, while open ones represent 
those on first (f c) copies. 



functional. 

In addition we have checked how the average values (A) of the respective 
eigenvalue distributions tends towards zero as the linear extension aL of the 
physical volume is growing. For this the lattice spacing a has been specified 
in physical units. As in Chapt. 4 we followed Ref. [NS02] to fix a. For 
P = 5.8 and 6.2 we used a~^=1.446 GeV and 2.914 GeV, respectively, using 
the Sommer scale tq = 0.5 fm. 

If the low-lying eigenvalues are supplemented with physical units it turns 
out that the average values of their distributions tend towards zero stronger 
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20 40 60 80 
conf 

Figure 6.2: The differences A^" — Af" of the lowest FP eigenvalues calculated on be and 
f c representatives for each gauge configuration are shown. From left to right the lattice 
sizes are 12*, 16* and 24* at /3 = 6.2 and 24* at /3 = 5.8. 



o?f{aL) 


C 


e 


xVndf 


(Ai) 


0.120(3) 


0.16(4) 


0.7 


(A2) 


0.165(4) 


0.24(5) 


1.8 


(A5) 


0.290(1) 


0.45(4) 


3.5 


(A2) - (Ai) 


0.045(2) 


0.47(9) 


0.4 


(A3) - (A2) 


0.051(1) 


0.88(8) 


0.2 


(A4) - (A3) 


0.033(1) 


0.62(33) 


2.0 


(A5) - (A4) 


0.037(1) 


0.89(1) 


0.003 



Table 6.2: The parameters C and e from fitting either the averages {\i)/a? or the dif- 
ferences of adjacent average values {Xi+i)/a'^ — {Xi)/a'^ of the corresponding eigenvalue 
distributions to the ansatz f{aL) — Ci/(aL)^+''. 



than by volume scaling proportional to l/(aL)^. In fact, using the ansatz 



to fit the data of {Xi) /a? for diff'erent (aL), always a positive e is found. The 
parameter of these fits are given in Table 6.2 and in Fig. 6.3 we show the data 
and the corresponding fitting functions. There one clearly sees, the low-lying 
eigenvalues not only approach zero, but also become closer to each other with 
increasing aL. The latter issue has been addressed by fitting the differences 
((Aj+i) — {\i))/a? of adjacent average values using the same ansatz Eq. (6.1). 
In Table 6.2 we give the parameter of such fits. 
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l/(aL) [fill-'] 



Figure 6.3: The average values {\i)/a? 
(scaled to physical units) of the eigenvalues 
\i (i = 1, 2, 5, 20) are shown vs. the inverse 
of the linear lattice extension aL. Only 
eigenvalues on be copies are shown. The 
lines represent fits to the data using the 
ansatz a~^(Ai) = Ci/(aL)^+'^' . 




0.02 0.04 0.06 
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Figure 6.4: The eigenvalue density p for 
be copies as a function of A estimated on 
12* and 16'' lattices at /? = 6.2 and on a 24^^ 
lattice for (3 = 6.2 and 5.8. Bin sizes have 
been chosen as small as possible for each 
lattice size. Symbols mark the middle of 
each bin and lines are to guide the eye. 



6.2.2 An estimate for the density of low- lying 
eigenvalues 

The eigenvalue density p(A) is of particular interest. At small A this quantity 
has been estimated here by 

i.e. the average number h of eigenvalues per gauge-fixed configuration within 
the interval [A, A + AA] divided by the bin size AA. For normalization the 
denominator = 81^ has been chosen, since the FP matrix is a x sparse 
symmetric matrix with A^ linearly independent eigenstates. The trivial zero 
modes would be described by a term 85(A) in p(A) (not shown). 

The estimates for the density p are shown in Fig. 6.4 for the different 
volumes used. The bin sizes have been reasonably adjusted for each volume 
separately. In Fig. 6.4 one clearly sees the eigenvalue density close to A = 
becomes steeper as a function of A as the physical volume becomes larger. It 
is remarkable that the increase going from (3 = 6.2 to /3 = 5.8 on a 24^^ lattice 
is larger than going from 12*^ to 24"^ at /3 = 6.2 fixed, although in both cases 
the physical volume is increased by a factor of about 16. 
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6.3 Eigenmode expansion of the ghost 
propagator 

Along with the calculation of the low-lying eigenvalues, the corresponding 
eigenvectors 4>{x) have been determined as well. These are of particular 
importance for the infrared behavior of the ghost propagator as it becomes 
clear from Eq. (3.37). In fact, if all eigenvalues Aj of the FP operator and 
the corresponding eigenvectors ^i{k) in momentum space were available the 
ghost propagator could be constructed out of them according to Eq. (3.37) 
and (3.36). Unfortunately, their determination for each configuration is nu- 
merically too demanding. 

However, restricting the sum in Eq. (3.37) to the n lowest eigenvalues 
and eigenvectors (n -C iV = 8V — 8), we can figure out to what extent these 
modes, i.e. the corresponding estimator Eq. (3.36) and (3.37) 



saturates the full ghost propagator G{(f). The latter is obtained, of course, 
independently for a set of momenta by inverting the FP matrix on a set of 
plane waves. See Sec. 4.3.2 for the corresponding data of G{q^). 

The degree of saturation is shown in Fig. 6.5 for the lowest q\ and the 
second lowest momentum available on different lattice sizes for /3 = 6.2. 
There the values of Gn{q^) have been presented relative to the values for the 
full propagator G(g^) in order to assess the saturation for different volumes. 
Since ^i{k) has been obtained by a fast Fourier transformation of the eigen- 
vector <pi{x), all lattice momenta k are available. Thus Gn{q^) refers to the 
average over all k giving raise to the same momentum q^. The full propagator 
values G{q^) at q\{k) and gf(/c), however, refer to the averages over lattice 
momenta k — ([1, 0], 0, 0) and to A; = (1, 1, 0, 0), respectively. 

Let us consider first the lowest momentum q\. We observe from Fig. 6.5 
that the approach to convergence differs, albeit slightly, for the three different 
lattice sizes. The relative deficit for n < 50 rises with the lattice volume. For 
n > 100 the rate on a 16^ lattice is even a bit larger than that on a 12^ lattice. 
We could not afford to get data for n > 50 on the 24^ lattice. However, for 
the 12^ and 16^ lattices the rates of convergence are about the same. For 
example, taking only 20 eigenmodes one is definitely far from saturation (by 
about 50%) whereas 150 to 200 eigenmodes are sufficient to reproduce the 



Gniq^k)) = {G{k\n))MC 



where 
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n 



Figure 6.5: The ratio of the truncated ghost propagator G„(g^) (in terms of the n lowest 
FP eigenmodes and eigenvalues) to the full estimate G{q^) (taken from [SIMPSOSd]) shown 
as a function of n for the lowest (qf) and second lowest (g|) momentum. The inverse 
coupling is /3 = 6.2 and the lattice size ranges from 12** to 24**. All data refer to be copies. 

ghost propagator within a few percent. In other words, the ghost propagator 
at lowest momentum on a 12^ (16^) lattice is formed by the lowest 0.12% 
(0.03%) of the eigenvalues and eigenfunctions of the FP operator. 

For the second lowest momentum g| the contribution of even 200 eigen- 
modes is far from being sufficient to approximate the propagator. 

6.4 The problem of exceptional configurations 

We turn now to a peculiarity of the ghost propagator at larger /3 of which 
has been reported already in [SIMPS05d; SIMP06]. It was also seen in an 
earlier SU (2) study [BIMMP04]. While inspecting our data we found, though 
rarely, that there are exceptionally large values in the Monte Carlo (MC) time 
histories of the ghost propagator at lowest momentum. 

In Fig. 6.6 such time histories are shown of the ghost propagator on f c and 
be gauge copies for the two smallest momentum realizations k = (1, 0, 0, 0) 
and k = (0,1,0,0). From left to right the panels are ordered in ascending 
order with the (physical) lattice sizes 12"^, 16^ and 24"^ at /3 = 6.2 and 24"^ at 
p = 5.8. 

As can be seen from this figure in the majority extreme spikes are re- 
duced (or even not seen) when the ghost propagator could be afforded to be 
measured on a better gauge copy (be) for a particular configuration. Further- 
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I2J ,1J = 6.2 m' ,d = 6.2 24' /.) = (i.2 24-" ;3 = 5.8 
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Figure 6.6: The MC time histories of the ghost propagator calculated on first (f c) and 
best (be) gauge copies on a 12*, 16^ and 24* lattice at /? = 6.2 and 24* lattice at /3 = 5.8. 
From left to right the corresponding runs listed in Table 6.1 and Table 4.1 are F-2, F-3, 
F-4/S-9 and F-l/S-2. The upper and lower panels show data for the lowest momentum 
realization k = (1, 0, 0, 0) and k — (0, 1, 0, 0), respectively. 

more, it is obvious that the exceptionality of a given gauge copy is exhibited 
not simultaneously for all different realizations of the lowest momentum. 
Consequently, to reduce the impact of such large values on the evaluation of 
a statistical average for the ghost propagator one should average, if possible, 
over all momentum realizations giving rise to the same momentum g^. 

We have tried to find a correlation of such outliers in the history of the 
ghost propagator with other quantities measured in our simulations. For 
example we have checked whether there is a correlation between the values 
of the ghost propagator G{k) as they appear in Fig. 6.6 and the lowest 
eigenvalue Ai of the FP operator. 

In Fig. 6.7 we show such correlation in a scatter plot for different lattice 
sizes at P = 5.8 and 6.2. There each entry corresponds to a pair [Xi,G{k)] 
measured on a given gauge copy of our sets of f c and be copies. It is visible 
in this figure, gauge copies giving rise to an extremely large MC value for the 
ghost propagator are those with very low values for Ai. This holds for the 
12"^ and 16*^ lattice. However, a very low eigenvalue is not sufficient to obtain 
large MC values for G{k) as can be seen in the same figure. It is possible, of 
course, that such gauge copies with extremely small eigenvalues would turn 
out to be exceptional for another realization of lowest momentum g^(/c) than 
those two we have used. This might explain why some configurations with 
extremely small lowest eigenvalues were not found to be exceptional with 
respect to the ghost propagator at k = (1, 0, 0, 0) and k = (0, 1, 0, 0). 

In the light of Eq. (3.37) it is not adequate to concentrate just on the 
lowest eigenvalues. Instead, one can monitor the contribution of a certain 
number of eigenvalues Aj and eigenmodes ^i{k) to the ghost propagator at 



6.4 The problem of exceptional configurations 



119 




50 - 



40 



30 



(3 = 6.2 24" _ 48 

□ 

□ n 

44 

□ — 



1%' *■ 

I," I". m, 

B 



0.004 0.008 
Ai 



36 



(3 = 5.8 24^ 



0.004 

Ai 



0.008 



Figure 6.7: Scatter plots of MC time history values of the ghost propagator G{k) at 
k = ([1,0], 0,0) vs. the lowest FP eigenvalue Ai are shown. The upper panels show data 
at /3 = 6.2 on a 12^ (left) and 16'* (right) lattice, the lower ones on a 24* lattice at f3 — 6.2 
(left) and P = 5.8 (right). Open symbols refer to f c gauge copies and filled symbols to be 
copies. 
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Figure 6.8: Scatter plot of the relative contribution of the truncated sums G{k\n) over 
eigenmodes (see Eq. (3.37)) to the full ghost propagator values G{k) versus G{k) for lattice 
momenta k = (1, 0, 0, 0) and k = (0, 1, 0, 0). From left to right the lattice sizes are 12'', 16* 
and 24* all for f3 = 6.2. Data for f c and be gauge copies have been plotted separately. 
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some momentum in question. Therefore, we have compared the truncated 
sums G{k\n) according to Eq. (3.37) with the MC history values of the full 
ghost propagator G. In fact, we show in the scatter plots in Fig. 6.8 the 
ratios G{k\n)/G{k) versus G{k) for n = 10 and for various lattice sizes. 

Obviously there is a strong correlation between the chosen group of low- 
lying modes and the MC time history values of the full ghost propagator. 
Indeed, if we consider values G{k) > 15 to be exceptional in the left-most 
panel (12^ lattice) we find that the contribution of the 10 lowest modes 
amounts to more than 75% of the actual value of the ghost propagator. On 
the opposite, for low G{k) values the main contributions come necessarily 
from the higher eigenmodes, while the 10 lowest modes contribute a minor 
part only. A similar but less dominant contribution of the 10 lowest modes 
is found for the time histories produced on larger lattices (16^, 24^). 

6.5 Localization properties of low- lying 
eigenvectors 

In recent years a good deal of attention has been directed to the localization 
properties of various operators (Dirac operator, covariant Laplacian) in the 
hope to understand more about confinement beyond heavy quark probes. 
The first quantity of interest is the inverse participation ratio (IPR). Given 
an eigenvector (f){x) the IPR is defined as 



Although a direct physical meaning for this quantity is lacking yet, it is 
a measure for the localization of an eigenvector. It enables us to distin- 
guish between eigenmodes with approximately uniformly distributed modu- 
lus squared |0(a;) p (IPR pa 1 ... 2) and more specific ones with a small number 
of sites X having large intensity |</)(x)p (IPR ~ O(IOO)) where they might be 
pinned down by special local gauge field excitations. Note, the eight (trivial) 
zero modes (A = 0) of the FP operator are constant and have IPR = 1. 

In Fig. 6.9 the relative distribution h of IPR values per gauge-fixed config- 
uration are shown, separately for certain groups of eigenstates. Again open 
(full) histogram bars refer to the distribution on f c (be) gauge copies. From 
this figure we learn that the majority of eigenvectors of the Faddeev-Popov 
operator is not localized independent of the choice of gauge copies. In any 
case, the rare cases of large IPR values have been found among the 10 lowest 
non-zero eigenmodes. This becomes more likely as the physical volume is in- 
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Figure 6.9: The relative distribution h of IPR values of the 10 (left), the 10 to 20 (middle) 
and the 20 to 50 (right) lowest eigenmodes are shown. Note there is a logarithmic scale 
for h{lPR). Each row corresponds to one pair of /3 and lattice size. Filled boxes refer to 
distributions on be gauge copies, while open ones correspond to f c copies. 



creased. So far we have no physical interpretation what causes the stronger 
localization in these rare cases. 
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Conclusions and outlook 



In this study we have focused on the infrared properties of SU (3) gluody- 
namics in Landau gauge using the framework of lattice QCD. We have 
tried to cover as much as possible the different aspects relevant within this 
context and to verify several assumptions made in recent years. In doing so, 
we hope to have provided a consistent analysis of several issues that have an 
effect on the infrared behavior of gluon and ghost propagators, and that, in 
turn, the behavior we have found satisfies necessary criteria for confinement 
which apply to QCD in Landau gauge. 

The method 

For this study we used the Wilson formulation of lattice QCD with and 
without dynamical clover-improved Wilson fermions. The gauge group was 
fixed to S'f/(3). The quenched gauge configurations were generated at the 
three values = 5.8, 6.0 and 6.2 using a variety of different symmetric and 
asymmetric lattice geometries. In the symmetric case we used the lattice 
sizes 12^, 16^, 24^, 32^ and 48^, whereas the asymmetric lattices were chosen 
to be of size 24=^ x 48, 32^ x 64, 16^ x 128 and 24^ x 128. To study unquenching 
effects we have analyzed gauge configurations provided to us by the QCDSF 
collaboration. These configurations were thermalized on a 24^ x 48 lattice in 
the presence of two fiavors of clover-improved Wilson fermions using three 
different pairs of (3 and k. The corresponding lattice spacings are comparable 
to that aX (5 — 6.0 in the quenched case. 

Both the quenched and unquenched configurations were transformed such 
that they satisfy the lattice Landau gauge condition. For gauge-fixing we 
used either over-relaxation or Fourier-accelerated gauge-fixing. A subset of 
our quenched configurations was gauge-fixed even more than once, always 
starting from a different random gauge copy of the initial (unfixed) configu- 
ration. This has allowed us to examine how the Gribov ambiguity affects the 
gluon and ghost propagators and the eigenvalue spectrum of the FP operator. 

The results 

In the following we give a summary of our results presented in previous 
chapters. Then we draw our conclusions and give recommendations for future 
studies. 
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The influence of the Gribov ambiguity 

We have demonstrated that the presence of Gribov copies systematically 
affects the ghost propagator at low momentum, whereas for the gluon prop- 
agator such an effect stays within error bars (see also [SIMPS05d]). To be 
specific: Measuring the ghost propagator and ignoring the Gribov ambiguity, 
the ghost propagator near the momentum = 0.2 GeV^ (1 GeV^) turns out 
to be overestimated by about 5% (2%) compared to an estimate obtained on 
a ensemble of best gauge copies. As best we have considered that gauge-fixed 
copy which gave rise to the largest gauge functional value for a particular 
gauge configuration. Our results corroborate previous findings for the SU{2) 
gluon and ghost propagators [Cuc97; BIMMP04; NF04a], but cast doubt on 
those for the SU{3) gluon propagator {j3 = 5.8, 12"^ lattice) in [SO04]. 

Additionally, our data are in favor of the picture promoted in [Zwa04]. 
According to this, continuum vacuum expectation values of correlation func- 
tions obtained from a functional integration over the fundamental modular 
region A are equal to those over the Gribov region fl. Gribov copies inside 
n should not affect expectation values in the continuum, because functional 
integrals are dominated by the common boundary of A and il. We have 
found some numerical evidence that the influence of Gribov copies on the 
ghost propagator decreases at the same (physical) momentum if the physical 
volume is enlarged. Note that very recently [BBMPM05] similar indications 
have been found for the SU{2) gauge group taking non-periodic Z{2) trans- 
formations into account. In order to eliminate the last doubts, a future 
study should continue and follow our f c-bc strategy (explained in the text) 
on lattice sizes larger than 24^. 

To make the study of the Gribov-copy dependence more complete, we 
have also shown that the Gribov ambiguity is reflected in the low-lying eigen- 
value spectrum of the FP operator. We have found that, on average, the 
low-lying eigenvalues extracted on be gauge copies are larger than those on 
f c copies. Thus better gauge-fixing (in terms of the gauge functional) has the 
tendency to keep gauge-fixed configurations slightly away from the Gribov 
horizon. 

Other systematic eff'ects on the propagators 

Using different lattices sizes at three different (3 values, we have tried to 
analyze the systematic effects on the gluon and ghost propagators of changing 
either the lattice spacing a or the physical volume V. We have found that 
for both, the gluon and ghost dressing functions, finite volume effects are 
clearly visible at volumes smaller than (2.2 fm)*^, which corresponds to a 
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16^ lattice at /3 = 5.8. The effect grows with decreasing momentum or 
decreasing lattice size. At larger volumes, however, the data for g > 1 GeV 
coincide within errors for the different lattice sizes. For q < 1 GeV we have 
found only small finite volume effects for both dressing functions at the lowest 
momentum. This is based on data obtained on the lattice sizes 24'^, 32^ and 
48^ using /3 = 5.8 and 6.0. Concerning discretization effects our study is 
only partial and limited to a region of intermediate momenta. In any case, 
fixing the physical volume to V ^ (2.2 fm)^ we have found that the gluon 
dressing function at constant physical momentum increases with decreasing 
the lattice spacing. A similar effect (beyond error bars) is not observable for 
the ghost dressing function. 

We have combined this investigation of lattice artifacts with an analysis 
of effects caused by asymmetric lattice geometries. We could demonstrate 
that the more asymmetric the lattice size has been chosen, the larger are 
the systematic errors induced by that. In fact, data obtained for the ghost 
(gluon) propagator on asymmetric lattices are less enhanced (suppressed) 
than those obtained on symmetric lattices geometries. The same effects have 
been reported recently for the gluon propagator in three-dimensional pure 
SU(2) gauge theory [CM06]. Therefore, extractions of an infrared exponent 
for the gluon propagator using only on-axis momenta on asymmetric lat- 
tices [SO05b; SOOoa] (without adapting the lattice spacing in the different 
directions to compensate this) should be taken with caution. 

Infrared behavior of gluon and ghost propagators 

Studying the momentum dependence of the gluon and ghost dressing func- 
tions, we have applied cuts on our data, namely the cylinder and the cone cut 
[LSWP99]. They have much reduced the lattice artifacts mentioned above. 
In fact, our data for the renormalized dressing functions surviving these cuts 
lie on smooth curves if considered as functions of the momentum. 

We have compared these renormalized data with recent solutions of trun- 
cated system of DSEs for the gluon and ghost propagators and have tried to 
extract the infrared exponents, and kq, for the gluon and ghost dressing 
functions, respectively. From fits with the corresponding power laws to our 
data we cannot confirm the relation kq = 2ko as expected from DSE studies 
in the continuum [vSAH97; vSHAQS]. Also, both exponents are found to be 
significantly lower than expected in [LvS02; Zwa02]. Such a conclusion was 
also drawn in other lattice studies [B"'"05b; B+OGa; FN04b]. 

Moreover, our fits suggest that power-behaved gluon and ghost propaga- 
tors are not the best description of the momentum dependence of our data, 
at least in the region of lower momenta considered here. The data for the 
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ghost propagator seem to depend logarithmically on the momenta in this 
region. This observation has been confirmed very recently in [B"^06a], even 
though in this reference also arguments for an infrared finite ghost dressing 
function have been put forward. 

Studying unquenching effects on the gluon and ghost propagators we have 
found that these effects are small for the ghost propagator, but are clearly 
visible for the gluon propagator at intermediate momenta. Concerning the 
infrared limit of both propagators, the influence of two fermion flavors seems 
to become less towards lower momenta. This all agrees with the findings in 
DSE studies (e.g. [FA03; FAC+06]). 

From our data, we cannot judge without doubt on the existence of an 
infrared vanishing gluon propagator. 

The running coupling constant and the ghost-gluon vertex 

In connection with the infrared behavior of the gluon and ghost dressing 
functions, we have determined the running coupling constant based on the 
ghost-gluon vertex. This coupling constant is found to match the RG- 
invariant two-loop expression considering data at large momenta. However, 
for g2 < 0.4 GeV^ the running coupling constant decreases with decreas- 
ing momentum. The same is observed considering data for the unquenched 
case. Also an influence of Gribov copies cannot be made responsible for this. 
Therefore, we cannot confirm an infrared fixed point for this coupling con- 
stant as it has been proposed in DSE studies [vSAH97; vSHA98] (see also 
[AFLE05]). This reflects once more the different infrared exponents we have 
found for the gluon and ghost propagators. 

In any case, our data are in qualitative agreement with recent studies of 
DSEs on a torus [R\R02; FA02; FGA06: FP06]. There a similar behavior 
for the gluon and ghost propagators and for the running coupling constant 
at low momentum has been presented. Also other lattice studies agree with 
our results for this coupling constant [FN04a; FN04b; B''"05b]. 

In order to verify that the assumption of a bare ghost-gluon vertex is valid 
beyond perturbation theory, we have calculated the corresponding renormal- 
ization constant Z\ using a MOM scheme with zero gluon momentum. Our 
data show that in this scheme Z\ is approximately equal to one within error 
bars for all momenta considered here. Only a slight deviation is visible in the 
interval 0.3 GeV^ < < 2 GeV^. Unquenching effects are not resolvable. 
We thus agree with a recent DSE study [SMWA05] where a semiperturbative 
calculation of Z\ has been presented using a MOM scheme with asymmetric 
(the same as used by us) and symmetric subtraction points. Furthermore, 
our results are in full agreement with those presented in [CMM04] for the 
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case oi SU{2). It is worthwhile to continue the lattice calculation of Zi using 
other renormalization schemes, for example with a symmetric subtraction 
point. 

Together with our data for the running coupling constant we can thus 
confirm that the product in Eq. (2.10) is renormalization-group invariant in 
the particular MOM scheme considered here and defines a nonperturbative 
running coupling constant which, as shown here, decreases monotonously 
with decreasing momentum < 0.4 GeV^. 

Results on the confinement criteria 

In addition we have presented results which support the Kugo-Ojima con- 
finement scenario [K079] to be realized for lattice QCD in Landau gauge. To 
confirm this, one has to show that the function = u{q'^)6°'^ defined via 

the Green's function in Eq. (2.17) [Kug95; AvSOl] has the zero-momentum 
limit: u{0) = —1. Alternatively, if this limit is realized then the (renormal- 
ized) ghost dressing function J{q^) must diverge in the same limit. This is 
because at zero momentum it holds that 1 -|- u{0) = 1/J(0) [Kug95]. 

In this thesis we have shown that the ghost dressing function seems to 
diverge at zero momentum which is in favor of the Kugo-Ojima confinement 
scenario and also satisfies the Gribov-Zwanziger horizon condition [Gri78; 
Zwa93; Zwa02]. Moreover, we have measured the function u{q^) at different 
momenta g^. To our knowledge a direct calculation of M(g^) has never been 
done before. Conclusions in other studies [FN04a; FN04b; WAOl; AvSWOl] 
were drawn on results obtained for the ghost dressing functions only. For 
the renormalization of u{q'^) we have developed a new ansatz using a min- 
imization process (see Sec. 5.2.2 for details). Our renormalized data for 
both the ghost dressing function and the function M(g^) are consistent with 
M(g^) — J~^{q^) approaching minus one in the limit of vanishing momenta, 
even though this convergence is very slow. Therefore, as the ghost dressing 
function diverges we expect u{q^) to reach minus one. Hence the Kugo-Ojima 
confinement criterion is realized. 

Based on our data obtained on larger lattice sizes we are also in the for- 
tunate position to present numerical evidence that not only the quenched 
but also the unquenched iS'f/(3) giuon propagator in Landau gauge violates 
reflection positivity explicitly. Thus, on one hand, we agree with other lattice 
studies for quenched SU{2,) [M087; BPS93; MMST93; FN04b] and for three- 
dimensional SU{2) [CMT05] gauge theory. On the other hand, we confirm 
presently available solutions to the corresponding DSE of the gluon propa- 
gator [ADFM04], even though our data suggest a different infrared exponent 
for the gluon propagator. 
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Since reflection positivity is violated by the gluon propagator and the 
Kugo-Ojima conflnement criterion seems to be satisfled, we agree with the 
conjecture that transverse gluon states are conflned by the quartet mecha- 
nism [K079]. 

Spectral properties of the FP operator 

We have also investigated the spectral properties of the FP operator and 
their relation to the ghost propagator in SU{3) Landau gauge [SIMP06]. As 
expected from Ref. [Zwa04] we have found that the low-lying eigenvalues are 
shifted towards A = and the eigenvalue density p(A) becomes a steeper 
rising function as the volume is increased. 

On average, the corresponding FP eigenmodes are not localized, however, 
a few large IPR values have been seen among the lowest eigenmodes. 

We could demonstrate that the ghost propagator at low momentum is 
dominated by the low-lying eigenvalues and eigenmodes of the FP operator. 
For example, the value of the ghost propagator at lowest momentum (on a 
12^ lattice at /3 = 6.2) can be estimated using a number of 200 low-lying 
eigenvalues and eigenmodes of the FP operator. In other words, a fraction of 
about 0.12% of the whole spectrum is sufficient to reconstruct the asymptotic 
result. For larger volumes the number of necessary eigenmodes seems to be 
somewhat larger. Also for the next higher momentum, saturation needs a 
much bigger part of the low-lying spectrum. 

Concluding remarks 

We hope to have presented a careful numerical study of different aspects of 
SU{3) Landau gauge gluodynamics, a subject which has been the focus of 
much attention in recent years. We have concentrated on the low momentum 
region and have clarified the infiuence of difi'erent systematic effects on the 
infrared behavior of gluon and ghost propagators in Landau gauge. We have 
shown that the momentum dependence of both propagators is consistent 
with different criteria for confinement. Considering also our results for the 
running coupling constant and for the ghost-gluon vertex we agree with the 
findings of other lattice studies and of DSE studies on a torus. However, 
there is some disagreement with DSE studies in the continuum. At present 
there seems to be no solution of this puzzle. 

We think, it is interesting to look in a future study also at the momentum 
dependence of other vertex functions and of the quark propagator in Lan- 
dau gauge. In particular, for clover-improved Wilson fermions this has not 
been done before. It is also worthwhile to perform a similar study using an- 
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other gauge condition, for example, the Coulomb gauge. Such investigations 
might provide further valuable information towards a full understanding of 
nonperturbative QCD (see e.g. [Zwa98; CZ02; ZwaOSa; GO03]). 

In the same way, the larger momentum region should be investigated 
further. Note that our study has focused only on the infrared momen- 
tum region. Since computing time was restricted, the amount of data for 
the ghost propagator and for the running coupling constant in the ultra- 
violet momentum region is not as large as in the low momentum region. 
It is worthwhile, however, to continue and to perform additional measure- 
ments at larger momenta considering different loop-expansions of the corre- 
sponding asymptotic form. In particular, an investigation of nonperturbative 
power corrections to the ghost and gluon propagators in Landau gauge (see 
e.g. [B+OOb; B+OOa; B+01; B+06b]) due to non-zero values of QCD conden- 
sates could be an interesting topic in future lattice studies. 
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Appendix A 



Some details on algorithms 

and performance 



e list algorithms and numerical libraries used in this study for the differ- 



V V ent purposes and we give references for further details. Then we report 
on some experiences we gained with gauge-fixing. In particular, we compare the 
performance of over-relaxation and Fourier-accelerated gauge-fixing. After this, 
we analyze gauge functional values obtained by following our f c-bc strategy and 
demonstrate that the final ranking of functional values is already visible at an in- 
termediate iteration step. We show that the inversion of the FP matrix can be 
accelerated by using a pre-conditioned conjugate gradient algorithm. 



A.l A note on the algorithms used 

For the generation of our quenched configurations we employed the hybrid 
over-relaxation algorithm. It has become standard for the simulation of pure 
gauge theory and is a combination of several, say Ngj}, micro- canonical over- 
relaxation steps [BW87; Cre87] and one heatbath [FH84; KPSo] step. In 
both steps a decomposition of SU{3) link variables into SU{2) matrices, as 
proposed by Cabbibo and Marin ARI [CM82], was applied. We have always 
used Nov — 4. Further details on these algorithms can be found, for examples, 
in the PhD theses [Kne99; Geh02]. 

All sets of our unquenched SU{3) gauge configurations were provided to 
us through the QCDSF collaboration which generated them using the hybrid 
Monte Carlo algorithm [DKPR87; GLT+87] with even-odd preconditioning 



Both the quenched and unquenched configurations were transformed such 
that they satisfy the lattice Landau gauge condition (see Eq. (4.1)). For 
gauge-fixing we have employed tw-o popular algorithms: over-relaxation (RLX) 
[MO90a] and Fourier-accelerated gauge-fi/xing (FAG) [D'^88]. For Fourier- 
accelerated gauge-fixing we used the implementation as introduced in [D^88]. 




[G+89]. 
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For future studies, however, we recommend to use the multigrid implementa- 
tion of FAG as proposed in [CM98]. This implementation is better suited for 
a parallel computing environment with distributed memory (see App. A. 2.1). 

Fourier transforms, for example of the gluon fields, were calculated using 
an algorithm that performs a Fast-Fourier transformation (FFT). For all 
FFTs, we have employed the FFTW-library [FJ98]. We can recommend this 
library (see also online: http://www.fftw.org). 

For observables which involve the inverse of the FP operator we applied 
the pre-conditioned conjugate gradient algorithm to solve the corresponding 
linear systems. As pre-conditioning matrix we used the inverse Laplacian 
operator A^^ with diagonal color substructure. This significantly has reduced 
the amount of computing time as it is discussed in more detail in App. A. 3. 

The eigenvalues and eigenmodes of the FP operator were calculated using 
the ARPACK package, PARPACK [LMSY]. For minimization of functions we 
used the new C++ implementation of the library MINUIT [JW04]. We can 
recommend both libraries too. 

A. 2 Experience report on lattice Landau 
gauge fixing 

In this section we report on some experiences we made with gauge-fixing. 
Since efficiency of gauge-fixing algorithms is an important issue in studying 
gauge-dependent quantities, this report might help in future studies. 

A. 2.1 Over-relaxation versus Fourier-accelerated 
gauge- fixing 

For the purpose of this study we have employed two different algorithms 
for fixing thermalized gauge configurations to Landau gauge. These are 
over-relaxation (RLX) [MO90a] and Fourier-accelerated gauge-fixing (FAG) 
[D+88]. Both algorithms are commonly used in the lattice community, but 
usually either of them is chosen. The reason why we used two is simple: 
Initially, we had started with over-relaxation, but later we checked whether 
computing time can be saved using FAG instead. Now this puts us in the 
fortunate position to provide an explicit cross-check on the performance of 
both algorithms. This might be interesting for future studies, in particular 
when writing applications for computing time. 

Even though both algorithms are quite different, they are both iterative 
in nature. After each iteration-cycle the functional Fu[g] (Eq. (3.14)) is 
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24^x48' 



323x64 



Figure A.l: The average number of iterations for both gauge-fixing using over- relaxation 
(rlx) and gauge-fixing using the Fourier-accelerated method (fag) is shown in the upper 
panel. The lower panel shows the corresponding CPU time in seconds (sum over all 
processors). 

increased a bit until a (local) maximum is reached. As stopping criterion 
not the functional itself, but the violation of transversality (see Eq. (3.15)) 
is used. In our implementation the iteration process stopped as soon as 
e < 10-1^ with 



was fulfilled at each lattice site x. The number of necessary iterations de- 
pends on the lattice sizes, but usually it also varies quite strongly for different 
configurations using the same lattice size. To get an impression about these 
variations have a look at Fig. A. 2. There 9 is shown versus the numbers of 
iterations for two different random gauge copies of the same gauge configu- 
ration. 

In any case, the mean iteration number Nuer and the mean iteration 
time Titer (until the final precision is reached) scale with the lattice size 
V = Lg X Lt according to 





Nn,,{V-a,b) 
Titer (V"; c, d) 



(A.2a) 
(A.2b) 



Ls and Lt denote the number of lattice points in spatial and temporal direc- 
tion, respectively. That this scaling holds for both algorithms can be seen in 
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A^iter 






T- 

iter 








a b 


X^/ndf 


c 




d 


X^/ndf 


rlx 


4.5(2) 0.53(1) 


1.6 


3.2(4) ■ 


10"^ 


1.55(1) 


5.6 


fag 


19(2) 0.38(1) 


3.0 


5.0(9) ■ 


10"^ 


1.51(2) 


8.1 



Table A.l: The parameters of the functions iter(F;a, 6) and cpu(F; c, c?) fitted to the 
average number of iterations and average CPU time, respectively. 



Fig. A.l. There in the upper (lower) panel we show A^iter (7]ter) as a function 
of the lattice size. The parameters from fits of the Ansatze in Eq. (A. 2) to 
the data are given in Table A.l. Obviously, FAG performs better on larger 
lattice sizes than RLX. 

The reason why in our implementation the actual amount of computing 
time for FAG is comparable to that of RLX is just because we have performed 
two Fast-Fourier transformations (FFTs) in each iteration cycle of FAG. Of 
course, in a parallel computing environment with distributed memory FFT 
routines usually do not scale well if the number of processors in increased. 
Therefore, for future studies we rather recommend to use the multigrid im- 
plementation of FAG as proposed by CUCCHIERI and Mendes [CM98]. 



A. 2. 2 A way to preselect best gauge copies 



We have seen above that fix- 
ing gauge configurations to Lan- 
dau gauge might become a CPU 
time intensive task on large lat- 
tices. This was even more inten- 
sive for us when we followed our 
f c-bc-strategy (see Sec. 4.1.3). 
For this, several random gauge 
copies for each gauge configura- 
tion U were gauge-fixed. Then 
we have selected the first gauge 
copy and that with the largest (fi- 
nal) functional value (Eq. (3.14)) 
to study the influence of Gribov 
copies on different observables. 
Obviously, it would be quite help- 
ful if the final ranking of func- 
tional values were known, without actually doing all the necessary itera- 




8000 



Figure A. 2: The values of F and 8 (see 
Eq. (3.14) and (A.l)) and are shown as functions 
of the iteration number for two sample gauge 
copies of the same configurations (/? = 6.0, 24**). 
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tions. 

Starting just from a ran- 
dom gauge copy of U we can- 
not forecast the final func- 
tional value, because those 
iterative gauge-fixing algo- 
rithms behave like a chaotic 
system under change of ini- 
tial conditions. For ex- 
ample, if the order of go- 
ing through the lattice is 
changed, a different func- 
tional value might been 
reached. Also checking the 
ranking of functional values 
after a fixed number of itera- 
tions is not very reasonable, 
since iteration processes are 
quite different as is illus- 
trated in Fig. A. 2. There the functional values F and the maximal violation 
of transversality 9 (Eq. (A.l)) are shown at intermediate iteration steps for 
two sample gauge copies. Although both iteration processes started from a 
(different) random gauge copy of the same configuration, the ways to conver- 
gence differ substantial. For one process, labeled as copy 15, the functional 
F reveals more often irregular jumps to larger values during the iteration loop 
than for the other process. This happens always in conjunction with sudden 
increases of 0. However, inspecting Fig. A. 3 such non-monotonous behavior 
is seen in the majority only for larger values of 9. Thus it seems that for 9 
below some threshold the final ranking of functional values is already fixed. 

From our study we know the final and a list of intermediate functional 
values for each gauge-fixed copy of U, because we actually have gauge-fixed 
each copy until the gauge condition was reached. Therefore, we can check now 
if the final ranking of functional values is already visible at an intermediate 
iteration state, i.e. at fixed, but larger values of 9. 

To be specific, in Fig. A. 4 we show the probability of whether that copy 
with the largest functional value at an intermediate value of 9 is that with 
the largest functional value (best copy) after gauge-fixing has finished. Of 
course, for 9 below the convergence criteria this holds and for 9 > 1 (not 
fixed) it does not. 

Looking at Fig. A. 4 we easily see that the probability of having found 
those copies which will result in the largest functional values increases with 



32'' 13 = 6.0 




2000 4000 6000 8000 10000 12000 



Figure A. 3: The local maximum of violation of 
transversality, 6, is shown as a function of the num- 
ber of iterations for over-relaxation (rlx) and Fourier 
accelerated gauge-fixing (fag). The lattice size is 32^ 
and P = 6.0. 
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lowering 9. For 9 ~ 10~^ the probability is almost 100%, independent of (3 
and the lattice size. In other words, to select the best gauge-fixed copy for 
each configuration U , it is enough to gauge-fix several random copies of U 
until the value 9 = 10""^ has been reached and then continue the process only 
for that with the largest intermediate functional value among all. In this way 
copies are singled out which will definitely not reach a larger functional value 
compared to others. If we had used this strategy we would have saved about 
40% (57%) of the total number of iterations for the 24^ lattice at /3 = 5.8 
(6.2), which is quite a lot. 
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Figure A. 4: In this figure we show the probabiHty of having found (at intermediate values 
of 8) those copies which will result in the largest functional values after the gauge-fixing 
process has converged. 



A. 3 Speeding up the inversion of the FP 
operator 

For the solution of the linear system Mcf) = t/>^ with symmetric matrix M, the 
conjugate gradient (CG) algorithm is the method of choice. Its convergence 
rate depends on the condition number, the ratio of largest to lowest eigenvalue 
of M. When all U^^^ = 1 obviously the FP operator is minus the Laplacian 
A with a diagonal color substructure. Thus instead of solving M(j) = i/}^ one 
rather solves the transformed system 

[MA-i](A(^)=Ve 
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CG 




PCG 


speed up 


lattice 


iter 


CPU[sec] 


iter CPU [sec] 


iter CPU[sec] 


8^ 


1400 


3.7 


570 2.4 


60% 35% 


16'^ 


3900 


240 


1050 130 


73% 46% 


32^ 


9900 


13400 


2250 3900 


77% 71% 



Table A. 2: The average number of iterations and CPU time per processor (PE) using the 
CG and PCG algorithm to invert the FP operator are given for different lattice sizes. All 
inversions have been performed at (3 = 5.8 with source S^'^ cxp(i k-y) where k = (1, 0, 0, 0). 
To compare the different lattice sizes 4 PEs have always been used. 

In this way the condition number is reduced, however, the price to pay is 
one extra matrix multiplication by A^^ per iteration cycle. In terms of CPU 
time this should be more than compensated by the reduction of iterations. 
The pre-conditioned CG algorithm (PCG) can be described as follows: 



initialize: 








7(0) = 


(p(O)^^(O)) 


start do loop: 


A; = 0,1,... 










r ('^■+1) = 










(-(/C^+i) < exit do loop 




^ + 7(fe) P 




end do loop 





Here (■, ■) denotes the scalar product. 

To perform the additional matrix multiplication with A^^ we used two 
fast Fourier transformations J-', due to (— A)^^ = JT^^ q^^{k) T . The perfor- 
mance we achieved is presented in Table A. 2. We conclude that on larger 
lattice sizes the reduction of iterations is about 70-75%, while the resulting 
reduction of CPU time depends on the lattice size. This is because we are 
using the fast Fourier transformations in a parallel CPU environment. If the 
ratio of used processors to the lattice size is small (see e.g. the data for 32"^ 
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lattice at this table), almost the same reductions of CPU time as for the 
number of iterations is achieved. 

Further improvement may be achieved by using the multigrid Poisson 
solver to solve Az ^''^ = r^^\ This method is supposed to perform better 
on parallel machines. Perhaps a further improvement is possible by using 
as pre-conditioning matrix — — A"-*^ — A~''^Mi A~^ + . . . which is an ap- 
proximation of the FP operator M — —A -|- Mi to a given order [Zwa94] (see 
also [FN04a]). However, the larger the order, the more matrix multiplications 
per iteration cycle are required. This may reduce the overall performance. 
We have not checked so far which is the optimal order. 
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